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Preface 


At  the  request  of  the  Defense  Atomic  Support 
Agency,  Oak  Ridge  National  Laboratory  has  under¬ 
taken  the  preparation  of  a  handbook  to  aid  engineers 
charged  with  the  responsibility  of  designing  shields  to 
protect  military  equipment  and  personnel  in  the  vicin¬ 
ity  of  a  nuclear  weapons  burst.  This  document  con¬ 
stitutes  the  fifth  chapter  of  the  Handbook  issued  thus 
far,  the  earlier  chapteis  being  Chapter  3,  “Methods  for 
Calculating  Neutron  and  Gamma-Ray  Attenuation,” 
Chapter  4,  “Neutron  and  Gamma-Ray  Albedos,” 
Chapter  5,  “Methods  for  Calculating  Effects  of  Ducts, 
Access  Ways,  and  Holes  in  Shields,”  and  Chapter  6, 
“Methods  for  Predicting  Radiation  Fields  Produced  by 
Nuclear  Weapons.” 

At  the  time  Chapters  3,  4,  and  S  were  published  it 
was  planned  that  the  Handbook  would  eventually  be 
bound  in  two  volumes,  with  Volume  I  serving  as  a 
textbook  and  ready  reference  and  Volume  II  serving  as 
an  engineering  guide.  Volume  1  would  include  an 
introductory  chapter,  a  chapter  describing  the  radiation 
fields  produced  by  nuclear  detonations,  and  Chapteis  3 
through  5.  Volume  II  would  present  two  or  more 
chapters  describing  engineering  methods  based  on  the 
more  sophisticated  techniques  described  in  Volume  I. 
During  the  preparation  of  Chapter  6,  however,  the  plan 
to  combine  the  individual  chapters  into  volumes  was 
abandoned  -  at  least  for  the  present  -  the  primary 
reason  being  that  it  would  complicate  the  later  publica¬ 
tion  of  updated  chapteis.  We  also  changed  the  order  of 
the  chapteis,  with  the  result  that  the  chapter  originally 
planned  as  Chapter  2  was  published  as  Chapter  6.  One 
reason  for  shifting  the  chapter  numbers  was  to  improve 
the  continuity  of  the  Handbook  by  putting  Chapter  6, 
which  describes  the  various  sources  of  radiation 
produced  by  a  nuclear  weapons  burst  and  presents 
techniques  for  calculating  the  transport  of  the  radia¬ 
tions  from  the  point  of  burst  to  the  surface  of  the 
shield,  immediately  preceding  the  engineering  chapters. 
Another  reason  was  that  an  additional  chapter  on  basic 
shielding  concepts,  the  present  one,  had  been  prepared 
and  it  would  logically  follow  the  introductory  first 
chapter.  With  this  reordering  the  basic  information 
generally  applicable  to  all  neutron  and  gamma-ray 


shielding  problems  is  presented  in  Chapters  2  through  5, 
and  the  information  speciilcally  api^icable  to  nuclear 
weapons  radiation  shielding  problems  is  presented  in 
the  chapters  beginning  with  Chapter  6.  The  engineering 
methods  dealing  with  the  design  of  shields  for  protec¬ 
tion  against  initial  weapons  radiation  will  be  divided 
into  two  chapters;  Chapter  7,  entitled  “Engineering 
Method  for  Designing  Initial  Radiation  Shields  for 
Blast-Hardened  Underground  Structures,”  which  will  be 
published  very  shortly,  and  Chapter  8,  “Engineering 
Method  for  Designing  Initial  Radiation  Shields  for 
Above-Ground  Structures.”  (Like  Chapter  6,  both 
Chapters  7  and  8  will  be  classified.)  It  now  appears  that 
methods  for  calculating  dose  rates  in  structures  may  be 
included  as  an  additional  chapter  in  the  Handbook. 
Whether  or  not  the  Handbook  will  include  a  chapter 
dealing  with  shields  to  protect  against  faUout  radiation 
is  at  present  uncertain.  If  a  fallout  chapter  is  published, 
it  will  rely  heavily  on  a  method  already  avidlable  and 
will  be  integrated  in  the  Handbook  only  so  that  all 
aspects  of  the  radiation  shielding  problems  will  have 
been  covered. 

In  order  to  prepare  this  Handbook,  it  has  been 
necessary  for  Oak  Ridge  National  Laboratory  to  obtain 
the  assistance  of  several  consultants  and  subcontractors. 
For  this  chapter,  for  example,  Paul  N.  Stevens,  a 
consultant  from  the  University  of  Tennessee,  together 
with  H.  Clyde  Claiborne  of  the  Laboratory,  prepared 
the  first  draft  with  which  the  editors  worked.  Other 
chapters  similarly  represent  a  cooperative  effort  of 
ORNL  staff  members  and  those  of  other  organizations. 

As  is  always  the  case  for  handbooks,  the  authors  and 
editors  must  rely  heavily  on  the  reviews  of  others  as  an 
aid  in  the  development  of  the  various  chapters.  The  list 
of  individuals  who  have  contributed  in  this  manner  Iras 
already  grown  very  large,  and  it  would  be  almost 
impossible  to  acknowledge  each  person  here.  However, 
there  are  always  reviewers  who  we  feel  have  made  such 
significant  contributions  as  to  warrant  individual 
acknowledgement.  With  respect  to  this  chapter,  we 
particulaily  wish  to  acknowledge  the  help  given  by  D. 
K.  Trubey  of  the  Laboratory,  who  spent  considerable 
time  reviewing  and  commenting  on  the  entire  chapter. 
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Appreciation  is  also  expressed  to  Maj.  F.  A.  Verser 
and  to  Lt.  Cols.  Charles  D.  Daniel  and  William  A. 
Alfonte,  who  as  past  DASA  Shielding  Project  Officers 
handled  the  early  administration  of  the  contract  and 
assisted  in  establishing  the  scope  of  the  Handbook. 
The  work  they  began  is  currently  being  ably  performed 
by  Maj.  R.  W.  Enz. 


Finally,  we  wish  to  thank  Mrs.  Viiginia  M.  Hamrick, 
who  assisted  in  the  editing  of  the  chapter,  and  Mrs. 
Virginia  Glidewell,  who  typed  and  helped  proofread  the 
many  drafts  which  are  always  necessary  precursors  to 
such  a  publication. 
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2.0.  Introduction 


Weapons  radiation  shielding  analysis  usually  involves 
three  major  steps:  describing  the  field  of  radiation 
against  which  protection  is  required;  determining  the 
extent  to  which  shielding  materials  will  attenuate 
pailicles  from  that  field ;  and  converting  the  number  of 
particles  that  penetrate  the  shield  into  a  dose  that  can 
be  related  to  a  physical  effect  or  to  a  biological 
response.  This  chapter  describes  the  basic  concepts 
underlying  these  steps.  For  those  who  plan  to  perform 
detailed  shielding  calculations  by  the  methods  described 
in  Chapters  3  through  6,  a  thorough  knowledge  of  these 
concepts  is  essential.  Such  knowledge  is  of  less  im¬ 
portance  to  those  who  plan  to  design  shields  merely  by 
applying  the  engineering  data  provided  in  Chapter  7  and 
subsequent  chapters,  but  it  would  be  extremely  helpful 
in  all  cases. 

The  first  step  in  the  analysis  -  describing  the 
radiation  field  -  begins,  of  course,  with  an  examination 
of  the  data  available  on  the  design  and  burst  conditions 
of  the  weapon  to  be  considered.*  The  resulting  descrip¬ 
tion  can  vary  from  a  simple  statement  of  the  total 
number  of  particles  of  a  particular  type  at  the  “point” 
of  detonation  to  a  detailed  specification  of  the  energies 
of  the  particles  and  their  locations  and  directions  at  all 
times  over  some  specified  volume.  It  is  obvious  that  in  a 
practical  situation  the  consideration  of  volume 
introduces  some  type  of  shielding  material  (for  ex¬ 
ample,  the  elements  in  the  weapon  or  the  atmosphere), 
and  the  effect  that  these  materials  have  on  the  radiation 
field  must  be  determined.  As  a  result,  the  first  two  steps 
of  the  analysis  become  inseparable,  the  two  together 
comprising  a  radiation  transport  problem.  The  custom 
is  to  divide  the  transport  problem  into  two  parts,  one  in 
which  the  attenuation  by  the  atmosphere  is  calculated 
and  another  in  which  the  attenuation  by  some  barrier  (a 
building  or  a  specially  designed  shield)  is  determined. 

The  various  methods  that  have  been  developed  for 
attacking  the  transport  problem  are  described  in  Chap¬ 
ter  3  of  this  Handbook,  and  the  application  of  some  of 


*The  interrelation  of  these  para’  icters  is  diseussed  in  detail 
in  Chapter  6. 


the  methods  is  demonstrated  in  subsequent  chapters. 
The  present  chapter  is  devoted  primarily  to  describing 
the  quantities  required  as  input  to  these  methods  and  to 
discussing  the  types  and  behavior  of  the  particular 
radiations  to  which  the  methods  are  applied  in  weapons 
radiation  shielding.  The  chapter  also  discusses  how  the 
results  of  the  transport  calculations  can  be  related  to  a 
physical  or  a  biological  effect,  that  is,  how  the  third 
step  of  the  analysis  can  be  accomplished. 

Discussions  of  the  quantities  used  to  describe  radia¬ 
tion  fields  are  presented  in  Sections  2.1  and  2.2.  The 
most  detailed  descriptions  are  given  in  seven¬ 
dimensional  phase  space  and  are  required  in  very 
sophisticated  transport  calculations.  However,  many 
transport  methods  cannot  handle  such  detail  and 
because  of  this,  plus  the  fact  that  the  detailed  descrip¬ 
tions  are  so  difficult  to  obtain,  descriptions  based  on 
fewer  parameters  have  come  into  common  use.  Also, 
some  special  descriptions  of  radiation  fields  that  are 
basic  to  many  of  the  mathematical  models  have 
emerged  and  these  forms  will  often  suffice. 

The  extent  to  which  shielding  materials  will  attenuate 
particles  from  a  radiation  field  depends  both  on  the 
type  and  energy  of  the  radiation  and  on  the  kind  of 
materials  comprising  the  shield.  Some  particles  will  pass 
through  the  shield  unimpeded,  and  others  will  collide 
with  the  nuclei  of  the  medium.  Particle-nucleus  colli¬ 
sions  can  result  in  a  change  in  the  direction  and  energy 
of  the  incident  particle  or  in  its  absorption,  the 
attenuation  within  the  shield  being  determined  by  the 
probabilities  for  the  various  types  of  interactions.  The 
quantities  representing  the  probabilities,  called  cross 
sections,  are  discussed  in  Section  2.3. 

The  formal  definitions  and  explanations  given  in 
Sections  2. 1-2.3  are  presented  in  terms  of  idealized 
particles  and  can  represent  any  type  of  neutral  radia¬ 
tion.  Insofar  as  shielding  against  nuclear  weapons  is 
concerned,  the  only  radiations  that  need  be  considered 
are  neutrons  and  gamma  rays.  While  both  these  radia¬ 
tions  can  be  regarded  as  neutral  particles,  there  are 
certain  obvious  differences  in  their  behavior.  For 
example,  the  speed  of  a  gamma-ray  “particle”  is 
constant  and  does  not  depend  on  its  energy,  whereas 
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the  speed  of  a  neutron  is  determined  by  its  energy. 
Thus  some  caution  must  be  exercised  when  specializing 
particle  concepts  to  gamma  rays. 

Neutrons  and  gamma  rays  can  undergo  many  different 
interactions  with  materials,  and  calculations  of  their 
transport  through  a  shield  can  be  quite  difficult.  The 
situation  is  complicated  by  the  fact  that  the  interaction 
of  a  neutron  can  produce  a  gamma  ray.  The  converse  is 
also  possible,  but  is  much  less  likely.  The  types  of 
interactions  that  neutrons  and  gamma  rays  can  undergo 
are  discussed  in  Sections  2.5  and  2.7  respectively,  and 
the  interactions  that  give  birth  to  these  radiations  are 
discussed  in  Sections  2.6  and  2.8  respectively. 

Not  all  quantities  that  can  be  determined  by  transport 
calculations  (for  example,  differential  flux)  are  easily 
measured  or  easily  related  to  a  physical  or  biological 
effect.  However,  a  quantity  that  has  come  to  be  known 
as  dose  and  that  is  determined  from  the  energy 
deposited  by  the  radiation  in  the  receiver  (for  example, 
in  an  instrument  or  in  a  human  body)  is  measurable, 
and  calculated  results  can  be  related  to  the  same 
quantity.  As  a  result,  whenever  a  physical  or  a 


biological  effect  is  being  considered,  dose  is  the 
quantity  that  is  used  to  specify  the  design  limit  on  the 
intensity  of  the  radiation  field. 

Unfortunately,  several  concepts  of  dose  exist,  and 
considerable  confusion  has  arisen  with  respect  to  their 
definitions.  Because  the  problem  has  not  yet  been  fully 
resolved,  all  viewpoints  are  presented  in  Section  2.9  of 
this  chapter  and  the  definitions  that  are  employed  in 
this  Handbook  are  specified.  Section  2.10  discusses 
allowable  dose  limits  and  the  probable  effects  of  such 
acute  doses  as  might  be  encountered  in  a  nuclear 
accident  or  in  warfare. 

Finally,  it  is  pointed  out  that  although  this  chapter 
was  prepared  primarily  to  describe  basic  concepts  and 
particle  interactions  of  especial  interest  to  the  designer 
of  weapons  radiation  shields,  its  contents  will  be 
equally  useful  to  those  engaged  in  reactor  shield  design 
since  the  radiations  involved  are  the  same.  The  informa¬ 
tion  will  also  be  useful  to  those  who  design  radiation 
shields  for  accelerators  or  spacecraft,  but  for  these 
purposes  it  is  incomplete  since  many  more  different 
types  of  radiations  must  be  considered. 


2.1.  Quantities  Used  to  Describe  Radiation  Fields 


A  complete  description  of  a  radiation  field  would 
specify  the  number  of  particles  of  a  given  type  and 
energy  which  at  a  given  time  exist  at  a  given  position 
and  travel  in  a  given  direction,  with  all  particle  types, 
energies,  positions,  and  directions  being  considered  at 
all  times.  However,  such  complete  descriptions  are 
seldom  required  -  except  in  the  most  sophisticated 
particle  transport  calculations  -  and  a  radiation  field 
cannot  be  measured  in  such  detail.  Therefore  particle 
populations  are  usually  described  in  terms  of  a  fewer 
number  of  parameters.  The  quantities  used  most 
frequently  in  shielding  analysis  are  listed  in  Table  2.1 
and  discussed  in  the  following  paragraphs.  The  names 
given  to  the  differential  quantities  have  not  been 
standardized,  and  occasionally  one  finds  that  a  term 
specifically  identified  as  being  differential  is  actually 
less  differential  than  a  companion  term  not  so 
designated.  While  in  the  first  column  of  Table  2.1  an 
attempt  is  made  to  specify  the  degree  of  differentiation 
in  each  case,  the  names  used  in  the  text  are  those  in 
keeping  with  general  usage. 


2.1.1.  PARTICLE  DENSITIES 

A  knowledge  of  tlie  particle  density  over  all  phase 
space  is  equivalent  to  a  complete  solution  of  a  particle 
transport  problem  and  comprises  more  information 
than  is  available  from  most  calculational  schemes  now 
in  use.  When  given  in  seven-dimensional  phase  space, 
which  consists  of  three  spatial  coordinates,  the  par¬ 
ticle’s  energy,  two  direction-defining  angles,  and  time, 
this  particle  density  is  defined  by 

n(r,E,^,t)dE  jn  =  the  number  of  particles  per  unit 
volume  at  space  point  f  and  time  l 
having  energies  in  dE  about  energy 
E  and  directions  in  t/H  about  the 
unit  direi  tion  vector  f2. 

The  particle  density  so  defined  is  doubly  differential  - 
in  energy  and  in  direction  -  and  less  detailed  forms  will 
often  suffice.  For  example,  one  may  use  the  steady- 


state*  particle  density  differential  in  energy  only, 
commonly  referred  to  as  differential  particle  density 
and  defined  by 

n(r,E)  dE  =  the  number  of  particles  per  unit  volume  at 
space  point  r  having  energies  in  dE  about  E 

=  f  n(r,E,n)dndE .  (2.1) 

J4jr 

Or  one  may  use  the  steady-state  total  particle  density, 
defined  as  the  number  of  particles  per  unit  volume  at 
space  point  r  and  given  by 

=  f  f"ni?,E,ti)  dE  da  =  r  n{},E)  dE .  (2.2) 
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2.1.2.  FLUX  DENSITIES 

Even  though  the  concept  of  particle  density  is 
basically  simple  and  has  a  unique  interpretation,  experi¬ 
ence  has  shown  that  the  Pux  density^  or,  as  it  is  more 
commonly  referred  to,  the  flux,^  serves  better  as  the 
dependent  variable  in  solutions  of  the  transport  equa¬ 
tion  (refer  to  Chapter  3,  Section  3.1).  The  flux  is 
related  to  the  particle  density  through  the  particle’s 
speed,  and  when  described  in  terms  of  seven¬ 
dimensional  phase  space  the  flux  is  given  by 

‘i>(r,£', a.t)  =  v  n{r,E.a. t)  ,  (2.3) 

where  v  is  the  particle’s  speed  and  corresponds  to  the 
energy  E.  (The  speed  is  the  scalar  magnitude  of  the 
particle’s  velocity  vector  P.) 


*The  stcady-slate,  or  time-independent,  condition  will  be 
denoted  in  phase  space  notation  by  dropping  the  time  symbol  /; 
for  example,  seven-dimensional  phase  space  becomes  six  di¬ 
mensional  and  in  phase  space  notation  is  given  by  (r.f.n). 

^Although  this  quantity  is  truly  a  density  and  the  Inter¬ 
national  Commission  on  Radiation  Units  and  Measurements 
recommends  the  use  of  the  term  flux  density,  the  simpler  term 
flux  is  ingrained  in  shielding  terminology  and  will  be  used  in 
this  and  other  chapters  of  this  Handbook. 
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Table  2.1.  Quantities  Used  to  Describe  Particle  Populations 
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The  flux  defined  in  the  preceding  paragraph  is  doubly 
differential  and  is  usually  referred  to  as  the  angular 
flux.  Greater  insight  into  the  use  of  the  angular  flux  as 
the  dependent  variable  in  mathematical  descriptions  of 
particle  transport  is  provided  by  its  interpretation  either 
as  the  tracklengths  traversed  per  unit  volume  and  time 
or  as  the  flow  of  particles  per  unit  area  and  time. 

The  tracklength  interpretation  of  angular  flux  follows 
from  the  observation  that  the  speed  of  an  individual 
particle  can  be  considered  as  its  scalar  tracklength  per 
unit  time.  The  product  of  particle  density  and  speed  is 
then  the  sum  of  the  tracklengths  traced  by  all  the 
particles  within  a  unit  volume  per  unit  time,*  in  which 
case  the  definition  of  the  angular  flux  would  be 

dt'dSl  =  the  total  tracklengths  traversed  per 
unit  volume  and  time  at  space 
point  r  and  time  t  by  particles 
having  energies  in  dE  about  energy 
E  and  directions  in  dO.  about  n. 

The  interpretation  of  the  angular  flux  as  a  flow  of 
particles  per  unit  area  and  time  is  closely  related  to  the 
concept  of  angular  current  (see  discussion  of  angular 
current  in  Section  2.1.4).  It  will  be  shown  in  the 
discussion  on  current  that  the  angular  flux  is  identical 
with  the  magnitude  of  the  current  vector  /  and  thus  can 
be  interpreted  as 

<tir,E,il,l)  dE  dn  =  the  flow  per  unit  area  and  time  at 
space  point  r  and  time  t  of  particles 
having  energies  in  dE  about  E  and 
directions  in  about  f2. 

When  the  transport  and  deposition  of  the  particle’s 
kinetic  energy  are  of  concern,  the  energy  flux  differ¬ 
ential  in  energy  and  angle  is  often  used.  This  quantity, 
referred  to  as  the  angular  energy  flux,  is  defined  by 

l{r,E,il.t)dEdn  =  ihe  energy  flow  per  unit  area  and 
time  at  space  point  Fand  time  t  due 
to  particles  having  energies  in  dE 
about  E  and  directions  in  dn  about 

n. 


♦Note  that  the  per-unit-time  units  of  the  flux  are  associated 
with  the  particle's  speed,  which  is  a  function  only  of  the  energy. 
However,  the  time  dependence  of  the  flux  is  a  consequence  of 
the  time  behavior  of  the  particle  density,  which  does  not  have 
time  units  even  though  the  time  symbol  t  is  included  in  the 
phase  space  notation  to  denote  a  dependence  on  time;  that  is, 
nif.E.ii.t). 


and  is  related  to  the  angular  flux  by 

l(r.E,n,t)  =  E^r,E.n.t).  (2.4) 

Many  calculational  models  employ  less  detailed  de¬ 
scriptions  of  the  flux;  for  example,  they  may  use  a 
description  in  which  steady  state  is  assumed  and  the 
angular  dependencies  are  removed  by  the  appropriate 
solid-angle  integration.  This  reduces  the  description  to 
four-dimensionab  phase  space  and  results  in  a  dependent 
variable  representing  the  particle  flux  differential  in 
energy.  Commonly  known  as  the  differential  flux,  this 
quantity  is  given  by 

<t>(rE,n)dn.  (2.5) 

•»4r 


Like  the  angular  flux,  the  differential  flux  can  be 
interpreted  in  terms  of  tracklength  per  unit  volume  and 
time  or  in  terms  of  the  number  of  particles  that  enter  a 
unit  sphere  per  unit  time.  In  the  latter  case,  the 
solid-angle  integration  can  be  regarded  as  a  summing  of 
particles  that  enter  a  sphere  of  unit  cross  section, 
regardless  of  their  directions  of  motion.  The  sphere  is  in 
effect  generated  by  the  rotation  a  circular  unit  area 
during  the  integration  over  47r  solid  angle  (see  Fig.  2.1). 
In  this  context  the  definition  for  the  differential  flux 
can  be  restated  as 

•PCr.E’)  dE  =  number  of  particles  having  energies  in  dE 
about  E  which  enter  a  sphere  of  unit  cross 
section  per  unit  time  at  space  point  F. 

While  this  definition  of  the  differential  flux  is  descrip¬ 
tive,  it  is  not  exact.  The  mathematically  rigorous 
definition  is 


<P{r.E)dE  = 


lim  N{E)dE 
A4  -♦■0  AA 


(2.6) 


which  implies  the  limit  process  A/1  -♦  0,  with  N(E)  dE 
denoting  the  number  of  particles  having  energies  in  dE 
about  E  which  enter  an  “incremental”  sphere  of  cross 
section  AA  per  unit  time. 

The  concept  of  the  incremental  sphere  is  the  best  way 
to  visualize  the  energy  flux  differential  in  energy. 
Referred  to  as  the  differential  energy  flux,  this  quantity 
may  be  defined  by 

I{r,E)  d£  =  the  energy  flow  per  unit  time  into  a  sphere 
of  unit  cross  section  at  space  point  r  due  to 
particles  having  energies  in  dE  about  E. 


UNCLASSIFIED 
ORNU-  DWG  69-I4306R 
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A/l 


(PROJECTED  AREA) 


Fig.  2.1.  Incremental  Sphere  Concept  of  Flux. 

The  differential  energy  flux  is  given  by 

lCr,E)  dE=  f  I(r.E,n)  dndE  =  E  dE .  (2.7) 

J4ir 

Other  quantities  used  are  the  total  flioc,  the  total 
energy  flux,  the  group  flux,  ihe  group  angular  flux,  and 
the  discrete  ordinates  flux.  The  total  flux,  defined  al¬ 
ternatively  as  the  total  particle  tracklength  per  unit 
volume  and  time  at  space  point  r  or  as  the  number  of 
particles  that  enter  a  sphere  of  unit  cross  section  per 
unit  time  at  space  point  7,  is  obtained  by  integrating 
the  differential  flux  over  all  energies: 

<P(r)=  f’‘M.f,E)dE .  (2.8) 

’'0 

Similarly,  the  total  energy  flux,  defined  as  the  total 
energy  flow  per  unit  time  into  a  sphere  of  unit  cross 
section  at  space  point  f,  is  obtained  by  an  integration  of 
the  differential  energy  flux  over  all  energies: 

I{r)=  I{f.E)dE  .  (2.9) 

•fo 


The  total  flux  has  only  limited  application  to  practi¬ 
cal  shielding  problems  because  of  the  strongly  energy- 
dependent  nature  of  the  particle’s  behavior.  A  more 
useful  approach  is  to  divide  the  total  energy  range  into 
NOG  energy  intervals,  called  energy  groups, 


and  to  define  the  group  flux  as  the  integral  of  the 
differential  flux  over  the  corresponding  energy  group: 

<|>g(r)  =  ^r.E)  dE,  (2.10) 

I 


with  the  obvious  constraint  that 

=  (2.11) 

NOG 

The  group  angular  flux  (group  tlux  differential  in 
angle)  has  a  similar  definition  and  is  obtained  by 
integrating  the  angular  tlux  over  a  specific  energy 
group: 

4>c(r.n)=  «l>(r./f,n)di'.  (2.12) 

’'Egt  1 


For  problems  that  involve  directional  symmetry,  the 
group  angular  flux  can  be  rewritten  in  terms  of  a  new 
angular  variable  p  -  f2-r/|r|,  the  direction  cosine: 

‘l>(;(f,p)  dp  =  <P(j  (r.n)  r/n  .  (2.13) 


The  group  angular  flux  can  then  be  defined  as  the  total 
particle  tracklength  per  unit  volume  and  time  at  space 
point  F  of  particles  with  energies  within  energy  group 
AEq  and  directions  defined  by  direction  cosines  which 
lie  in  dp  about  p. 

A  mathematical  treatment  which  evaluates  the  angu¬ 
lar  flux  at  discrete  angles,  called  discrete  ordinates 
forms  the  bases  for  many  detailed  numerical  solutions 
of  the  transport  equation.  A  form  of  discrete  ordinates 
could  be  simply  group-angular-flux  values  that  corres¬ 
pond  to  given  spatial  positions  and  particle  directions. 
However,  a  more  useful  form  is  obtained  by  integrating 
the  group  angular  flux  over  specified  segments  of  p 
space, Apo=Pj  D=  1,2  ..  .NOA, 


AEq  =Eg- Eg^^,*  G=  1,2  . . .  NOG, 


g.d(^=  ^ 


(2.14) 


*The  subscripts  g  and  fi-l  refer  to  the  upper  and  lower  limits 
respectively  of  the  Cth  energy  group,  and  G=  1  corresponds  to 
the  highest  energy  group.  An  alternate  convention  would  as¬ 
sociate  G=1  with  the  lowest  energy  group  and  the  subscripts 
g4|  and  g  would  then  correspond  to  the  upper  and  lower 
energy  limits. 


The  "mean  value”  approximation  to  Eq.  2.14  is  given 
by 
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where  is  a  mean  value*  for  the  direction  cosfne  over 
the  increment  The  discrete  ordinate  fluxes 

are  those  values  of  ^(f)  which  correspond 
to  the  spatial  meshpoints  r,: 

^l,G,D  •  (2.16) 


2.1.3.  FLUENCES 

A  useful  measure  of  total  exposure  to  a  flux  for 
applications  involving  energy  deposition  or  dose  is  the 
integral  quantity  called  fluence.  Fluence  is  defined  by 
the  International  Commission  on  Radiation  Units  and 
Measurements  as  the  quotient  of  hN  divided  by  AA, 
where  tsN  is  the  number  of  particles  that  enter  a  sphere 
of  cross-sectional  area  A4  and  the  A’s  imply  a  limiting 
process.'  This  definition  is  equivalent  to  regarding 
fluence  as  a  time-integrated  flux  over  some  specified 
time  interval.  As  such,  the  fluence  can  be  written  as 

=  ri 

^t)dt,  (2.17) 

I 

where  Ar  corresponds  to  some  specified  time  interval 
and  '!>(/)  can  be  any  one  of  the  several  kinds  of  fi  x 
described  in  Section  2.1.2.  For  example,  the  energy 
fluence  is 

F{f)=  l{f.t)dt .  (2.18) 

•''I 

(NOTE:  The  symbolism  for  particle  fluences  and 
fluxes  has  been  confusing  in  the  past,  with  the  symbols 
nvt  and  <!>  being  used  for  both.  In  order  to  avoid 
confusion  in  this  chapter.  <l>  is  reserved  for  particle 
fluxes  and  when  particle  fluences  are  intended,  the 
symbol  is  used.) 


*tip  may  be  adjusted  to  give  the  equality:  for  well-behaved 
functions  the  closer  tip  is  to  the  true  mean,  the  better  the 
approximation. 

tThc  following  subscript  notation  is  used  throughout  this 
section  and  in  Chapter  3:  subscripts  /.  G,  and  D  denote 
functions  whose  values  are  associated  with  the  /th  space  inter¬ 
val,  the  f/th  energy  group,  and  the£)th  angular  interval,  respec¬ 
tively;  /  and  i-i-l  refer  to  a  function  evaluated  at  the  lower  and 
upper  limits  of  the  /th  space  interval;  gs-i  and  g  refer  to  a 
function  evaluated  at  the  lower  and  upper  limits  of  the  Gth 
energy  group;  d*\  and  d  refer  to  a  function  evaluated  at  the 
lower  and  upper  limits  of  the  //th  angular  interval. 


2.1.4.  CURRENT  DENSITIES 

The  characteristic  property  of  the  current  variable,  is 
its  close  relation  to  the  convective  (leakage)  effects  in 
the  theoretical  description  of  particle  transport.  The 
most  general  form  of  the  current  variable  is  differential 
both  in  energy  (or  speed)  and  in  angle.  Called  the 
angular  current  density  or,  more  frequently  the  angular 
current,  this  quantity  is  symbolized  by7(rv,n,/)and  is 
defiled  as  the  directed  flow  per  unit  a:ea  (normal  to 
the  n  direction)  and  time  at  the  space  point  rand  time 
/  of  particles  having  speeds  in  dv  about  v  and  directions 
in  dSl  about  SI. 

The  relation  between  current  and  particle  density  can 
be  established  by  considering  that  (1)  the  product  of  v 
and  the  particle  density  can  be  regarded  as  a  vector  sum 
of  the  individual  co-directional  velocity  vectors  (v), 
yielding  the  resultant  vector  J  =nv,OT  (2)  the  (y  X  dt 
X  dA  X  n)  particles  contained  within  the  volume 
element  shown  in  Fig.  2.2  will  all  exit  through  the 
deferential  area  dA  within  differential  time  dt  if  v  = 
vf2.  These  models  can  be  expressed  mathematically  in 
terms  of  the  angular  current  as 

J(r.v,Sl,t)  dvdSl  =  Slv  n(r,v,Sl,t)  dv  dSl ,  (2.19) 

and  since  the  particle’s  kinetic  energy  is  a  function  of 
its  speed,  Eq.  2.19  can  be  rewritten  as 

J{r.E.Sl,t)  dE dVi  =  Vi  V  n(f,E.Vi,t)  dE dH  ,  (2.20) 

where  J{r,E, SI, t)  dE  =  J(r,v,Sl,t)  dv  and  n(f,E,Vi,t)  dE  = 
n{r,v,Sl,t)  dv.  Then  when  v  n{r,E,Sl,t)  is  identified  as 
the  angular  flux,  the  relation  between  the  angular  flux 
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Fig.  2.2.  Particle  Flow  Concept  of  Current.  The  particles 
contained  within  the  volume  ejement  will  all  exit  through  the 
differential  area  dA  within  dt'dv  =  vn. 


and  the  angular  current  pointed  out  in  Section  2.1.2  is 
obtained: 

j(r,E,  n,t)dEdn  =  n  ^(r.E.  n,r)dEdn.  (2.2 1 ) 

Other  less  detailed  forms  of  the  angular  current  are 
the  group  at^lar  current  jQ{r, SI)  and  the  total  angular 
current  J(r,Sl),  which  are  obtained  by  integrations  over 
an  energy  group  and  all  energies,  respectively; 

Jcir.Sl)  J{r.E.Sl)  dE  =  fi  ;  (2  22) 

Vi 

7(r,n)  =  / ”7(r.£.n)  dE  =  n  <t>(r.Si)  •  (2.23) 

•^0 

It  is  apparent  that  the  angular  current  variables  all  have 
essentially  the  same  simple  relation  with  the  corres¬ 
ponding  angular  flux  variables,  because  the  energy 
integrations  are  performed  directly  on  the  flux  vari¬ 
ables.  For  example,  in  the  case  of  the  group  angula: 
current, 

Tmr.E.n)dE 

V' 

=  ^r,E.n)  dE  =  n  •P(j(r.n) .  (2,24) 

Vi 


The  integral  of  the  angular  current  over  all  directions 
(47r  solid  angle)  constitutes  a  vector  summation,  and  the 
resultant  vector  is  regarded  as  the  net  current,  often 
called  simply  current. 

The  net  current  differential  in  energy  only,  referred 
to  as  the  differential  net  current,  is  defined  by 


Jir.E.t)  dE  =  lhe  net  flow  per  unit  area  and  time  at 
space  point  r  and  time  t  of  particles  having 
energies  in  dE  about  E,  where  the  unit 
area  is  normal  to  the  direction  of  the 
resultant  vector  J{r,E,t) 


=  r  J{r.E,Sl,t) 

JAn 


dSldE 


which  is  the  angular  flux-weighted  vector  summation  of 
the  unit  vectors  SI  over  47r  solid  angle.  The  group  net 
current  JQif)  and  total  net  current  J{r)  are  4Tr  solid- 


angle  integrations  of  the  group  angular  current  and  total 
angular  current,  respectively: 

4(0  =  x.  y^(r,n)  dSl  =J^  SI  4'G(r,f2)  dSl  ;  (2.26) 
7(r)=r  7ir,n)dSl=  f  n<iir.Sl)dn.  (2.27) 

Jin  J^it 

The  flow  of  n-directed  particles  across  an  arbitrarily 
oriented  differential  area  is  a  necessary  concept  in  the 
description  of  the  directed  flow  of  particles  in  terms  of 
a  specific  coordinate  system  and  can  be  related  to  the 
angular  current  J{Sl)*  by  consideration  of  Fig.  2.3, 
where  the  direction  vector  n  is  normal  to  the  differ¬ 
ential  area.  The  number  of  H-directed  particles  crossing 
the  differential  area  dA  per  unit  time  is  equal  to 
J(Si)  (n  dA).  A  scalar  current  i„(f2)  that  describes  the 
flow  of  the  r2-directed  particles  per  unit  area  normal  to 
the  direction  n  is  defined  as 

J„in)dA  -Ji^yindA) :  (2.28) 

it  follows  that 

j„(a) = n-7(n)  =  n-«  *(0) = coi  6  «t(n) ,  (2.29) 

where 

n  =  the  unit  vector  corresponding  to  an  arbitrary 
direction, 

n  =  a  coordinate-identifying  subscript,  for  ex¬ 
ample,  n  =  x  when  n  =  /, 

J„(S2)  =  the  flow  of  H-directed  particles  per  unit  area 
(normal  to  the  direction  n)  and  time, 

<l>(n)  =  the  angular  flux  v^iable  corresponding  to  the 
angular  current  J(Sl). 

Note  that  y„(S2)  is  a  scalar  quantity  but  is  uniquely 
related  to  the  coordinate  system  through  the  direction 
unit  vector  n.  The  corresponding  vector  current  is  given 
by 

7„(n)  =  n y„(n)  =  n  7{n)  h .  (2.30) 

wtwre  the  vector  J„{n)  is  the  component  of  the  vector 
J(Sl)  with  respect  to  the  n  direction  and  J„{Sl)  is  the 

*ln  this  further  discussion  of  current,  the  notation  7(n)  will 
be  used  to  denote  any  of  the  angular  currents  and  J  the 
corresponding  net  current;  for  example,  7(n)  may  re,sresent 
7(r.£n,f),  JQlr.n.t),  or  7(f.n)  and  7  may  represent  j(r.E,t), 
J(f.r),  or7(r). 


t 
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Fig.  2.3.  Schematic  of  Particle  Flow  Acroaa  an  Arbitrarily 
Oriented  Surface. 


where 

J*  =  the  flow  per  unit  area  (normal  to  the  n  direc¬ 
tion)  and  time  of  particles  having  positive 
direction  cosines  relative  to  the  n  direction 

J~  =  the  flow  per  unit  area  (normal  to  the  «  direc¬ 
tion)  and  time  of  particles  having  negative 
direction  cosines  relative  to  the  n  direction 

The  partial  currents  do  not  occur  explicitly  in  most 
theoretical  formulations  of  particle  transport; however, 
the  boundary  conditions  of  many  problems  of  interest 
are  best  expressed  in  terms  of  the  partial  current 
variable. 


projection  of  the  vector  7(12)  on  the  n  axis  and  also  the 
magnitude  of  the  component  7„(12). 

The  three  components  of  7(J1)  in  cartesian  coordi¬ 
nates  are  given  by 

7,(12)  =  77^(12) =7«lKn)(7n), 


7/12)  =77^(12)  =/<l>(12)  (712), 

7j(12)  =  k  7^(12)  =  k  4>(12)  (it  H) . 

(2.31) 

and  7(12)  is  equal  to  the  vector  sum  of 
components: 

the  three 

7(12)  =  7,(0)  +  7^(0)  +  7j(0) 

=  77,(0) +  77v(0)  +  fc  7,(0). 

(2.32) 

The  net  flow  of  particles  across  a  unit  area  normal  to 
the  h  direction  is  obtained  by  integrating  7„(0)  over  all 
directions: 

J„=  f  7,,(0)  70  =  r  n-7(0)  70  =  n  J 

.  (2.33) 

Further,  7„  can  be  represented  in  terms  of 
currents*  J*  and  7"  such  that 

partial 

•’n  ''n  '’n  ’ 

(2.34) 

2.1. 5.  ADJOINT  FLUX 

The  adjoint  flux  does  not  correspond  in  the  usual 
sense  to  a  real  physical  quantity.  Instead  it  is  the  flux 
associated  with  the  adjoint  equation  when  that  equa¬ 
tion  is  interpreted  to'  describe  the  transport  of  hypo¬ 
thetical  radiation  particles,  sometimes  called  “adjunc- 
tons.”  The  rationale  of  the  concept  of  adjunctons 
follows  at  least  in  part  from  the  similarity  of  tikC  ?djoint 
equation  to  the  Boltzmann  transport  equation.t  Inte¬ 
gral  forms  of  these  two  equations  are  given  by  Eqs.  2.35 
and  2.36  respectively,  and  a  coupling  between  the  two 
is  provided  by  Eq.  2.37: 


<P*(P)  =  f  K{P  -*  P')  <t>*(P’)  dP'  +  f{P) , 

(2.35) 

X(/’)  =  / X(^')  Kip'  ^P)dP'^  S{P) , 

(2.36) 

X  =  fs{P)  <!>*(/>)  TP  =/ xiP)  fiP)  dP  , 

(2.37) 

where 

<1>*(P)  =  the  adjoint  flux,  defined  as  the  ex¬ 
pected  value  of  the  present  and  future 
contributions  to  the  effect  of  interest 

by  a  radiation  particle  leaving  a  col¬ 
lision  at  a  point  whose  coordinates  in 
phase  space  are  denoted  by  P, 


♦The  partial  currents  are  usually  defined  to  be  positive 
numbers,  with  their  difference  (J*  -  7^  determining  the  sign  of 
7„. 


'I  Derivations,  relations,  and  other  useful  and  interesting 
properties  of  the  various  forms  of  the  Boltzmann  equation  are 
given  in  Chapter  3  of  this  Handbook. 
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K{P'  -*  P)  dP  =  the  probability  that  a  particle  experi¬ 
encing  an  event  at  a  point  P  in  phase 
space  will  experience  another  event  at 
a  point  whose  phase  space  coordinates 
lie  within  dP  about  P, 

/^  =  the  payoff  function,  which  is  the 
probable,  or  expected,  direct  contri¬ 
bution  to  the  effect  of  interest  by  a 
particle  leaving  a  collision  at  point  P  in 
phase  space, 

X{P)  dP  =  the  expected  number  of  particles  leav¬ 
ing  collisions  at  points  whose  phase 
space  coordinates  lie  within  dP  about 
P,  \(P)  being  called  the  “emergent 
particle  density,” 

S(P)dP  =  the  expected  number  of  initial 
(source)  events  which  occur  within  the 
phase  space  coordinates  dP  about  P, 

X  =  the  effect  of  interest  which  can  be 
calculated  if  either  (the  adjoint  flux 
and  source  function)  or  (the  emergent 
particle  density  and  payoff  function) 
are  known . 

The  adjoint  equation,  which  mathematically  defines 
the  characteristics  of  the  adjunctons,  has  often  been 
referred  to  as  the  “backward”  equation  because  the 
movement  of  adjunctons  is  backward  to  that  of  real 
particles.  That  is,  they  originate  at  a  point  in  phase 
space  associated  with  the  detector  and  move  in  the 
direction  of  phase  space  that  corresponds  to  the  real 
source.  This  behavior  requires  that  in  a  collision  the 
adjuncton  gain  energy,  which  contradicts  the  usual 
energy  conservation  principles  unless  prior  to  the 
collision  the  target  nucleus  has  the  kinetic  energy  it 


would  have  acquired  in  a  real  collision  process  and  the 
particle  has  the  final  energy  of  the  adjuncton.  This 
precludes  the  assignment  of  physical  properties  to  the 
adjunctons  that  would  resemble  those  of  real  particles. 

The  lack  of  physical  properties  is  not  a  serious 
shortcoming,  however,  since  adjunctons  are  not  useful 
as  radiation  particles  per  se.  Rather,  the  adjoint  flux 
which  may  be  mathematically  evaluated  for  some  point 
in  phase  space  has  the  precise  interpretation  of  “value 
function.”  The  value  function  has  been  defined  to  be  a 
measure  of  the  present  and  future  contributions  of  a 
real  particle  at  some  point  in  phase  space  to  an  effect  of 
interest  and  as  such  has  been  used  to  bias  the  Monte 
Carlo  analysis  of  the  behavior  of  real  particles  (neutrons 
and  gamma  rays).  This  is  accomplished  by  modifying 
the  basic  stochastic  processes  so  that  event  densities 
that  contribute  most  strongly  to  the  effect  of  interest 
are  sampled  most  frequently.  The  transmission  proba¬ 
bility  used  in  the  Monte  Carlo  analysis  is  modified  by 
using  the  value  function  as  an  importance  function, 
with  the  particle’s  “weight”  being  adjusted  to  preserve 
parity  with  the  true  analog  calculation. 

Calculation  of  the  answer  of  interest  through  Monte 
Carlo  solution  of  the  adjoint  problem  avoids  the 
difficulties  associated  with  estimation  to  a  point  detec¬ 
tor,  provided,  of  course,  that  the  source  has  a  more 
agreeable  configuration.  Also,  the  adjoint  concept 
provides  a  fundamental  basis  for  perturbation  analyses. 

Since  the  adjoint  equation  is  very  similar  to  the  usual 
Boltzmann  transport  equation,  it  is  not  surprising  that 
methods  developed  to  solve  the  latter  will  also  apply  to 
the  former.  For  example,  adjoint  fluxes  have  been 
calculated  using  the  techniques  of  diffusion  theory, 
Monte  Carlo,  and  ransport  theory.  Consequently,  there 
should  be  no  serious  practical  or  fundamental  diffi¬ 
culties  in  calculating  the  adjoint  flux. 


2.2.  Special  Radiation  Fields 


Certain  radiation  fields  are  basic  to  many  of  the 
mathematical  models  that  describe  particle  transport, 
and  they  possess  functional  simplicities  that  sometimes 
permit  analytical  solutions.  In  this  section  four  such 
flelds  are  presented  in  terms  of  the  various  flux  and 
current  variables.  They  are  a  monodirectional  beam,  an 
isotropic  flux  field,  a  point  source,  and  an  isotropic 
plane  source. 

2.2.1.  MONODIRECTIONAL  BEAM 

The  monodirectional  beam,  shown  schematically  in 
Fig.  2.4,  consists  of  a  one-dimensional  system  of 
particles  that  are  all  moving  in  the  same  direction  Hq.  It 
is  described  in  most  detail  in  terms  of  the  angular 
current: 

1^,E, fi)  dEdfL^fl  Jo  F{E)  6(n,no)  dE dn  ,  (2.38) 

where  Jq  is  the  magnitude  of  the  curr^t  (particles 
cm"^  sec“‘)  for  a  unit  area  normal  to  the  f2o  direction, 
F{E)  is  the  normalized  energy  spectrum  with  the 
property 

rF{E)dE=\,  (2.39) 

0 

and  5(n,f2o)  is  an  angular  delta  function*  with  the 
property 
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/  g(n)6(n,f2o)rfn=g(f2o),  (2.40) 

471 

a  mathematical  operation  that  aligns  all  direction 
vectors  with  f2o  • 

When  the  monodirectional  beam  is  specified  to  be 
monoenergetic,  which  it  frequently  is  for  calculational 
convenience,  the  Dirac  delta  function  with  the  property 

/^”  m  6(£  -  £o)  dE^HEo)  (2.41) 

may  be  used  to  describe  the  energy  “spectrum”  of  the 
beam.  That  is, 

J(r,E,fi)dEdfi 

*  n  /o  6(£’  -  Eo)  6(?2,no)  dE  dVt .  (2.42) 

Integration  of  Eq.  2.42  over  all  directions  (47r  solid 
angle)  yields  the  differential  net  current 

j(r,E)=  f  Jij,E,n)dfL 

477 

=JoF{_E)  f  n6(n,no)<ff2 

4»r 

=  fio  Jo  E{E)  ,  (2.43) 

and  a  subsequent  integration  over  all  energies  gives  the 
total  net  current 

=  7(7,  E)dE 

=  fio  Jo  S’  EXE)  dE  =  noJo.  (2.44) 


Fig.  2.4.  Monodirectional  Beam  Incident  on  the  Shield. 
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*A  delta  function  is  simply  a  mathematical  notation  to  indi¬ 
cate  that  the  function  is  non-zero  at  only  one  particular 
value  of  the  argument. 
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The  angular  flux*  for  the  beam  is  equal  to  'he  scalar 
magnitude  of  the  angular  current  and  is  given  by 

<I)(r>.S2)  =  /o  SCf  -  Ao)  6(n.i2o)  .  (2.45) 

The  differential  flux  is  obtained  by  integrating  the 
angular  flux  over  47r  steradians, 


and  the  total  flux  is  the  differential  flux  integrated  over 
all  energies, 

‘l>o(7)  =  /"‘l>o(7.A0t/A\  (2.50) 

The  angular  current  is 


d<r.£)  =  /  <t>(r.E.n)dn 

4jr 

=  /oS(£-£o).  (2.46) 

and  the  total  flux  is  obtained  by  integrating  the 
differential  flux  over  all  energies, 


j(7,A',n)  =  n  <Kr,A',n)  =  ,  (2.51 ) 

47r 

and  the  differential  net  current  is  obtained  by  inte¬ 
grating  the  angular  current  over  all  directions.  For  an 
isotropic  flux  field  the  differential  net  current  is 
characteristically  zero: 


<lir)=f“<P(r.F)c/F  =  yo. 


(2.47) 


J(r.£)=  J(r,£.ii)dn 


The  beam  concept  is  primarily  used  to  specify  a 
source  configuration  at  a  surface  boundary  and  as  such 
is  usually  regarded  as  a  boundary  condition.  The  partial 
current  with  respect  to  the  positive  x  direction  is 

j;  =  (Hoi)  Jo  .  (2.48) 

The  partial  current  in  the  negative  direction, cannot 
be  used  as  a  boundary  condition  in  this  situation 
because  the  reflected  particles  would  also  contribute  to 
the  partial  current  in  the  negative  direction,  which  is 
unknown  prior  to  the  solution  of  a  problem  with  an 
incident  beam  source  configuration. 


2.2.2.  ISOTROPIC  FLUX  FIELD 

The  radiation  fleld  for  many  physical  situations  can 
be  regarded  as  being  isotropic  with  respect  to  the 
angular  distribution  of  the  particle  velocity  vectors;  that 
is,  motion  in  any  direction  is  equally  probable.  The 
angular  flux  for  this  condition  is  simply  the  differential 
flux  4>o(r,£)  divided  by  4n  steradians. 


*o{r.£) 

4ir 


(2.49) 


*The  flux,  usually  referred  to  as  the  “incident  flux,"  would 
not  be  the  flux  existing  in  real  situation  since  backscattering 
would  contribute  to  the  flux  at  a  surface. 


=  /  f2Jf2  =  0.  (2.52) 

4jr  4Tr 

Under  the  concept  of  an  isotropic  flux  field  the 
positive  partial  cunent  is  sometimes  used  to  specify  an 
equivalent  surface  source  configuration.  Consider,  for 
example,  a  unit  area  normal  to  the  direction  vector  h  as 
shown  in  Fig.  2.5,  The  scalar  angular  cuirent  that 
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Fig.  2.5.  Orientation  of  n-Directed  Particles  with  Respect  to 
a  Unit  Aica  Nonnal  to  the  Direction  Vector  E 
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describes  the  flow  of  f2-directed  particles  in  terms  of  a 
unit  area  normal  to  the  n  direction  is  given  by 

J„(r,EM)  =  n-j(r,ESi)  =  «  «  .fl) .  (2.53) 


2.2.3.  POINT  SOURCE 

The  general  volumetric  source  term  in  seven-dimen¬ 
sional  phase  space  is  defined  as 


and  the  positive  partial  current  that  conesponds  to  the 
n  direction  is  obtained  by  appropriate  integrations: 

n'n^Q  0 


=  /  _  j”  n-n^^^clEdO. 


47r 


Arr 


cos  0(27r  sin  d  dd) 


(2.54) 


It  can  be  seen  that  the  total  flux  <I>o(r),  which 
corresponds  to  a  47r  solid-angle  integration  of  the 
angular  flux,  is  exactly  four  times  as  large  as  the 
corresponding  positive  partial  current. 

In  many  real  situations  the  angular  flux,  although 
anisotropic,  can  be  regarded  as  being  isotropic  in  the 
forward  half-space  (or  in  27r  angular  directions),  ’^or 
this  case  the  positive  partial  cuitent  would  have 
essentially  the  same  magnitude  as  a  positive  partial 
current  in  an  isotropic  flux  field  would  have.  As  a  result 
it  has  become  standard  practice  in  shielding  calculations 
to  define  a  new  total  flux,  referred  to  as  the  27r  total 
flux*  and  given  by 


S(r,E,Sl,t)drdEdn  =  the  number  of  particles  having 
energies  in  dE  about  ^and  direc¬ 
tions  in  dn  about  S2  that  are 
emitted  per  unit  time  at  the  time 
t  from  the  differential  volume  df 
at  the  space  point  r. 

A  useful  specialization  of  the  general  source  term  is  the 
angular  point  source,  which  is  given  by 

SQ(ro,E,El,t)dE  dn  =  X\\e  number  of  particles  having 
energies  in  df  about  E  and  direc¬ 
tions  in  dn  about  SI  that  are 
emitted  per  unit  time  at  time  t 
from  a  point  source  located  at 
the  space  point  To . 

The  angular  point  source  is  related  to  the  general  source 
term  by  the  identity 


So(ro,E,Sl,t) 

=  f  f  S(r,E.n',t)  5(n,n')  dr  dSl’ ,  (2.57) 
•'4rr  ■'all  space 

where  the  general  source  term  is  specialized  for  the 
point-source  configuration  by  inclusion  of  the  Dirac 
delta  function 


=  f  <P(r.S2)dS2.  (2.55) 

”  2n 

where  the  27r  solid-angle  integration  is  accomplished  in 
the  forward  half-space.  Since  for  the  isotropic  flux  field 
the  4n  total  flux  is  exactly  two  times  the  2ir  total  flux, 
the  positive  partial  current  can  be  expressed  as 

y>)  =  «l'^(7)/2.  (2.56) 

From  the  above  it  is  obvious  that  care  must  be 
exercised  in  applying  this  formula,  which  was  derived 
on  the  basis  of  a  J^r)  source  condition,  since^ither  the 
2ff  total  flux  4>J''(r)  or  the  47r  total  flux  <I>j,(r)may  be 
specified  for  a  typical  problem.  Th£  conespoi^ng 
partial  current  would  be  either  4>§”(r)/2  or  4>Q(r)/4, 
and  a  factor-of-2  error  could  result  unless  it  is  clear 
which  total  flux  is  specified. 


*The  2n  total  flux  is  also  called  the  incident  flux. 


Sir,E,Sl'.t)  =  S(r  -  ro)So(E,Sl'.r) ,  (2.58) 

and  So(E,Sl',t)  is  an  energy-,  angle-,  and  time-dependent 
source  intensity  factor.  For  the  usual  problem  the 
assumption  of  energy-angle  separability  is  permissible, 
and  the  source  intensity  factor  can  be  expressed  as 

So(E,n.l)  =  So{l)F(E)g(Sl ) ,  (2.59) 

where  So{t),  the  source  intensity  factor,  is 

=  f  f  So(E,si.t)dEdn,  (2.60) 

4ir  •'O 

and  the  other  terms  are  defined  as 

5(r  -  To)  =  the  Dirac  delta  function  which  locates  the 
point  source  at  To , 
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/•’(f)  =  normalized  energy  spectrum, 

=  normalized  directional  probability  density 
function  having  the  property ^"4^  dil 
=  1, 

g(j2)c/S2  =  fraction  of  particles  having  directions  in 
dil  about  S2, 

5(n,f2')  =  angular  delta  function  which  aligns  the 
directions  of  the  emitted  particles  with 
the  direction  S2, 

f2  =  direction  vector  which  lies  along  the 
vector  7  7o  as  shown  in  Fig.  2.6. 

Evaluation  of  the  integrals  in  Eq.  2.S1  shows  that 
So{?oJ:.n,t)  =  So{t)F{t^g{Tl)  ■  (2.61) 


tiiat  the  various  mathematical  descriptions  of  the  point 
source  presented  here  apply  to  both  steady-state  and 
lime-dependent  source  conditions.  For  the  steady-state 
problem  the  phase  space  notation  becomes  (r,2;',S2).  and 
an  emission  rate  constant  in  time  is  presumed. 

Many  source  configurations  can  be  represented  by  a 
point  source  or  a  simple  arrangement  of  point  sources; 
for  example,  isotropic  line  and  plane  sources  consist  of 
appropriate  distributions  of  isotropic  point  sources. 
Also,  the  kernel  technique,  known  as  the  method  of 
Green’s  functions,  has  for  its  basis  the  point  kernel,  the 
point  kernel  being  the  solution  to  a  problem  with  a  unit 
point  source.  In  summary,  the  point  source  is  a  most 
elementary  and  versatile  source  configuration  and  has 
considerable  application  to  complex  source  configura¬ 
tions  such  as  the  detonation  of  a  nuclear  weapon. 


For  isotropic  emission  ^^2)  dil  =  dill4n,  in  which  case 

So{roM.n.!)  =  ^  F(E) .  (2.62) 

47r 

The  differential  point  source  is  obtained  by  integrating 
Eq.  2.61  over  all  directions.  Thus 

Soi7o.t\t)  =  Soit)F(H)  .  (2.6.1) 

Similarly,  the  total  number  of  particles  emitted  per  unit 
time  at  time  l.  5o(2),  is  obtained  by  integrating  the 
differential  point  source  ever  all  energies.  Finally,  if 
total  emission  of  particles  is  desired,  the  time- 
dependent  source  intensity  factor  Sgit)  is  integrated 
over  some  specific  interval  of  time.  It  should  be  noted 


2.2.4.  ISOTROPIC  PLANE  SOURCE 

An  fsoiropic  plane  source,  depicted  in  Fig.  2.7.  can  be 
regarded  as  an  arbitrary  distribution  of  unit  isotropic 
point  sources,  and  in  six-dimensional  phase  space  is 
defined  as 

.S(.voxV.r./i'.i2)c//'.’t/12  =  the  number  of  particles  having 
energies  in  dl-^  about  /•.’  and 
directions  in  dil  about  S2  that 
are  emitted  per  unit  time  from 
a  unit  area  located  at  .vo.,i',:  and 
normal  to  the  direction  vector  fi 

=  S(y.:)-^^dF,/n  .  (2.64) 

4;r 
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Fig.  2.6.  Point  Source  Specified  in  General  r  Space. 
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Fig.  2.7.  Urientition  of  Isotropic  Plane  Source  with  Respect 
to  the  Rectangular  Coordinate  System, 


17 


where 

.Vo  =  location  of  the  source  plane  along  tiie  x  axis, 

5(_v.:)  =  planar  distribution  of  unit  isotropic  point 
sources. 


— —  =  unit  angular  point  source  energy  distribution 
for  isotropic  emission. 

For  an  isotropic  plane  source  having  a  uniform 
distribution  of  point  sources, 50’,r)  =  S  (particles  cm"* 
sec ),  and  the  expression  for  the  source  becomes 

S{xo.t:,h)JEdh=S— ciEdH  .  (2.65) 

47r 


azimuthal  symmetry;  that  is, 


<J)(xo.£.f2)  =  (n  «)‘'  S(xo.E.n) 


4;r  (n-n) 


(2.66) 


The  utility  of  ^Kvo.f.n)  is  sometimes  limited  by  the 
singularity  experienced  as  (il-n)  approaches  zero.  The 
term  S(jro.£'.n)  has  proved  to  be  the  more  useful  form, 
and  integrations  over  direction  fl  and  energy  E  lead  to 
physically  meaningful  quantities. 

The  differential  form  for  the  plane  source  is  obtained 
by  integrating  Eq.  2.64  over  all  directions. 


Tliis  source  is  of  considerable  importance  in  shielding 
work;  for  example,  it  is  used  to  represent  radiation 
fields  such  as  those  due  to  radioactive  fallout. 

Note  that  S(xo,E,ii}  is  based  on  a  unit  area  that  is 
normal  to  the  .v  axis.  The  angular  flux  that  corresponds 
to  the  isotropic  plane  source  has  for  its  basis  a  unit  area 
normaNo  the  fl  direction  and  therefore  differs  from 
S(xo.E,Tl)  by  a  factor  of  (flvi)  '  or  cos  0  with 


S{xo.E)=  f  S(xoE.n)dn=SE(E)  ,  (2.67) 

and  the  total  emission  (particles  cm"*  sec"')  is 
obtained  by  the  subsequent  integration  over  all  en¬ 
ergies, 

5(xo)  =  /  S{xo,E)dE  =  So. 

0 


(2.68) 


2.3.  Quantities  Used  to  Describe  Radiation  Interactions 


In  describing  the  transport  of  radiation  through 
matter,  the  interactions  that  the  various  particles  may 
have  with  the  nuclei  of  the  medium  must  be  considered. 
The  probability  of  the  occurrence  of  a  given  inter¬ 
action,  which  will  vary  with  the  type  of  particle,  its 
energy,  and  the  material  it  traverses,  is  usually  given  in 
terms  of  microscopic  or  macroscopic  cross  sections. 
These  cross  sections  are  described  in  the  following 
paragraphs,  along  with  their  relation  to  reaction  proba¬ 
bilities  and  reaction  rates. 

2.3.1.  MICROSCOPIC  CROSS  SECTIONS 

The  microscopic  cross  section,  designated  by  the 
symbol  o,  can  be  considered  to  be  the  effective  area  that 
a  nucleus  presents  to  an  oncoming  particle,  regardless  of 
the  particle  direction.  That  is,  a  particle  will  interact 
with  a  nucleus  with  a  probability  of  unity  when  the 
particle  is  incident  on  the  area  and  with  a  probability  of 
zero  when  it  misses  the  area.  The  units  usually 
associated  with  a  are  cm^/atom  or  barns/atom,  where  I 
bam*  =  10'^^  cm’.  Because  of  the  wave-like  behavior 
of  the  particles  and  the  complex  nature  of  the  nucleus, 
the  magnitude  of  a  in  general  bears  no  direct  relation 
with  the  physical  size  of  the  nucleus,  but  is  strictly  a 
function  of  the  properties  of  the  atom  and  the  energy 
and  type  of  the  interacting  particle. 

2.3.2.  REACTION  PROBABILITIES 

If  a  nucleus  having  a  cross  section  o  is  contained 
within  some  characteristic  volume  of  material, 
which  has  a  projected  area/l^,  then  the  probability  that 
a  particle  incident  on  the  volume  will  interact  with  that 
particular  nucleus  is  the  ratio  of  the  effective  area  o  to 
the  projected  area  (see  Fig.  2.8fl).  It  follows  then 
that  ojAg  is  the  reaction  probability  per  nucleus 
associated  with  the  area  A^  and  that  it  can  be 
interpreted  as  the  average  number  of  reactions  oc¬ 
curring  in  Vf,  per  incident  particle. 


*The  term  barn  is  attributed  to  have  resulted  from  a  humor¬ 
ous  remark  to  the  effect  that  10'^*  cm*  looked  as  big  as  a  barn 
to  a  neutron. 


If  the  nuclide  density  N  (nuclides/cm*)  is  uniform 
throughout  the  total  number  of  nuclei  contained 
within  is  given  by  NV^  and  the  reaction  probability 
associated  with  the  total  number  of  nuclides  is  equal  to 
the  linear  sum  of  the  reaction  probabilities  associated 
with  the  individual  nuclei  provided  that  the  character¬ 
istic  volume  is  optically  thin  (no  self  shielding).  The 
reaction  probability  now  associated  with  the  character¬ 
istic  volume  is  NV^{olA^). 

This  summation  of  probabilities  is  shown  schemat¬ 
ically  in  Fig.  2.8/).  The  optically  thin  condition  is 
assured  if  the  characteristic  volume  is  differential  in 
size.  (The  characteristic  area  A^  would  also  be  differ¬ 
ential  in  size.)  Otherwise,  the  condition  shown  in  Fig. 
2.8c  might  result,  in  which  case  Uie  effective  cross 
section  per  nucleus  would  be  decreased  because  of  the 
self-shielding  effect.  This  follows  since  the  reaction 
probability  would  no  longer  be  a  simple  linear  sum¬ 
mation  of  the  individual  probabilities.  This  can  easily  be 
seen  for  the  example  of  an  opaque  volume;  increasing 
the  thickness  has  no  further  effect. 


2.3.3.  MACROSCOPIC  CROSS  SECTIONS 

The  reaction  probability  associated  with  the  optically 
thin  characteristic  volume  can  be  rearranged  and 
written  as  {aNy,  VjA^).  If  VJA^.  is  recognized  as  the 
characteristic  differential  distance  dS^,  traveled  by 
particles  within  and  if  oN  is  identified  as  the 
macroscopic  cross  section  then  2k/5p  is  the  reaction 
probability  associated  with  the  differential  tracklength 
dS^.  It  follows  that  a  formal  definition  for  the 
macroscopic  cross  section  £  is  the  reaction  probability 
per  differential  tracklength,  or,  presuming  that  differ¬ 
ential  tracklength  is  implicit  in  £  and  interpreting 
reaction  probability  as  reactions  per  incident  particle, 
the  definition  for  £  becomes  the  number  of  reactions 
per  particle  tracklength.  However,  the  units  for  £  are 
usually  given  as  reciprocal  centimeters  (cm"' ). 


^The  symbol  ^  is  used  for  gamma  rays. 
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There  is  another  approach  in  defining  the  macro¬ 
scopic  cross  section.  In  particle  transport  problems  large 
numbers  of  particles  are  involved,  and  so  it  is  sometimes 
convenient  to  describe  the  process  in  terms  of  averages 

UNCLASSIFIED 
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CHARACTERISTIC  VOLUME  (  cm^  ) 
HAVING  A  PROJECTED  AREA  (cm2) 


O  -  (MICROSCOPIC 
^  CROSS-SECTION) 


Cr 


'PARTICLE 


(<7) 


or  expected  values.  It  can  be  assumed  that  the  particle 
has  a  certain  probability  of  undergoing  an  interaction 
with  a  nuclide  within  a  distance  dx.  This  average  or 
expected  value,  called  the  macroscopic  cross  section 
£  in  units  of  inverse  length,  depends  on  the  particle 
energy  and  the  density  and  nuclear  properties  of  the 
medium.  Thus  in  penetrating  a  distance  dx  in  a 
medium,  a  particle  will  suffer  Sr/jc  average  number  of 
interactions  per  unit  time. 

In  order  to  examine  the  concept  more  closely, 
consider  the  simple  model  of  a  beam  of  particles 
normally  incident  on  an  infinite  slab  of  some  material. 
Let  <l>(x)  be  the  uncollided  flux  at  thickness  x.  The 
number  being  removed  in  dx  is 

-r/«lKx)  =  ‘l>(^)2Jx,  (2.69) 


which  integrates  to 

4>(x)  =  <l>oe--'' .  (2.70) 

where  «l>o  is  the  incident  flux.  The  probability  that  a 
particle  will  suffer  an  interaction  in  a  distance  less  than 
L  is 


,  _,.zi 

Jo  <l>0  •'0 

(2.71) 

where  is  the  probability  density  function  for  an 

interaction.  The  mean  distance  a  particle  will  travel  (the 
mean-free  path  X)  before  it  suffers  an  interaction  is 

\=f"xle--^dx=^.  (2.72) 

•'o  2, 


Fig.  2.8,  Relationship  Between  Microscopic  Cross  Sections 
and  Reaction  Probabilities. 

(a)  A  nucleus  having  a  microscopic  cross  section  a  is 
contained  within  a  characteristic  volume  of  material  which 
has  a  projected  area  A  particle  is  entering  and  the 
probability  that  it  will  interact  with  the  nucleus  is  given  by  the 
ratio  olAf..  This  ratio  can  be  interpreted  as  the  number  of 
reactions  per  incident  particle  per  nucleus. 

(/>)  Several  nuclei,  each  having  a  microscopic  cross  section  o 
are  contained  within  the  volume  V^,  which  is  sufficiently  thin 
that  the  o's  will  not  overlap.  The  probability  that  the  incident 
particle  will  interact  with  a  nucleus  in  is  given  by 
(olAc)NVf.,  where  N  is  the  number  of  nuclides  per  cm^  in 
(c)  The  characteristic  volume  is  sufficiently  thick  for  the 
o’s  to  overlap.  In  this  case  the  self-shielding  effect  would 
prevent  the  reaction  probability  from  being  a  simple  linear 
summation  of  the  individual  probabilities. 


If  the  foregoing  discussion  is  accepted  as  establishing 
the  definition  of  a  macroscopic  cross  section,  which  is  a 
function  of  both  the  properties  of  the  nuclide  (nucleus 
only  for  the  neutron)  and  the  particle  energy  and  type, 
the  microscopic  cross  section  a  can  be  defined  as  simply 
£/A',  where  N  is  the  nuclear  density  and  a  represents 
the  cross-sectional  area  per  nucleus  seen  by  the  incident 
particles. 

It  should  be  clear  on  the  basis  of  the  foregoing  dis¬ 
cussion  that  £  should  not  be  used  as  a  probability 
density  function  to  determine  reaction  probabilities 
except  in  special  cases.  Even  though  the  reaction  proba¬ 
bility  for  a  particle  normally  incident  on  a  differential 
volume  dx  within  the  slab  is  given  by  Srfx,  the 
probability  that  a  reaction  will  occur  within  the 
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optically  thick  slab  is  not  simply  /  ^J.x  =  As  was 

■'o 

shown  above,  the  reaction  probability  for  the  slab  is 
given  by  (1  -  e  which  approaches  2/.  only  for 
the  special  case  of  an  optically  thin  slab  since  ( 1 
e  forSZ.  <\* 

2.3.4.  REACTION  RATES 

The  concept  of  a  macroscopic  cross  section  as 
described  in  the  preceding  section  suggests  that  particles 
which  move  through  matter  will  undergo  interactions  at 
a  rate  proportional  to  the  total  particle  tracklength,_a 
quantity  that  is  closely  related  to  the  particle  tlux  d>(r) 
(see  Table  2.1).  The  macroscopic  cross  section  S  is  in 
fact  the  required  proportionality  constant  between  the 
reaction  rate  and  the  particle  flux;  the  reaction  rate  per 
unit  volume  at  the  space  point  r  is  given  by 

/?(7)  =  2«J>(r) ,  (2.73) 


iXA')  =  macroscopic  cross  section  for  the  reac¬ 
tion,  which  depends  on  the  particle’s 
energy  as  well  as  the  nuclear  character¬ 
istics  of  the  medium, 

=  differential  particle  flux  density. 

The  reaction  rate  fur  an  energy  interval,  or  group 
{AEq  =  Eg  Eg+  1 ),  is  simply  related  to  the  differential 
reaction  rate  by 

t  1 

=  /^7  m:')<iir.tr)JE .  (2.75) 

1 

This  leads  to  a  definition  of  the  reaction  rate  in  terms 
of  an  energy-averaged  macroscopic  cross  section 
(group  cross  section)  and  tlie  group  flux; 


the  dimensional  correctness  of  which  is  verified  by 
noting  that 

reactions  _  reactions  particle-cm 
cm^-sec  particle-cm  cm^-sec 

The  differential  reaction  rate  is  related  to  the  differ¬ 
ential  particle  flux  and  is  given  by 


where 


I _ 

/'■*  >l>{f,E)JE 
1 


R{f,E)  =  £(£1  Mf,E) ,  (2.74)  and 


(2.76) 


(2.77) 


where 

R(FE^dE  =  number  of  reactions  at  space  point  f  per 
unit  volume  and  time  which  involve 
particles  having  energies  in  dE  about  the 
energy  E, 


*Another  and  more  definitive  description  of  the  opticaliy 
thin  condition. 


=  'Kr£)tZ£'.  (2.78) 

^■g*  I 

The  concept  of  the  reaction  rate  in  its  many  forms 
is  fundamental  to  the  formulation  of  the  various 
equations  that  analytically  describe  particle  transport 
(for  example,  the  Boltzmann  transport  equation)  and  to 
the  determination  of  useful  information,  such  as  energy 
deposition  and  biological  dose  from  known  or  calcu¬ 
lated  values  of  the  particle  flux  density. 


I 


2.4.  Radiations  Produced  in  Nuclear  Weapons 


As  pointed  out  in  the  introduction  to  this  chapter, 
several  types  of  radiations  are  produced  by  nuclear 
detonations,  but  generally  the  only  radiations  that  the 
shield  designer  must  be  concerned  with  are  neutrons 
and  gamma  rays.  It  is  the  nature  and  classification  of 
these  radiations  that  are  discussed  in  the  following 
paragraphs.  The  processes  whereby  neutrons  and 
gamma  rays  interact  with  materials  are  discussed  in 
Sections  2.5  and  2.7  respectively,  and  the  processes 
whereby  they  are  produced  are  described  in  Sections 
2.6  and  2.8  respectively. 


Outside  a  nucleus  the  neutron  is  not  a  stable  particle 
and  decays  with  a  half  life  of  12.8  min.  However,  this 
effect  is  always  neglected  in  neutron  transport  problems 
because  of  the  extremely  short  transport  lifetimes  (in 
general  much  less  than  1  sec). 

Ine  theory  of  wave  mechanics  suggests  that  all 
moving  particles,  including  neutrons,  have  wave  proper¬ 
ties.  The  deBroglie  wavelength  h/p  can  be  assigned  to  a 
neutron,  and  if  its  kinetic  energy  is  less  than  the  energy 
equivalent  of  its  rest  mass,+  that  is,  if  E  <939.5  MeV, 
the  nonrelativistic  wavelength  is  given  by 


2.4.1.  NEUTRONS 

The  presence  of  a  liighly  penetrating  electrically 
neutral  radiation  had  already  been  experimentally 
demonstrated*  by  several  groups  of  investigators  when 
Chadwick^  in  1932  established  the  existence  of  a 
neutral  particle  having  a  mass  nearly  equal  to  that  of  a 
proton.  Since  then,  this  neutral  particle,  which  has 
become  known  as  the  neutron,  has  proved  to  be  not 
only  a  key  to  nuclear  weaponry  but  also  a  bane  to  those 
persons  charged  with  providing  protection  against 
nuclear  radiation  fields. 

Neutrons,  like  protons,  are  constituents  of  the 
nucleus,  and  both  are  referred  to  as  nucleons.  The 
ability  of  a  neutron  to  transfonn  into  a  proton  -  and 
the  converse  -  suggests  that  protons  and  neutrons  may 
be  simply  different  states  of  a  single  particle,  and  a  full 
understanding  of  the  relation  between  the  two  is  one  of 
the  current  problems  in  nuclear  physics.  Neutrons  can 
be  regarded  as  fundamental  particles  since  they  are 
readily  identifiable  by  the  following  physical  properties: 


Charge 
Rest  mass 
Spin 

Magnetic  dipole  moment 

Statistics 

Intrinsic  parity 


0 

I.0UH98  atomic  mass  units 

Vr 

1.913  nuclear  magnetons 
I'enni 
I'vcn 


*This  radiation  was  produced  when  beryllium  was  bombarded 
by  alpha  particles. 


2.8600  X  10*” 
1,112 


(2.79) 


where 

X  =  deBroglie  wavelength  (cm), 

E  ~  neutron  energy  (MeV), 

/i  =  Planck’s  constant  (6.624  X  10*^’  erg-sec,  or 
4.132  X  10*'*  eV-sec), 

p  =  momentum  of  the  neutron. 

The  neutron’s  wavelength  is  inversely  proportional  to 
its  velocity  and  provides  valuable  insist  into  the 
general  nature  of  neutron  interactions  with  matter. 
Themial  neutrons  have  wavelengths  of  the  same  order 
as  the  atomic  spacings  in  crystals.  Therefore,  when 
thermal  neutrons  interact  with  a  crystal,  diffraction 
patterns  in  accordance  with  Bragg’s  law  are  produced. 
(The  main  difference  between  the  diffraction  of  neu¬ 
trons  and  of  gamma  rays  is  that  neutrons  interact  with 
the  nuclei  whereas  gamma  rays  interact  primarily  with 
the  orbital  electrons.)  Neutrons  with  energies  from  1  to 
100  MeV  have  wavelengths  of  the  order  of  nuclear 
dimensions,  and  neutrons  with  relativistic  energies 
greater  than  10  GeV  have  wavelengths  comparable  with 


^The  conversion  factor  that  relates  a  particle's  rest  mass  in 
atomic  mass  units  (amu)  to  the  energy  equivalent  of  its  rest 
mass  is  931  MeV/amu.  The  mass  of  a  neutron  is  1,00898  amu, 
which  when  multiplied  by  931  MeV/amu  yields  939.5  MeV  for 
its  energy  equivalent. 


21 


22 


the  spacing  of  nucleons  within  the  nucleus,  which 
provides  another  tool  for  exploring  the  structure  of  the 
nucleus.* 

Neutrons  are  classified  according  to  energy,  with  the 
divisions  occurring  somewhat  naturally  because  of 
behavior  peculiar  to  particular  energy  ranges.  In  some 
cases  the  divisions  are  somewhat  arbitrary,  and  they  can 
be  ambiguous  since  the  energy  structures  overlap  and 
differences  in  opinion  exist  as  to  the  proper  termi¬ 
nology.  However,  classification  is  helpful  in  discussing 
neutron  behavior.  Some  of  the  classifications  used  are 
given  below. 

Slow  Neutrons.  Neutrons  with  energies  from  zero  to 
about  I  or  2  keV  are  sometimes  referred  to  as  slow 
neutrons.  This  is  a  broad  classification,  however,  and  is 
not  in  general  use  since  it  does  not  serve  as  a  very 
precise  description  of  the  energy  state  of  low-energy 
neutrons.  In  reactor  physics  and  technology  jargon, 
the  term  “slow  neutrons”  refers  to  thermal  neutrons. 

Thermal  Neutrons.  Thermal  neutrons  are  in  thermal 
equilibrium  with  the  atoms  (or  molecules)  of  the 
medium  in  which  they  are  present.  A  particular  themial 
neutron  may  gain  or  lose  energy  in  any  one  collision 
with  the  nuclei  of  the  medium,  but  there  is  no  net 
energy  change  for  a  large  number  of  neutrons  diffusing 
in  a  nonabsorbing  medium.  For  this  condition  the 
kinetic  energy  of  the  neutrons  is  distributed  statistically 
according  to  the  Maxwell-Boltzmann  equation  as  de¬ 
rived  from  the  kinetic  theory  of  gases.  The  equation  for 
the  spectral  distribution  is 


2m 


e-EjkT  i;\l2  ^ 


(2.80) 


where 

dn  =  n{E)  dE  -  number  of  neutrons  having  energies  in 
the  energy  interval  dE  about  E, 

n  =  total  number  of  neutrons  in  the  system, 

E  »  neutron  energy  (eV). 

A:  =  Boltzmann  constant  (8.61  X  10"’  eV  per^K), 

T  =  absolute  temperature  (°K). 

It  can  be  shown  that,  when  the  Maxwell-Boltzmann 
distribution  applies,  the  energy  corresponding  to  the 
most  probable  neutron  velocity  (2200  m/sec  at  20° C)  is 
given  by  kT  (0.0252  eV),  the  average  kinetic  enerr,y  of 


*Additional  details  on  the  properties  of  neutrons  are  available 
in  most  texts  on  nuclear  theory  (see,  for  example,  a  book  by 
Evans^). 


themial  neutrons  is  %kT,  and  tlic  most  probable 
neutron  energy  is  %  kT. 

For  absorbing  media  the  true  energy  distribution  of 
the  neutrons  deviates  from  the  Maxwell-Boltzmann 
distribution  because  neutron  absorptions  tend  to  in¬ 
crease  with  decreasing  energy,  which  causes  an  increase 
in  the  average  energy.^  This  tendency  is  illustrated  in 
Fig.  2.9.  The  spectral  distribution  shown  in  Fig.  2.9a 
exhibits  very  little  distortion,  which  is  characteristic  of 
poorly  absorbing  media,  while  the  distribution  shown  in 
Fig.  2.96  clearly  demonstrates  the  effects  of  the 
preferential  absorptions  at  the  lower  neutron  energies. 
This  phenomenon  is  called  hardening  and  can  in  part  be 
compensated  for  by  the  use  of  an  “effective  neutron 
temperature”  T,,,  which  is  defined  as  a  number  that, 
when  used  in  the  Maxwell-Boltzmann  distribution,  gives 
the  best  least-squares  fit  to  the  tme  distribution  as 
determined  by  a  Monte  Carlo  calculation.  The  effective 
neutron  temperature  T„  has  been  found  to  be  approxi¬ 
mately  related  to  tlie  environmental  temperature  in 
the  following  way:'* 


T  =  T 
‘  n  '  m 


1  + 


1.84  2:^(A:r) 


=  7',,,  11+0.46 A]  ,  (2.81) 


where 

^a(kT)  ~  macroscopic  absorption  cross  section  evalu¬ 
ated  at  the  energy  kT, 

=  macroscopic  scattering  cross  section  for 
thermal  neutrons, 

5  =  average  logarithmic  energy  decrement  per 
collision, 

42^ 

• 

Cold  Neutrons.  -  Neutrons  with  an  average  energy 
less  than  that  of  thermal  neutrons  are  sometimes 
refened  to  as  cold  neutrons.  They  are  of  no  concern  in 
shielding  calculations. 

Epithermal  Neutrons.  -  Epithemial  neutrons  are  the 
neutrons  in  a  system  whose  energies  exceed  those  of  a 
Maxwell-Boltzmann  distribution  for  the  effective  tem¬ 
perature  of  the  neutron-supporting  medium. 

Resonance  Neutrons.  -  Neutrons  having  energies  in 
the  range  between  1  and  100  eV  are  frequently  called 
resonance  neutrons.  Many  nuclei  exhibit  strong  absorp¬ 
tion  of  neutrons  at  well-defined  eneigies  in  this  energy 
range,  and  these  absorptions  are  commonly  referred  to 
as  resonance  absorptions. 
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UNCLASSIFIED 
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Fig.  2.9,  Comparison  of  the  Neutron  Energy  Disbributions  with  Maxwell-Boltzmann  Distribution  in  (a)  Poorly  Absorbing 
Medium  and  (A)  Highly  Absorbing  Medium.  (From  Coveyou  et  al,  ref.  4.) 


Intermediate  Neutrons.  -  Neutrons  in  the  energy 
range  1  keV  to  about  0.5  MeV  are  often  referred  to  as 
intennediate  neutrons.  The  slowing  down  and  removal 
aspects  of  neutron  transport  primarily  involve  inter¬ 
mediate  neutrons,  and  their  behavior  strongly  affects 
the  design  of  high-performance  shields. 

Fast  Neutrons.  -  Neutrons  having  energies  between 
0.5  MeV  and  about  15  MeV  are  regarded  as  fast 
neutrons.  Most  neutrons  produced  by  fission  and  fusion 
reactions  are  fast  neutrons. 

Very  High  Energy  and  Ultra  High  Energy  Neutrons.  - 
Neutrons  referred  to  as  having  very  high  energy  usually 
fall  in  the  15-  to  50-MeV  range  and  tliose  referred  to  as 


having  ultra  high  energy  fall  in  the  range  >  50  MeV. 
These  neutrons  are  either  cosmic  radiations  or  by¬ 
product  particles  produced  during  accelerator  opera¬ 
tions.  They  are  not  of  interest  in  weapons  radiation 
shielding,  but  are  included  here  to  complete  the 
neutron  classifications. 


2.4.2  GAMMA  RAYS 

During  early  investigations  of  radioactive  materials 
three  types  of  radiation  were  observed  that  were 
arbitrarily  identified  as  alpha,  beta,  and  gamma  rays. 
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Gamma  rays,  unlike  the  alpha  and  beta  rays,  which 
were  found  to  carry  positive  and  negative  charges 
respectively  and  to  be  easily  stopped,  were  found  to  be 
electrically  neutral  and  highly  penetrating.  It  was 
postulated  that  gamma  rays  were  actually  electro¬ 
magnetic  radiations  similar  to  X  rays,  and  this  was 
confirmed  experimentally  in  1914  by  E.  Rutherford. 

According  to  classical  theory,  electromagnetic  radia¬ 
tion  consists  of  transverse  electric  and  magnetic  fields 
that  propagate  as  waves  at  the  speed  of  light.  They  may 
exist  in  regions  that  are  void  of  matter  and  contain  no 
electric  charges  or  currents.  Within  the  framework  of 
classical  theory  a  continuous  flow  of  energy  is  associ¬ 
ated  with  the  propagation  of  the  electromagnetic  wave, 
and  there  is  no  upper  or  lower  limit  to  the  amount  of 
energy  that  can  be  transported  by  electromagnetic 
radiations. 

While  wave  theory  gives  a  good  description  of  the 
diffraction,  reflection,  and  interference  of  electro¬ 
magnetic  radiations,  it  is  quite  inadequate  for  describing 
phenomena  such  as  nuclear  emission,  the  photoelectric 
effect,  Compton  scattering,  and  bremsstrahlung.*  How¬ 
ever,  the  quantum  theory,  which  was  formulated  by 
Max  Planck  in  1900,  does  provide  a  satisfactory 
explanation  for  such  phenomena.  This  theory  suggests 
that  electromagnetic  radiations  are  emitted  as  discrete 
bundles  of  energy  of  size  hv,  where  Planck’s  constant  h 
is  equal  to  6.624  X  10'^’  ergs/sec  and  v  is  the 
frequency.  These  pseudo-particles  are  commonly  re¬ 
ferred  to  as  photons.  They  have  energies  of  E  =  hv  and 
possess  the  properties  of  momentum  (p  =  hu/c)  and 
relativistic  mass  (M  =  hv/c^ ). 

Wave  theory  and  quantum  theory  seemingly  con¬ 
tradict  each  other  in  many  ways,  but  in  1928  Max  Born 
argued  that  electromagnetic  radiation  was  indeed 
quantum  in  nature  and  that  the  classical  wave  theory 
could  be  interpreted  as  the  probability  of  finding  a 
photon  at  a  particular  point  in  space  and  time.  It  is 


*These  phenomena  are  discussed  in  Section  2.7. 


within  the  context  of  quantum  theory  that  electro¬ 
magnetic  radiations  are  treated  in  this  chapter. 

In  the  years  since  1914  it  has  been  established  that 
electromagnetic  radiations  originate  whenever  transi¬ 
tions  take  place  between  two  energy  levels  of  nuclei  or 
between  two  electronic  levels  of  an  atom.  The  fre¬ 
quencies  and  wavelengths  of  the  radiations  are  de¬ 
termined  by 

r  =  hv  =  hc/\,  (2.82) 

where  the  level  width  1'  is  the  difference  in  the 
excitation  energy  of  the  two  levels.  Since  the  level 
widths  are  greater  for  nuclear  transitions  than  for 
electronic  transitions,  the  radiations  accompanying  nu¬ 
clear  transitions  are  characteristically  of  higher  energy 
than  those  accompanying  electronic  transitions.  For 
historical  reasons  the  radiations  produced  by  nuclear 
transitions  are  usually  referred  to  as  gamma  rays  and 
those  produced  by  atomic  transitions  are  called  charac¬ 
teristic  X  rays.  Typical  of  the  latter  are  the  soft  X  rays 
(fluorescence  radiation)  produced  as  a  result  of  the 
photoelectric  effect  discussed  in  Section  2.7.  Two  other 
types  of  electromagnetic  radiations  are  bremsstrahlung 
and  annihlation  radiation.  Bremsstrahlung  radiation 
arises  as  a  result  of  the  acceleration  or  deceleration  of 
free  charged  particles,  and  annihilation  radiation  is 
produced  when  a  positron  and  an  electron  combine, 
both  particles  being  annihilated  in  the’  process.  (Posi¬ 
trons  are  produced  by  the  pair-production  process 
described  in  Section  2.7.2). 

From  the  foregoing  it  is  seen  that  electromagnetic 
radiations  are  classified  by  origin;  that  is,  (I)  gamma 
rays  accompany  nuclear  transitions,  (2)  X  rays  are 
emitted  in  atomic  transitions,  (3)  bremsstrahlung  radia¬ 
tion  results  from  the  acceleration  or  deceleration  of  free 
charged  particles,  and  (4)  annihilation  radiation  is 
emitted  when  a  positron  combines  with  an  electron.  As 
it  turns  out,  and  as  is  discussed  in  Section  2.7, 
bremsstrahlung  and  annihilation  radiations,  as  well  as 
soft  X  rays,  are  not  usually  highly  penetrating  and  the 
transport  of  these  radiations  is  usually  ignored  in 
shielding  design . 


2.5.  Interactions  of  Neutrons 


Neutrons  are  attenuated  by  interactions  with  the 
nuclei  of  the  elements  comprising  the  medium  that  they 
are  traversing.  The  neutron  interactions  most  important 
to  shielding  are  elastic  scattering,  inelastic  scattering, 
radiative  capture,  and  charged-particle  reactions.  The 
relative  importance  of  these  interactions  is,  at  least  in 
part,  determined  by  their  respective  cross  sections,  the 
magnitude  of  which  varies  with  the  neutron  energy  and 
the  element  (see  Section  2.3). 

The  relation  of  the  various  cross  sections  is  shown  in 
Fig.  2.10.  From  this  figure  it  is  apparent  that  the  total 
cross  section  for  a  neutron  interaction  can  be  thought 
of  as  being  the  sum  of  the  elastic-scattering  cross 
section  and  the  nonelastic  cross  section.  The  two 
classifications  distinguish  between  those  interactions  in 
which  the  neutron  and  nucleus  collide  as  two  perfectly 


elastic  spheres  (sometimes  referred  to  as  billiard-ball 
collisions)  and  all  other  interactions.  The  total  cross 
section  can  also  be  thought  of  as  the  sum  of  the 
total  scattering  cross  section  and  the  total  absorption 
cross  section.  In  the  latter  interpretation  the  scattering 
cross  section  includes  both  elastic  and  nonelastic 
scatterings,  and  the  absorption  cross  section  includes 
all  processes  that  cause  the  product  nucleus  to 
differ  either  in  atomic  number  or  in  atomic  weight  from 
the  target  nucleus.  Because  of  the  change  in  the  target 
nucleus,  absorption  processes  are  usually  referred  to  as 
nuclear  reactions  rather  than  mreractions,  although  the 
two  terms  are  often  used  interchangeably. 

The  accuracy  (or  reliability)  of  a  given  calculation  can 
be  strongly  influenced  by  the  quality  of  the  cross- 
section  data  used.  Despite  their  obvious  importance, 
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Fig.  2.10.  Relation  of  Neutron  Cross  Sections.  (Note;  In  some  neutron  transport  calculations  all  neutron-nucleus  interactions  in 
which  neutrons  leave  the  target  nucleus  are  considered  to  be  scattering  processes,  in  which  case  the  n,2n  and  fission  processes  would 
be  scattering  processes,  as  indicated  by  the  dashed  lines.  In  most  shield  attenuation  calculations  neither  process  is  important, 
however.) 
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detailed  neutron  cross  sections  are  not  presented  in  this 
document  because  of  the  vastness  and  complexity  of 
the  information.  Much  cross-seciion  data  are  readily 
available  in  a  number  of  publications.^"^  Also  available 
is  a  comprehensive  and  periodically  updated  index'**  to 
the  latest  literature  on  microscopic  neutron  data. 

2.5.1.  ELASTIC  SCATTERING 

As  mentioned  above,  the  elastic  scattering  of  neutrons 
with  nuclei  can  be  viewed  as  the  collision  of  perfectly 
elastic  spheres.*  This  process  is  characterized  by  the 
conservation  of  kinetic  energy  and  momentum  and  is 
exactly  described  by  the  methods  of  classical  me¬ 
chanics.  The  various  theoretical  analyses  are  usually 
performed  in  the  center-of-mass  (c.m.)  frame  of  ref¬ 
erence,  which  assumes  that  the  center  of  mass  of  the 
neutron-nucleus  system  remains  at  rest.  An  external 
observer  (e.g.,  one  performing  a  physical  measurement) 
would  view  the  collision  process  in  the  laboratory  frame 
of  reference.  Theoretical  results  are  more  useful  when 
expressed  in  terms  of  the  laboratory  system  coordi¬ 
nates,  but  the  theoretical  relations  are  more  easily 
derived  in  the  c.m.  system. 

If  it  is  assumed  that  the  target  nucleus  in  the  labora¬ 
tory  system  is  stationary,  the  change  in  kinetic  energy 
experienced  by  the  neutron  is  uniquely  related  to  its 
change  in  direction  by  the  following  equation: 


where 


(2.84) 


a  =  |(.4  I  )KA  t  I )]  ^ 


(2.85) 


Tlie  cosine  of  the  scattering  angle  in  the  laboratory 
system,  symbolized  by  p,  is  simply  related  to  the  cosine 
of  the  scattering  angle  in  the  c.m.  system  by 


I  -I-  r}A 

^  {A^  +  2r)/l  +  1  )'/^ 


(2.86) 


If  it  is  further  assumed  that  the  scattering  in  the  c.m. 
system  is  isotropic  (an  excellent  assumption  for  hy¬ 
drogen),  then  the  following  relations  hold  for  the 
elastic-scattering  process: 

1 .  The  probability  that  a  neutron  will  undergo  ji 
change  of  direction  into  a  differential  solid  angle  t/f2 
corresponding  to  a  conical  element  defined  by  the 
scattering  angles  0  and  0  +  dO  in  the  c.m.  system  is 

dU  1 

p(0)d0= — =-mOdO.  (2.87) 

4rr  2 


2.  The  fraction  of  all  elastic-scattering  collisions  that 
involve  neutrons  of  energy  Hq  scattered  into  the  energy 
interval  dE  about  the  final  energy  E  is  given  by 


E 


(A^  +2JJ4  +  1), 


(2.83) 


p{E;Eo )  dE  = 


dE 

£■0(1  -  a)' 


(2.88) 


where 

Eq  =  initial  neutron  energy  in  the  laboratory  system, 
E  =  final  neutron  energy  in  the  laboratory  system, 

A  =  ratio  of  the  atomic  mass  of  the  nuclide  to  the 
mass  of  the  neutron, 

17  =  cos  0,  the  cosine  of  the  scattering  angle  in  the 
c.m.  system. 

The  minimum  value  of  E/Eq,  which  corresponds  to  the 
maximum  loss  in  energy  experienced  by  a  neutron  in  a 
single  collision,  occurs  for  a  head-on  collision,  0  =  rr  or 
I)  =  -  I.  The  minimum  final  neutron  energy  is  given  by 


3.  The  average  final  energy  of  all  neutrons  of  initial 
energy  £’0  is  given  by 

f  =  /  ^°  Ep{E:Eo)dE  =  (( 1  +  a)/2]  £„  -  (2.89) 

OSq 

4.  The  average  value  of  the  cosine  of  the  scattering 
angle  in  the  laboratory  system,  an  important  parameter 
in  neutron  transport  calculations,  is  given  by 


5.  The  average  logarithmic  energy  decrement  per 
collision  is  defined  by 


*Thi.s  description  applies  only  to  “potential”  elastic  scatter¬ 
ing,  which  is  the  type  discussed  in  this  section.  Another  kind  of  — —  „ 

elastic  scattering,  called  capture  scattering,  is  very  much  like  ^  =  ®  ln—p(EEo)dE=  1  |  **  *^  ** 

inelastic  scattering  in  that  it  involves  the  formation  of  a  com-  E  ’'otEo  E  '  I  -  a 

pound  nucleus  and  the  subsequent  emission  of  a  neutron  with 
nearly  its  original  kinetic  energy  (see  later  discussion). 


I 


(2.91) 
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An  interesting  property  of  ^  is  its  independence  of  the 
initial  energy  of  the  neutron.  The  quantity  ^  can  also 
be  interpreted  as  the  average  change  in  lethargy  per 
collision,  where  the  lethargy  u  is  defined  as 

i/  =  ln§,  (2.92) 

E 

and  the  average  change  in  lethargy  per  collision  is  given 
by 

M  -  «o  =  In  ^  -  In  ^  =  In  ^  =  I .  (2.93) 

t  t.Q  L 

The  lethargy  is  a  useful  alternate  energy  variable  in  the 
description  of  the  slowing-down  process  for  neutrons. 

For  most  nuclides  the  isotropic  assumption  generally 
gives  a  good  description  of  the  elastic-scattering  process 
for  neutrons  in  the  thermal-to-keV  energy  range; 
however,  it  is  invalid  for  high-energy  (MeV)  neutrons  or 
for  interactions  in  which  chemical  binding  energies  are 
important.  Thus  in  shielding  calculations,  which  usually 
involve  high-energy  neutrons,  the  isotropic  assumption 
cannot  be  made  safely  except  for  hydrogen,  for  which 
the  scattering  is  virtually  isotropic  in  the  c.m.  system 
for  neutron  energies  up  to  14  MeV. 

In  a  more  general  treatment  of  the  scattering  process, 
which  can  apply  to  inelastic  as  well  as  to  elastic 
scattering,  the  cross  sections  for  scattering  are  repre¬ 
sented  by  an  expansion  in  a  series  of  Legendre 
polynomials: 

aJE)  „  2m  +  I 

a,(£.t7)=-,-  L— ;-A,„(£l/',„(r,).  (2.94) 
where 

a^iE.Tj)  =  energy-dependent  differential  cross  section 
for  the  elastic  scattering  of  neutrons  of 
energy  E  through  an  angle  in  the  c.m. 
system  whose  cosine  is  tj  per  unit  solid 
angle, 

4, „(£)  =  energy -dependent  Legendre  coefficients 

(except  that  4o  =  1  always), 

=  Legendre  polynomials. 

The  energy-dependent  elastic-scattering  cross  section 
a^(E)  is  obtained  by  integrating  the  differential  elastic¬ 
scattering  cross  section  a^(E,Ti)  over  all  possible  values 
of  rj: 

ag(E)  =  2n  f  ‘  a^iE.ri)  Jrj .  (2.95) 

Representation  of  the  cross  section  as  an  expanded 
series  of  orthogonal  functions  is  necessary  in  the 
spherical  harmonics  apd  *9,,  solutions  of  the  Boltzmann 


equation  (see  Chapter  3).  Other  orthogonal  functions 
can  be  used  -  and  have  been  -  but  the  Legendre  poly¬ 
nomials  generally  give  best  fits  and  are  universally  used. 

Note  that  the  first  term  (zeroth  order)  of  Eq.  2.94 
gives  the  angular  cross  section  for  isotropic  scattering 
OglAu.  The  coefficients  4^(£’)  obtained  by  fitting 
angular  cross-section  data  are  reported  in  both  labora¬ 
tory  and  c.m.  systems.  Usually  the  c.m.  coefficients 
must  be  transformed  to  the  laboratory  system  for 
transport  calculations.  As  the  atomic  number  increases, 
however,  the  difference  becomes  small. 

The  elastic<ollision  process  is  an  important  energy 
degradation  mechanism  and  is  most  important  for  the 
light  elements,  particularly  hydrogen.  Efficient  neutron 
shields  always  contain  hydrogenous  materials. 


2.5.2.  NONELASTIC  REACTIONS 

All  reactions  that  are  not  elastic  scatterings  are  called 
nonelastic  reactions,  and  the  sum  of  their  cross  sections 
is  referred  to  as  the  nonelastic  cross  section.  In 
experimental  determinations  of  neutron  cross  sections 
in  the  energy  range  3  to  1 5  MeV,  the  nonelastic  cross 
section  is  the  quantity  that  is  actually  measured.  Or,  in 
other  words,  the  quantity  measured  is  the  total  cross 
section  minus  the  elastic-scattering  cross  section. 

Many  models  are  available  for  describing  nonelastic 
neutron  reactions  with  nuclei,  but  the  one  most 
frequently  used  is  the  Bohr  model.  In  this  model  the 
interaction  is  envisioned  as  beginning  when  the  neutron 
comes  within  the  range  of  the  nuclear  forces  of  the 
nucleus  and  as  terminating  when  the  products  of  the 
interaction  depart  from  this  region.  In  the  interim  the 
combination  of  the  neutron  and  target  nucleus  is  called 
simply  the  compound  nucleus.  It  is  convenient  to 
consider  the  interaction  as  occurring  in  two  distinct 
stages: 

1.  A  compound  nucleus  is  formed  when  the  incident 
neutron  loses  its  identity  and  becomes  assimilated  into 
the  unstable  new  system.  This  event  brings  to  the  newly 
formed  compound  nucleus  an  excitation  energy  that 
consists  of  the  available  portion*  of  the  neutron’s 
kinetic  energy  of  motion  and  its  binding  energyt  with 


*Duc  to  conservation  of  momentum,  some  of  the  kinetic 
energy  of  the  incident  neutron  must  appear  as  translational 
energy  at  the  compound  nucleus  and  is  therefore  unavailable  as 
excitation  energy. 

fThc  binding  energy  per  nucleon  is  the  average  energy 
required  per  nucleon  to  separate  the  nucleus  into  its  constituent 
particles.  Refer  to  subsequent  discussions  of  binding  energy 
under  "Nuclear  Reactions.” 
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respect  to  the  compound  nucleus.  The  excitation 
energy  is  statistically  shared  by  all  tlie  nucleons. 

2.  The  unstable  compound  nucleus  separates  into  its 
products  when  one  or  more  of  the  nucleons  acquire 
sufficient  energy  to  escape,  that  is,  when  the  nucleons 
acquire  energy  that  is  equal  to  or  greater  than  their 
binding  energy.  In  some  cases  tlie  unstable  compound 
nucleus  may  give  off  excitation  energy  in  the  form  of 
gamma  radiation.  In  either  case  the  resulting  product 
nucleus  may  itself  be  unstable  and  undergo  radioactive 
transformation,  but  in  all  cases  the  end  result  is  a 
system  of  stable  products.  A  schematic  representation 
of  the  two  stages  of  neutron  interaction  with  a  nucleus 
is  shown  in  Fig.  2.1 1.  (The  example  given  in  th,*  figure 
is  specifically  for  the  inelastic-scattering  proiess  de¬ 
scribed  below.) 

The  compound  nucleus  concept  is  being  replaced  by 
more  elaborate  models;  however,  it  will  continue  to  be 
useful  in  providing  an  easily  understood  representation 
of  nonelastic  neutron  interactions  with  nuclei. 

The  nonelastic  neutron  reactions  are  divided  into 
scattering  processes  and  absorption  processes,  tlie  ab¬ 
sorption  processes  often  being  referred  to  as  nuclear 
reactions.  They  include  radiative  capture,  charged- 
particle  reactions,  (n,2n)  reactions,  and  fission,  insofar 
as  shielding  attenuation  is  concerned,  the  (n,2n)  and 
fission  reactions  are  usually  unimportant  unless  ma¬ 
terials  such  as  beryllium  and  depleted  uranium  are  used. 
The  fission  reaction  is,  of  course,  very  important  as  a 
neutron  source  in  weapons  detonations  and  as  the 
energy  source  in  the  production  of  electric  power  (see 
Section  2.6). 
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Fig.  2.11.  Schematic  Repieaentation  of  a  Nonelastic  Inter¬ 
action  of  a  Neutron  with  a  Nucleus  (Inelastic  Scattering). 


Inelastic  Scattering 

In  inelastic  scattering  tiic  incident  neutron  and  target 
nucleus  fonii  a  compound  nucleus  from  which  a 
neutron  is  subsequently  emitted.  If  the  kinetic  energy 
of  the  incident  neutron  is  greater  than  the  lowest 
energy  level  of  the  target  nucleus,*  the  nucleus  may  be 
left  in  an  excited  state  following  the  neutron  emission. 
This  results  in  the  emitted  neutron  having  substantially 
less  energy  than  the  incident  neutron,  thereby  con¬ 
tributing  strongly  to  the  energy  degradation  process  and 
therefore  to  tlie  attenuation  of  neutrons  in  shields.  The 
excited  nucleus  then  returns  to  its  ground  state  by  the 
emission  of  gamma  rays.  An  energy  level  diagram  of  this 
process  is  shown  in  Fig.  2.1 2. 

Determining  the  detailed  cross  sections  for  inelastic 
scattering  is  a  complicate  I  and  difficult  task,  whether 
accomplished  by  experiment  or  by  theory.  Conse¬ 
quently,  the  inelastic-scattering  cross  sections  for  many 
elements  are  unknown  and  a  large  fraction  of  those 
used  in  shielding  calculations  are  little  more  than  a  best 
guess.  However,  there  arc  some  generalizations  that  can 
be  made; 

1.  The  cross  section  usually  increases  with  increasing 
neutron  energy. 

2.  The  cross  section  usually  increases  wiui  increasing 
atomic  weight. 

3.  The  threshold  for  the  interaction  in  iiglit  elements 
is  generally  several  MeV,  and  the  gamma  rays  emitted 
are  correspondingly  harder  than  those  produced  in 
interactions  with  lower  energy  thresholds. 

4.  The  tlireshold  for  the  interaction  in  heavy  elements 
is  lower,  sometimes  much  lower,  than  for  that  in  light 
elements. 

Nuclear  Reactions 

The  temi  nuclear  reaction  is  often  reserved  for  those 
interactions  in  which  the  product  nucleus  differs  either 
in  atomic  niimber  or  atomic  weight  from  tlie  target 
nucleus.  The  type  of  reaction  that  occurs  depends  on 
several  factors,  the  most  important  being  the  nuclear 
structure  of  the  target  nucleus  and  the  energy  of  the 
incident  neutrons.  Some  reactions  are  possible  for 
neutrons  of  any  energy,  while  others  have  threshold 
energies.  In  some  cases  more  than  one  de-excitation 
process  is  possible,  with  the  relative  probabilities  of 
each  being  given  by  their  respective  reaction  cross 
sections. 


*lf  this  condition  is  not  satisfied,  the  capture  scattering  form 
of  ela.stic  scat;-  ing  may  occur. 
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Fig.  2.1 2.  Fneigy  Diagrim  for  the  Inelastic  Scattering  of  a  Neutron  with  a  Target  Nucleus. 

•'Given  by  931  MeV  |(A/ y  +  3/„)  .Wy|  ^  7  MeV,  where  and  My  we  the  masses,  in  atomic  mass  units,  of  the  target 

nucleus,  the  neutron,  and  the  compound  nucleus,  respectively, 

^Given  by  KEf  |3f;y/(3/,v  Aff|)|.  Due  to  the  conservation  of  momentum,  some  of  the  kinetic  energy  of  the  incident  neutron 

must  appear  as  translational  energy  of  the  compound  nucleus  and  is  therefore  unavailable  as  excitation  energy. 

•■Given  by  931  (My  A/ a')  '^  ‘>31  MeV. 

•'Given  by  Ai'j  |(A/,v  +  A/„)/A/,vl . 

••Given  by  Ai’,  |Af,v/(Af.v  +  A/„)  J  -  Af  j  |(Af y  +  • 
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The  equation  for  a  reaction  between  a  target  nucleus 
and  a  neutron  with  essentially  zero  kinetic  energy  is 
given  by 

+  >'/>  +  (;.  (2.%) 

where 

=  target  nucleus  with  atomic  num¬ 
ber  Z|  and  mass  number^, 

oH*  =  incident  neutron  with  atomic 

number  0,  mass  number  I,  and 
essentially  zero  kinetic  energy, 

\  =  compound  nucleus  in  its  excited 
state, 

(Zi  Zj)  ^  “  product  nucleus, 

z^y^  ~  eniitted  particle, 

Q  -  nuclear  reaction  energy. 

Equation  2.96  requires  that  a  balance  exist  with  respect 
to  both  the  atomic  number  and  the  mass  number  and  is 
useful  in  analyzing  nuclear  reactions. 

The  nuclear  reaction  energy  Q  is  determined  by  the 
change  in  mass  experienced  by  the  individual  com¬ 
ponents  of  the  reaction.  The  original  equation  can  be 
written  as  a  total  energy  equation.  (All  masses  are 
expressed  in  terms  of  their  energy  equivalent,  with  I 
amu  being  equal  to  931  MeV.)  The  explicit  expression 
for  Q  is  given  by 

e  =  931  (%  +  A/„-(A/y+A/v)|  .  (2.97) 

where 

Q  =  nuclear  reaction  energy  (MeV), 

=  mass  of  the  target  nucleus  (amu), 

M„  =  mass  of  the  incident  neutron  (amu). 

My-  mass  of  the  product  nucleus  (amu). 

My  =  mass  of  the  emitted  particle  (amu). 

It  is  noted  that  if  Q  is  positive  the  reaction  is  exoergic 
and  neutrons  with  essentially  zero  kinetic  energy  can 
produce  the  reaction.  If,  on  tlie  other  hand,  Q  is 
negative,  the  reaction  is  endoergic  and  the  neutron  must 
possess  suffleient  kinetic  energy  to  make  up  the  deficit. 
It  can  be  shown  that  the  threshold  energy  required  of 
the  neutron  is  given  by 


£7.= 


Mv 


Q. 


(2.98) 


This  follows  from  tlie  law  of  the  conservation  of 
momentum  when  applied  to  the  formation  of  the 
compound  nucleus,  since  it  requires  that  some  of  the 
neutron's  kinetic  energy  must  appear  as  translational 
energy  of  the  compound  nucleus  and  is  therefore 
unavailable  as  excitation  energy. 

As  pointed  out  above,  those  processes  that  fall  in  the 
category  of  nuclear  reactions  all  contribute  to  the 
absorption  cross  section  (see  Fig.  2.10)  and  include 
radiative  capture,  (n,2/i)  reactions,  charged-particle  re¬ 
actions,  and  fission.  Of  these  processes,  radiative  cap¬ 
ture  is  the  most  important  experienced  by  low-energy 
neutrons  in  sliields.  Following  radiative  capture,  the 
complete  de-excitation  of  the  compound  nucleus  is 
accomplished  by  the  emission  of  gamma  radiation.  The 
binding  energy  of  the  neutron  in  the  compound  nucleus 
resulting  from  slow-neutron  capture  is  given  by 

F.),=93l  |(3/v+3/„)  My].  (2.99) 

in  this  case  the  total  excitation  energy  is  equivalent  to 
the  binding  energy  and  is  about  8  MeV. 

The  other  nuclear  reactions  become  increasingly  more 
important  as  the  energy  of  the  neutron  increases, 
because  other  forms  of  de-excitation  become  more 
probable.  However,  epithennal  radiative  captures  do 
occur  and  the  gamma  rays  thus  produced  may  con¬ 
tribute  significantly  to  or  even  dominate  the  overall 
shielding  problem. 

Because  of  limitations  in  basic  nuclear  tlieories  and 
the  scarcity  of  experimental  data,  the  standard  practice 
has  been  to  assume  that  the  energy  spectra  of  the 
epithermal  radiative  capture  gamma  rays  are  the  same  as 
those  due  to  thennal  capture.  Yost  and  Solomito' ' 
have  shown  that  this  assumption  is  inappropriate  for 
many  materials  and  tliat  the  exclusive  use  of  thermal- 
capture  spectra  can  introduce  considerable  error  in  the 
subsequent  ganuna-ray  transport  problem  and  may,  in 
fact,  be  the  most  questionable  of  the  nuclear  data  used 
in  the  calculation.  Yost'^  developed  a  method  for 
determining  capture  gamma-ray  spectra  as  a  function  of 
neutron  energy  (capture  states)  and  wrote  a  machine 
code  to  generate  capture  gamma-ray  spectra  with  a 
desired  group  structure  resulting  from  capture  of 
neutrons  of  each  energy  group.  The  code  is  based  upon 
an  appropriate  nuclear  model  involving  the  spin  and 
parity  of  the  nucleus. 

The  (n,2n)  reaction  involves  the  emission  of  two 
neutrons  and  its  occurrence  is  usually  limited  to 
neutrons  with  high  incident  kinetic  energies  -  greater 
than  10  MeV.  The  product  nucleus  is  an  isotope  of  the 
target  nucleus  and  is  frequently  radioactive.  The  thresh- 
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old  energy  for  the  in,2/i)  reaction  is  approximately 
equal  to  the  binding  energy  per  nucleon  of  the 
compound  nucleus  and  is  related  to  the  threshold 
energy  for  the  photoneutron  (y.n)  reaction  by  the 
relation 


M  y  +  3/,, 

''r(n.2n)~ 


(2.100) 


For  certain  very  light  or  very  heavy  nuclides  the 
thresholds  for  the  (n,2n)  reaction  are  sometimes  much 
less  than  10  MeV.  For  example,  JH,  3  Li,  and  ’B  have 
thresholds  of  3.34  MeV,  6.2  MeV,  and  1.85  MeV 
respectively,  and  ^S’hi,  ^2oTh,  and  ^qjUhave  thresh 
olds  of  7.4  MeV,  6.44  MeV,  and  6.0  MeV  respectively. 
The  (n,2n)  production  of  neutrons  is  usually  ignored  in 
shield  design;  however,  it  is  possible  that  this  effect 
could  be  important  for  shield  systems  containing 
significant  amounts  of  certain  materials  (beryllium,  for 
example). 

In  charged-particle  reactions  the  de-excitation  of  a 
compound  nucleus  formed  by  neutron  capture  is 
accomplished  by  the  emission  of  a  charged  particle, 
either  a  proton  or  an  alpha  particle.*  This  type  of 
reaction  is  most  probable  in  tlie  MeV  range;  however, 
there  are  several  important  charged-particle  reactions 
for  slow  neutrons.  Consider  the  following  two  methods 
for  producing  tritium: 

^  He  +  in  -  I?  He)  *  -*  f  H  +  1 H  +  0.764  MeV  K.  E. 


‘  ?B  +  in  -  I*  i  B]  •  ^  I U*  +  ; He  +  2.30  MeV  K.E. 
iLi*- J  Li +  0.478  MeV  7 

The  *°B(n,a)  reaction  is  useful  for  the  detection  of 
slow  neutrons  because  of  its  large  thermal  cross  section. 
Also,  '“B  serves  as  an  important  constituent  for  shields 
designed  to  suppress  high-energy  capture  gamma-ray 
production. 

In  fission  reactions  the  compound  nucleus  separates 
into  fission  fragments,  concomitantly  releasing  neutrons 
and  gamma  rays.  This  process  rarely  occurs  except  in  a 
few  of  the  very  heavy  elements,  which,  with  the 
exception  of  depleted  uranium,  are  not  used  in  shields. 
(See  further  discussion  of  fission  reaction  in  Section 
2.6.) 

2.5.3,  NEUTRON  TRANSPORT  CROSS  SECTION 
The  neutron  transport  cross  section  is  defined  as 

=  a^{d)co&0dn,  (2.101) 

where 


a,  =  total  microscopic  cross  section, 

0^(0)  *  differential  microscopic  cross  section  for  elas¬ 
tic  scattering  through  a  given  angle  6  per  unit 
solid  angle, 

6  =  scattering  angle  in  the  laboratory  frame  of 
reference. 


Li  +  -*  (I  Li]  *  -  ?  H  +  t  He  +  4.785  MeV  K.E. 

Another  important  charged-particle  reaction  involves 
the  '°B  isotope  of  boron: 

*An  alpha  particle  is  Hc^. 


The  reaction  rate  <t>  may  be  regarded  as  an  estimate 
of  the  rate  at  which  neutrons  lose  their  forward 
momentum.  In  this  context  the  transport  cross  section 
has  been  used  as  an  approximation  to  the  energy- 
dependent  removal  cross  section  (see  Chapter  3,  Section 
3.8.2).  Quite  often,  o,,  appears  naturally  in  the  results 
of  a  derivation'  ^  which  has  other  formal  bases. 


2.6.  Neutron  Production  Processes 


Except  for  the  small  contributions  made  by  such 
reactions  as  the  (n,2n),  («,3w),  and  (7/1)  reactions,  the 
neutrons  emitted  in  the  detonation  of  a  nuclear  weapon 
are  produced  directly  by  the  fission  process  or  the 
fusion  process.  In  the  fission  process  the  resulting 
neutrons  cover  a  wide  range  of  energies,  whereas  in  the 
fusion  process  the  emissions  are  essentially  mono- 
energetic,  with  the  energies  being  determined  by  the 
type  of  fusion  reaction. 

2.6.1.  FISSION  REACTION 

The  details  of  the  fission  process  are  very  complex 
and  only  a  brief  description  is  presented  here.*  As  in 
other  nonelastic  neutron  interactions  (see  Section 
2.5.2),  fission  begins  with  the  formation  of  a  compound 
nucleus  in  an  excited  state.  If  the  excitation  energy 
exceeds  a  critical  energy,  then  fission  may  occur.^ 
Calculated  values  of  the  fission  threshold  energy  are 
shown  in  Table  2.2  for  a  wide  range  of  atomic 
weights.**  It  is  apparent  that  only  the  very  heavy 
nuclides  (A  >  230)  have  reasonably  low  threshold 
energies. 

Experimental  values  of  the  fission  thresholds  for  the 
nuclear  fuel  nuclides  of  practical  interest  are  given  in 
Table  2.3  (ref.  14).  The  significance  of  a  negative 
threshold  energy  “’U,  and  ^^’Pu)  is  that 

neutrons  with  essentially  zero  kinetic  energy  can 
produce  the  fissioning  of  these  nuclides  (thermal 
fission).  All  the  other  nuclides  (”^Th, 

and  ^^’Np)  can  undergo  fast  fission  only. 

When  fission  occurs,  the  highly  unstable  compound 
nucleus  almost  always  splits  into  two  fission  fragments. 
Splitting  into  three  or  more  fragments  is  theoretically 


^Although  this  discussion  is  limited  to  the  fission  of  the  heavy 
nuclei  by  neutrons,  it  is  pointed  out  that  any  nucleus  except 
hydrogen  may  undergo  fission  if  its  level  of  excitation  is 
sufficiently  high  and  if  this  excitation  can  be  imparted  to  the 
target  nucleus  by  its  interaction  with  gamma  radiation  (photo- 
fission)  or  with  a  variety  of  high-energy  particles. 

'I'The  alternative  to  the  fission  mode  of  de-excitation  is  the 
emission  of  gamma  radiation,  which  returns  the  .;ompound 
nucleus  to  its  ground  state,  Le.,  effectively  radiative  capture. 


possible,  but  such  an  occurrence  is  rare  and  is  generally 
ignored.  The  nucleons  are  more  firmly  bound  within 
their  respective  fission  fragments  than  they  were  in  the 
original  fissile  nuclide,  thereby  liberating  the  very  large 
energy  of  fission  (about  200  MeV  per  fission  for  the 
very  heavy  nuclides).  Fissioning  usually  results  in  the 
prompt  emission  of  several  neutrons,  the  average 
number  released  per  fission  being  a  function  of  the 
neutron  energy  and  the  fissile  material.  For  example, 
fissioning  of  and  ^^’Pu  by  thermal  neutrons 

produces  about  2.5  and  2.9  neulrons/fission  respec- 


Table  2,2.  Computed  Neutron  Fission  Thresholds 
as  a  Function  of  Nuclear  Mass* 


Mass  No., 

Fission  Threshold 

A 

(MeV) 

16 

18.5 

60 

48 

100 

47 

140 

62 

200 

40 

236 

~5 

‘From  ref.  14,  p. 

.  107. 

Table  2.3.  Experimental  Neutron  Fission  Thresholds 

of  Nudear  Fuel  Nudides 

Target 

Compound 

Fission  Threshold 

Nucleus 

Nucleus 

(MeV) 

“*Th 

233Th 

1.3 

234^, 

<0 

235y 

0.4 

<0 

■>“u 

0.8 

238u 

239y 

1.2 

*”Np 

”*Np 

0.4 

»’Pu 

”°Pu 

<0 

‘From  ref.  14,  p.  108. 
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lively.  In  addition,  some  of  the  fission  fragments 
contain  nuclides  .Wth  neutrons  in  excess  of  those 
required  for  nuclear  stability.  These  unstable  nuclides 
have  half-lives  up  to  about  1  min,  and  the  neutrons 
emitted  are  called  delayed  neutrons.*  (It  is  these 
neutrons  that  make  control  of  a  nuclear  reactor 
relatively  easy.  If  all  neutrons  were  prompt,  every 
reactivity  excursion  would  be  prompt  critical,  and  a 
reactor  would  be  difficult  to  operate  safely.)  Also, 
because  the  fission  fragments  are  neutron-rich,  they 
undergo  radioactive  decay,  emitting  primarily  beta  and 
gamma  radiations. 

The  mass  numbers  of  the  fission  fragments  versus  the 
corresponding  fission  yields  are  shown  in  Fig.  2.13  for 
the  thermal  fission  of  ’^*U.  and  ^^®Pu  (ref.  15). 
The  fission  yield  is  defined  as  the  percentage  of  the 
total  fissions  yielding  products  of  a  given  mass  number. 
As  the  energy  of  the  incident  neutron  increases,  the 
probability  of  producing  two  fission  fragments  of 
nearly  the  same  size  increases  and  the  characteristic  dip 
of  the  thcrmal-fission-yield  curves  at  a  mass  number  of 
about  1 20  disappears. 

The  energy  spectrum  of  fission  neutrons  has  been 
experimentally  detennined  for  ”*U,  tiie  nuclide  used 
most  frequently  as  a  nuclear  fuel,  and  is  most  con¬ 
veniently  represented  by  formulas  known  as  the  fission 
spectra  of  Watt,'*  of  Cranberg  et  al..'’’  and  of  Los 
Alamos.'  * 

The  formula  for  the  fission  spectruir  as  determined 
by  Watt  is 

mi  =  0.484  sinh  y/^:  ,  (2.102) 

where  N(hl  is  tlie  fraction  of  neutrons  emitted  per 
fission  with  an  energy  in  MeV,  per  unit  energy  range. 
This  fomiula  is  regarded  as  adequate  for  the  energy 
range  0.075  to  17  MeV,  which  includes  the  energy 
range  of  importance  to  shielding  (3  to  17  MeV).  Table 
2.4  presents  values  of  mi<  together  with  two  other 
useful  quantities,  F{E)  and  1  -  F{hl''  is  the 
fraction  of  neutrons  with  energies  greater  than  the 
energy  h  and  is  given  by 

/•TA’)  =  f  °°  NiE")  JH' .  (2.1 03) 

The  quantity  1  -  F{E)  is  the  fraction  of  neutrons  with 
energies  less  than  the  energy  E  and  is  given  by 


*l‘or  thermal  fissioning  of  and  the  delayed 

neutron  fraction  is  about  0.77n  and  0.27.  respectively. 


\-F{E)  =  Ni^)  dE’ .  (2.104) 

The  formula  for  the  fission  spectrum  obtained  by 
Cranberg  et  al. '  ’  is 


A^(£)  =  0.453  sinh  (2.105) 

The  measurements  on  which  this  formula  is  based 
involved  an  energy  range  of  0.18  to  12.0  MeV  and 
represent  no  substantial  deviation  from  Watt’s  earlier 


Table  2.4.  Watt  Prompt-Neutron  Spectrum  Due  to 
Fiaaon  of  ^  ^  *  U  by  Thetmal  Neutrons'' 


i 

(MeV) 

NiE) 

(MeV-') 

\-F(E) 

F{E) 

0.1 

0.2023 

0.014 

0.9860 

0.2 

0.2676 

0.377 

0.9623 

0.3 

0.3068 

0.663 

0.9337 

0.4 

0.3300 

0.987 

0.9013 

0.5 

0.3450 

0.1321 

0.8679 

0.6 

0.3528 

0.1669 

0.8331 

0.7 

0.3557 

0.2027 

0.7973 

0.8 

0.3542 

0.2385 

0.7615 

0.9 

0.3513 

0.2729 

C.7271 

1.0 

0.3446 

0.3082 

0.6918 

1.1 

0.3368 

0.3425 

0.6575 

1.2 

0.3271 

0.3759 

0.6241 

1.3 

0.3179 

0.40/4 

0.5926 

1.4 

0.3073 

0.4383 

0.5617 

1.6 

0.2841 

0.4983 

0.5017 

1.8 

0.2604 

0.5631 

0.4369 

2.0 

0.23/1 

0.0028 

0.3972 

2.5 

0.1839 

0.7073 

0.2927 

3.0 

0.1384 

0.7876 

0.2124 

3.5 

0.1026 

0.8476 

0.1524 

4.0 

0.07472 

0.8914 

0.1086 

4.5 

0.05381 

0.9192 

0.0808 

5.0 

0.03847 

0.9463 

0.0537 

5.5 

0.02729 

0.9624 

0.0376 

6.0 

0.01916 

0.9739 

0.0261 

6.5 

0.01336 

0.9821 

0.0179 

7.0 

9.316(-3)*’ 

0.9876 

0.0124 

7.5 

6.436(-3) 

0.9916 

0.0084 

8.0 

4.423(-3) 

0.9943 

0.0057 

8.5 

3.034(-3) 

0.9960 

0.0040 

9.0 

2.076(-3) 

0.9973 

0.0027 

9.5 

1.418(-3) 

0.9982 

0.0018 

10.0 

9.53(-4) 

0.9988 

1.21 1(-3) 

11.0 

4.384(-4) 

0.9994 

5.501(-4) 

12.0 

1.994(-4) 

0.9998 

2.473(-4) 

13.0 

8.953(-5) 

0.9999 

1.103(-4) 

14.0 

3.997(-5) 

1.0 

4.875(-5) 

15.0 

1.771(-5) 

1.0 

2.148(-5) 

“From  ref.  16. 

'’Read:  9.316  X  10■^ctc. 
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Table  2J,  Cranbeig  e/  aL  Prompt-Neutron 
Spectrum  Due  to  Fission  of  ^  U  by 
Thermal  Neutrons'^ 


E 

(MeV) 

V(£-) 

(MeV**) 

F{E) 

G(E) 

(MeV) 

0.00 

0.000 

1.000 

1.981 

0.25 

0.290 

0.948 

1.973 

0.50 

0.347 

0.867 

1.942 

0.75 

0.358 

0.779 

1.887 

1.00 

0.347 

0.690 

1.810 

1.25 

0.325 

0.606 

1.715 

1.50 

0.298 

0.528 

1.608 

1.75 

0.268 

0.457 

1.493 

2.00 

0.239 

0.394 

1.374 

2.50 

0.184 

0.2884 

1.138 

3.00 

0.138 

0.2082 

0.919 

3.50 

0.102 

0.1486 

0.726 

4.00 

0.0738 

0.1050 

0.563 

4.50 

0.0528 

0.0736 

0.430 

5.00 

0.0375 

0.0512 

0.324 

5.50 

0.0263 

0.03544 

0.241 

6.00 

0.0184 

0.02439 

0.178 

6.50 

0.0127 

0.01669 

0.1 3C 

7.00 

8.77(-3)'’ 

0.01138 

9.41(-2) 

7.50 

6.0 1(  3) 

7.72(-3) 

6.76(-2) 

8.00 

4.10(-3) 

5.22(-3) 

4.83(-2) 

8.50 

2.79(-3) 

3.52(-3) 

3.43(-2) 

9.00 

1.88(-3) 

2.364(-3) 

2.42(-2) 

9.50 

1.27(-3) 

1.583(-3) 

1.70(-2) 

10,00 

8.56(-A) 

1.058(-3) 

l.l87(-2) 

10.50 

5.74(-4) 

7.05(-4) 

8.25V-3) 

11.00 

3.84(-4) 

4.686(-4) 

5.72(-3) 

!  1..50 

2.56(-4) 

3.108{-4) 

3.947(-3) 

12.00 

1.70(-4) 

2.058(-4) 

2.714(-3) 

13.00 

7.48(-5) 

8.97(-5) 

1.271(-3) 

14.00 

3,26(-5) 

3.88(-5) 

5.88(-4) 

15.00 

1.41  (-5; 

1.67(-5) 

2.69(-4) 

16.00 

6.07(-6) 

7.10(-6) 

1.22(-4) 

17.00 

2.59(-6) 

3.00(-6) 

5.48(-5) 

18.00 

1.10(-6) 

1.26(-6) 

2.44(-5) 

“From  ref.  18,  p.  48. 
'’Read:  8.77  X  10'^.  etc. 


results.  Values  of  N(E)  and  F{E)  as  calculated  from 
Cranbcrg’s  spectrum  are  presented  in  Table  2.5,  along 
with  values  of  G(E),  which  is  the  energy  per  fission 
carried  by  neutrons  having  energies  greater  than  E  and 
is  given  by 

G(E)=  f°°  E' N{E')dE' .  (2.106) 


The  formula  developed  by  Los  Alamos’®  for  the 
fission  spectrum  in  the  energy  range  below  9  MeV  is 

yV(£)  =  0.77£’‘/2  e-o-7  76^  .  (2.107) 

This  equation  is  a  simpler  approximation  to  the 
experimental  data  and  is  used  primarily  for  reactor  core 
calculations.  The  error  associated  with  the  formula  is 
less  than  12%  for  the  specified  energy  range.  For  the 
energy  range  4  to  14  MeV  the  Los  Alamos  formula  has 
been  reduced  to  the  simple  exponential 

N(E)=  \J5  .  (2.108) 

The  error  associated  with  this  formula  is  always  less 
than  1 5%  and  is  less  than  7%  for  energies  from  5  to  1 3 
MeV. 

2.6.2.  FUSION  REACTION 

The  fusion  reaction  involves  the  nuclear  combination 
of  two  light  nuclei  having  sufficient  kinetic  energy  to 
overcome  the  repulsive  force  that  exists  between  nuclei. 
Because  the  mean  binding  energy  per  nucleon  for 
certain  light  nuclei  is  less  than  that  for  nuclei  of 
intermediate  mass  (essentially  the  same  change-in-mass 
phenomenon  which  characterizes  the  fission  process), 
the  fusing  process  of  these  nuclei  is  accompanied  by 
considerable  energy  release  and  by  the  emission  of 
nucleons. 

The  fusion  processes  of  practical  interest  involve  the 
light  nuclei  deuterium,  tritium,  and  ®He.  Under  proper 
conditions  several  combinations  of  these  nuclei  can 
undergo  fusion,  the  reactions  being  represented  in 
equation  form  as  follows: 

?D  +  ?D  ->■  |He(0,82  MeV)  +  in(2.45  MeV; 

?D  +  ?D^?T(1.01  MeV)+  |H(3.02  MeV) 

?D  +  ?T->-tHe(3.5  MeV)  +  i«(l4.l  MeV) 

?  D  +  J  He  ->■  5  He(3 .6  Me V)  +  I  H(  1 4 .7  MeV) 

It  is  the  2.45-  and  14.1 -MeV  neutrons  produced  in  the 
D-D  and  D-T  reactions  utilized  in  nuclear  detonations 
that  are  primarily  of  concern  in  shield  design. 

Unlike  in  tlie  fission  process,  no  prompt  gamma  rays 
are  emitted  in  the  fusion  process  and  neither  are  the 
products  radioactive.  Bremsstrahlung  radiation  (see 
Section  2.4.2)  is  produced  in  conjunction  with  charged- 
particle  deflections  by  the  Coulomb  field  of  other 
charged  particles,  but  such  radiation  does  not  con¬ 
tribute  significantly  to  the  shielding  problem. 


2.7.  Interactions  of  Gamma  Rays 


There  are  many  possible  processes  by  which  gamma 
rays  may  interact  with  matter;  however,  the  attenuation 
of  gamma  rays  in  a  shield  is  completely  dominated  by 
three  processes:  the  photoelectric  effect,  pair  produc¬ 
tion,  and  Compton  scattering.  The  three  processes  are 
described  below,  as  are  a  number  of  processes  that  are 
usually  neglected  in  shielding  calculations. 

2.7.1.  PHOTOELECTRIC  EFFECT 

The  photoelectric  effect  is  characterized  by  the 
complete  absorption  of  the  gamma  ray.  The  gamma  ray 
interacts  directly  with  an  orbital  electron,  which  is 
initially  bound  in  an  atom.  The  momentum  is  conserved 
by  the  recoil  of  the  entire  atom,  thereby  allowing  a 
complete  loss  of  energy  by  the  gamma  ray.  Conse¬ 
quently,  the  more  tightly  bound  electrons  have  the 
greatest  probability  of  absorbing  the  gamma  ray.  It  has 
been  shown*’  that  about  80%  of  the  photoelectric 
absorptions  involve  the  A'-shell  electrons.  The  orbital 
electron  (photoelectron)  is  ejected  from  the  atom 
having  a  kinetic  energy  T  =  where  Bg  is  the 

binding  energy  of  the  ejected  electron.  The  remainder 
of  the  energy,  which  is  equal  to  Bg,  appears  as  very  soft 
characteristic  X  rays,  usually  referred  to  as  fluorescence 
radiation. 

The  photoelectric  psocess  does  not  lend  itself  to 
theoretical  analysis  and  most  cross-section  data  for  this 
interaction  are  empirical,  with  theoretical  methods 
being  employed  for  interpolation  and  extrapolation  of 
the  data.  The  qualitative  dependence  of  the  probability 
for  the  photoelectric  effect  is  given  by  the  approxi¬ 
mation 


w^iere  E^  is  the  incident  photon  energy  and  Z  is  the 
atomic  number.  Equation  2.109  shows  that  the  photo¬ 
electric  effect  is  important  for  low-energy  gamma  rays 
and  large-atomic-number  materials,  such  as  lead.  Table 
2.6  gives,  for  the  elements  of  greatest  importance  to 
shielding,  the  energy  at  which  the  photoelectric  effect 


Table  2,6.  Gamma-Ray  Energies  at  Which 
Photoelectric  Effect  Provides  One-Half 
the  Total  Absorption  Coefficient* 


Atomic  No., 

Z 

Element 

Gammu-Kay 

Energy 

(MeV) 

A-Shell 

Ionization 

Energy 

(MeV) 

1 

Hydrogen 

10"^ 

1.4  X  10‘* 

4 

Beryllium 

0.011 

2.2  X  10^ 

6 

Carbon 

0.016 

2.8  X  10'“ 

8 

Oxygen 

0.025 

5.2  X  10"^ 

13 

Aluminum 

0.046 

1.5  X  10"^ 

20 

Calcium 

0.079 

4,0  X  10'^ 

26 

Iron 

0.11 

6.9  X  10'^ 

42 

Molybdenum 

0.195 

2.0  X  10'^ 

50 

Tin 

0,25 

2.9  X  10'^ 

74 

Tungsten 

0.42 

6.06  X  10'* 

82 

Lead 

0.50 

8.8  X  10'* 

92 

Uranium 

0.62 

1.16  X  lO"' 

•From  ref,  18,  p.  144, 


provides  one-half  the  total  absorption  coefficient,  along 
with  the  ionization  energy  of  the  K  shell. 

The  fluorescence  radiation  associated  with  the  photo¬ 
electric  effect  is  significant  only  when  tlie  interaction 
occurs  in  heavy  elements.  Even  then  the  range  of  the 
emitted  X  rays  is  comparable  to  that  of  the  photo- 
electroas,  and  their  contribution  to  the  penetrating 
dose  is  so  small  that  they  can  be  assumed  to  be 
absorbed  at  their  point  of  origin. 

2.7.2.  PAIR  PRODUCTION 

The  pair-production  process  has  a  threshold  energy  of 
1.02  MeV  and  becomes  increasingly  more  probable  as 
the  gamma-ray  energy  increases  above  1.02  MeV.  In  this 
process  the  gamma  ray  is  completely  absorbed  and  an 
electron-positron  pair  appears  in  its  place.  Pair  produc¬ 
tion  must  occur  within  the  electrostatic  field  of  a 
charged  particle.  The  momentum  is  conserved  in  the 
process  by  the  recoil  of  the  charged  particle,  thereby 
making  possible  the  complete  absorption  of  the  gamma 
ray.  Either  a  nucleus  or  an  atomic  electron  can  provide 
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the  necessary  electrostatic  field;  however,  due  to  the 
greater  charge  of  a  nucleus,  the  interactions  are  pri¬ 
marily  nuclear. 

The  sum  of  the  kinetic  energies  of  the  electron  pair  is 
equal  to  the  initial  gamma-ray  energy  less  the  energy 
equivalent  of  their  rest  masses  and  the  kinetic  energy  of 
the  recoil  nucleus: 

7’+  +  7’_=£'^-2moc' -T^ecoii-  (2.110) 


where 

=  5.80  X  10'^®  cm*/nucleus  ,  (2.113) 

/HoC*  =  energy  equivalent  of  the  rest  mass  of  an 
electron 
=  0.5110MeV, 


The  kinetic  energy  acquired  by  the  recoil  nucleus  can 
be  neglected  because  of  the  relativ  '  mass  of  the 

nucleus.  If,  however,  the  process  >  ‘he  field  of 

an  atomic  electron,  considerable  er  -insferred  to 

the  recoil  electron  and  a  triplet  is  p  The  triplet 

consists  of  the  newly  created  elect.  pair  (one 
positron  and  one  electron)  and  of  the  recoil  electron. 
Because  of  momentum  considerations,  the  threshold 
energy  for  triplet  production  is  2.04  MeV. 

Every  positron  created  in  pair  production  loses  its 
kinetic  energy  by  ionization  and  then  ultimately  inter¬ 
acts  with  an  ordinary  electron,  converting  the  electron 
and  itself  into  electromagnetic  radiation  called 
annihilation  radiation.  Thus  for  each  positron,  two 
0.51 -MeV  gamma  rays  are  emitted  in  opposite  direc¬ 
tions.  These  secondary  gamma  rays  are  most  important 
for  high-energy  gamma-ray  sources  (greater  than  6 
MeV)  and  medium-to-heavy  elements.  They  are  usually 
assumed  to  be  absorbed  at  the  point  of  formation; 
however,  in  the  more  refined  Monte  Carlo  and  S„ 
calculations,  transport  of  the  0.51 -MeV  annihilation 
gamma  rays  can  be  included. 

Aside  from  the  effects  of  screening,*  the  nuclear 
pair-production  cross  sections  are  theoretically 

and  experimentally  proportional  to  Z^.  Analytical 
expressions  are  available  only  for  two  special  cases: 


nto  =  electron  rest  mass 

=  0.5488  X  10‘®  amu  =9.108  X  10"^®  g, 
c  =  speed  of  light  in  a  vacuum 
=  2,99793  X  in“>  cm/sec. 

Table  2.7  gives  the  gamma-ray  energies  at  which  pair 
production  provides  one-half  the  total  absorption  coef¬ 
ficient  for  elements  of  greatest  interest  in  shielding 
calculations. 

Electronic  pair  production  can  occur  when  the 
gamma-ray  energy  exceeds  2.04  MeV.  The  analysis  of 
this  effect  is  difficult  because  of  screening  by  other 
electrons  and  by  the  nucleus.  However,  the  ratio  of  the 
total  cross  section  for  Z  electrons  of  a  given  atom, 
Zoppj;,  to  the  cross  section  for  nuclear  pair  production, 
Oppfj,  can  be  written  as 


OppN  CZ  ' 


(2.114) 


Table  2.7,  Gamma-Ray  Eneipea  at  Which  Pair 
Production  Provides  One-Half  the  Total 
Absorption  Coefficient* 


(1)  If  screening  is  neglected  and  nioc^  <  <  137 

ntoC^  Z“'/^,  then 


Opu  “  OqZ 


28,  /2E\  218 
9  UocV  27 


(2.111) 


(2)  If  complete  screening  is  considered  and  E  >  137 
nioc^ 


o,,p  =  OoZ^lfln(183Z-tP)-^ 


,  (2.112) 


Atomic  No., 

Z 

Element 

Gamma-Ray 

Energy 

(MeV) 

1 

Hydrogen 

78 

4 

Beryllium 

35 

6 

Carbon 

28 

8 

Oxygen 

20 

13 

Aluminum 

15 

20 

Calcium 

12 

26 

Iron 

9.5 

42 

Molybdenum 

7.5 

50 

Tin 

6.5 

74 

Tungsten 

5.2 

82 

Lead 

5.0 

92 

Uranium 

4.8 

*The  positive  charge  of  the  nucleus  can,  in  effect,  be  reduced 
by  the  cloud  of  negatively  charged  orbital  electrons. 


*From  ref.  18,  p.  145. 
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where  the  constant  C  depends  upon  and  varies  from 
about  2.6  at  =  6.5  MeV  to  about  1.2  at  E^  =  100 
MeV. 

The  effect  of  the  electronic  pair  production  can  be 
included  in  a  total  pair-production  cross  section  as 

Opp  =  OppN  +  ^Oppj.  =  OppN  ^1  '  (-••15) 

where  1/CZ  =  0  <  2.04  MeV. 

2.7.3.  COMPTON  SCATTERING 

The  Compton  scattering  process,  named  for  the 
'experimentalist  A.  H.  Compton,  can  be  described  as  an 
elastic  collision  between  the  incident  gamma  ray  and  a 
free  electron.*  Under  these  conditions  the  kinetic 
energy,  as  well  as  the  momentum,  is  conserved,  which 
results  in  a  precise  relation  between  the  loss  of 
gamma-ray  energy  and  its  change  in  direction.  The 
functional  representation  of  this  relation  is  given  by 

E - - ,  (2.116) 

1  -t  — ^  ( 1  -  cosfl) 

where 


Eo  =  initial  gamma-ray  energy, 

E  =  final  gamma-ray  energy, 

6  =  angular  change  experienced  by  the  gamma 
ray, 

WoC*  =0.5 110  MeV, 

mo  =0.5488  X  10'^  amu  =9.108  X  10'^*  g, 
c  =  2.99793  X  lO'”  cm/sec. 

A  plot  of  Eq.  2.1 16,  which  is  referred  to  as  Compton’s 
equation,  is  shown  in  Fig.  2.14. 

Another  useful  form  of  Compton’s  equation  is  ob¬ 
tained  if  the  gamma-ray  energy  E  and  wavelength  X  are 
expressed  in  units  of  the  electron  rest -mass  energy  (mc^ 
=  0.51 10  MeV)  and  the  Compton  wavelength  (h/mc  = 
0.02426  A)  respectively.  In  these  units  the  wavelength 
and  energy  variables  are  simply  reciprocals  of  each 
other  and  Compton’s  equation  assumes  the  form 

X-Xo  =  l-cosfl,  (2.117) 

where  X  is  the  final  wavelength. 


*As  the  energy  of  the  incident  gamma  ray  increases,  even  the 
tightly  bound  orbital  electron  behave  as  free  electrons. 
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Fig.  2.14.  Reduction  in  Gamma-Ray  Energy  by  Compton 
Scattering,  (fo  =  initial  energy =  reduced  energy.) 


The  differential  cross  section  o(0)  for  Compton 
scattering  is  defined  as  the  cross  section  per  electron  for 
scattering  through  a  given  angle  0  per  unit  solid  angle. 
Expressed  in  Thomson  units,  this  cross  section  is  given 
by  the  Klein-Nishina  forinula^^*** 


\f)‘nEl  \E  Eo 


(2.118) 


'i'This  form  of  the  equation  does  not  include  polarization 
effects  that  are  due  to  scattering.  However,  this  effect  is  of  no 
consequence  to  practical  problems  in  radiation  transport. 
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where  1  Thomson  unit  =  (8;r/3)  {e^jmc^Y  -  0.665 
bam.  The  dependence  of  the  differential  cross  section 
on  the  scattering  an^e  d  for  various  initial  ene^es  is 
illustrated  in  Fig.  2.15.  It  is  noted  that  a(d)  is 
symmetric  only  for  the  nonrelativistic  limit  £  *♦  0.  As 
the  initial  energy  increases,  the  cross  section  becomes 
more  and  more  sharply  peaked  in  the  forward  direction. 
Consider  an  initial  energy  E  =  2.5  MeV;  the  cross 
section  decreases  to  one-half  its  initial  value  at  a 
scattering  angle  of  about  20°  and  then  drops  to  about 
3%  of  its  initial  value  for  6  =  180°.  This  characteris¬ 
tically  forward  scattering  of  the  high-energy  gamma 
rays  strongly  influences  the  nature  of  the  deep- 
penetration  problem. 

The  total  cross  section  for  Compton  scattering  per 
electron  is  obtained  by  integrating  the  differential  cross 
section  over  47r  solid  angle: 

r  0(d)  c/n.  (2.119) 

^  4ff 


where  dO,  =  Iv  sind  dd.  Substituting  Eq.  2.1 18  into  Eq. 
2.1 19  and  integrating  provides  the  following  expression 
for  the  total  cross  section: 


2£’(1  +E) 
1  +  2E 


-In  (1  +  2£) 


4 


2E 


In(l  +  2E)- 


1+3E  ] 

(1  +2£)^J 


(2.120) 


where  the  units  of  the  gamma-ray  energy  E  are  given  in 
terms  of  the  electron  rest-mass  energy  (moC^  =  0.51 10 
MeV).  For  small  gamma-ray  energies  (E  <  1  nioc^)  the 
value  of  Oq(E)  approaches  one  Thomson  unit  and  for 
large  gamma-ray  energies  it  approaches  (%£)  (In  2E  + 
%). 

The  cross  section  related  to  the  average  fractional 
energy  carried  off  by  the  scattered  gamma  ray  is 
defined  as 
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Fig.  2.15,  Differential  Klein-Nishina  Cross  Sections. 
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which  on  integration  yields 


ln(H-2£)  2(1  +E)i2E^  -  2E  -  \) 
.  £•*  £*(l+2£)’ 


8£* 


3(1  +  2£)^J 


.  (2.122) 


The  other  cross  section  of  importance  relates  to  the 
average  fractional  energy  lost  by  a  gamma  ray  suffering 
a  Compton  collision;  the  fraction  increases  mono- 
tonically  as  the  gamma-ray  energy  increases.  This  cross 
section  represents  the  energy  absorption  or  energy 
deposition*  cross  section  and  is  defined  by 

OcP=r”  ^^(Ke)dn,  (2.123) 


*Thc  temi  energy  deposition  was  suggested  by  A.  B.  Chilton 
of.  the  University  of  Illinois  as  a  replacement  for  the  term 
energy  absorption  since  the  ICRU  has  designated  the  energy 
absorption  coefficient  to  mean  an  essentially  similar  quantity 
but  with  correction  terms  to  account  for  subsequent  transport 
of  the  energy  alter  having  made  an  absorption  collision.  A 
discussion  of  these  considerations  is  given  in  ORNL-RSIC-16 
(ref.  21).  One  must  be  careful  in  the  use  of  such  correction 
terms,  which  should  be  consistent  with  the  transport  model 
assumed  or  energy  may  not  be  conserved. 
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which  is  also  simply  expressed  as 

o^{E)  =  a,,(h:)-a^.ih-).  (2.124) 

It  is  noted  that  the  various  Compton  cross  sections 
are  usually  given  on  a  per  electron  basis  since  the 
scattering  occurs  with  free  electrons.  The  cross  sections 
are  simply  multiplied  by  the  atomic  number  Z  to  obtain 
the  atomic  cross  section.  More  detailed  descriptions  of 
the  Compton  cross  sections  are  presented  else¬ 
where.* 

2.7.4.  NEGLECTED  PROCESSES 

As  pointed  out  ;tbove,  the  tlir^e  processes  described  in 
the  preceding  paragraphs  are  usually  the  only  gamma- 
ray  interactions  considered  in  shielding  attenuation 
calculations  since  all  other  processes  make  only  negli¬ 
gible  contributions.  Some  of  the  interactions  that  are 
ignored  are  described  below. 

Rayleigh  Scattering 

Of  the  neglected  processes  Rayleigh  scattering  is  the 
most  important  and  of  some  slight  significance  at  low 
gamma-ray  energies,  but  it  is  still  less  probable  than 
Compton  scattering  by  1  or  2  orders  of  magnitude  for 
energies  greater  than  0.1  MeV.  Rayleigli  scattering  is 
elastic  coherent  scattering  by  tightly  bound  electrons, 
and  the  process  is  more  probable  with  heavy  elements. 
In  the  energy  range  of  interest  to  shielding,  Rayleigh 
scattering  results  in  small  changes  in  direction  and  in 
essentially  no  energy  loss.  The  net  result  of  coherent 
scattering  is  to  make  the  angular  distribution  of 
low-energy  scattered  gamma  rays  more  nearly  isotropic. 
Tabulations  of  total  attenuation  coefficients  for  practi¬ 
cal  design  calculations  generally  do  not  include  the 
Rayleigh  scattering  cross  section,  which  neglect  tends  to 
produce  a  slight  overestimate  of  the  gamma-ray  flux  at 
field  points  removed  from  the  source. 

Thomson  Scattering 

Thomson  scattering  is  quite  similar  to  Rayleigh 
scattering,  with  the  interaction  involving  the  nucleus 
rather  than  atomic  electrons.  The  effect  is  very  small 
and  is  truly  negligible. 

Nuclear  Resonance  Scattering 

Nuclear  resonance  scattering  involves  the  excitation 
of  a  nuclear  level  by  the  complete  absorption  of  an 
incident  gamma  ray  followed  by  the  re-emission  of  the 


excitation  energy.  The  probability  for  nuclear  reso¬ 
nance  scattering  is  extremely  small  and  its  effect  would 
be  to  increase  the  isotropy  of  tlie  gamma-ray  flux. 

Photodisintegration 

Photodisintegration  of  nuclei  is  possible  if  the 
gamma-ray  energy  exceeds  the  separation  energy  of  a 
nucleon.  Except  for  the  ’Be(7,;i)  and  ^H(7.«)  reac¬ 
tions,  this  effect  is  confined  to  gamma  rays  with 
energies  in  excess  of  8  MeV.  The  cross  sections  for 
photodisintegration  arc  usually  negligible  when  com¬ 
pared  with  those  for  Compton  scattering  and  pair 
production. 

2.7.S.  GAMMA-RAY  CROSS  SECTIONS 

The  total  gamma-ray  interaction  cross  section  for  the 
usual  shielding  calculation  is  given  by 

o,  =  0,,,.;  +  o,,,.  +  ZO(.  ,  (2.125) 

where 

Of  =  total  gamma-ray  interaction  cross  section  per 
atom, 

0|,|.;  =  photoelectric  microscopic  cross  section  per 
atom, 

Oj,,,  =  pair-production  microscopic  cross  section  per 
atom, 

a^.  =  total  Compton  cross  section  per  electron, 

Z  =  atomic  number. 

More  useful  macroscopic  variations  of  the  total  gamma- 
ray  cross  section  are  the  linear  attenuation  coefficient 
fi,  given  in  units  of  cm'*  and  defined  as 

pA'o 

M  =  — .  (2.126) 

where  p  is  the  density  (g/cm^),  A^o  is  Avogadro’s 
number  =  6.023  X  10*^  nuclei/g-atom,  and  the  mass 
attenuation  coefficient  is  designated  by  p/p  in  units  of 
em’/g. 

The  gamma-ray  energy  deposition  coefficient  is 
usually  taken  as 

pTVo  / 

^1,=— (°PF.  +  Opp  +  ^OtO.  (2.127) 
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where  is  the  Compton  energy-deposition  cross 
section  per  electron.  This  relation  presumes  that  the 
photoelectric  and  pair-production  reactions  are  pure 
absorptions;  other  variations  are  sometimes  used  (see 
Section  2.9.3).  These  cross  sections  are  often  reported 
as  mass  energy  absorption  coefficients,  tijp.  Extensive 


tables  for  the  cross  sections  of  these  gamma-ray 
interactions,  as  well  as  others,  have  been  prepared  by 
Plechaty  and  Terrall^*  for  elements  1  through  100  and 
gamma-ray  energies  ranging  from  0.001  to  100  MeV, 
and  later  data  for  the  energies  1  keV  to  1  MeV  have 
been  published  by  McMaster  et 


2.8.  Gamma-Ray  Production  Processes 


It  was  pointed  out  in  Section  2.4.2  that  the  electro¬ 
magnetic  radiations  of  primary  concern  in  shielding  are 
the  nuclear  transition  gamma  rays.  G;unma  rays  pro¬ 
duced  by  a  weapons  burst  are  due  to  a  number  of 
sources,  with  only  one,  tlie  fission  reaction,  a  direct 
result  o*'  the  detonation.  All  other  sources  are  produced 
either  by  the  interaction  of  fission  or  fusion  neutrons 
with  various  materials  or  by  the  radioactive  decay  of 
the  fission  products.  The  important  sources  are  listed 
below. 


2.8.1.  FISSION  REACTION 

The  fission  reaction  utilized  in  weapons  detonations  is 
always  accompanied  by  the  emission  of  prompt  gamma 
lays,  whose  energies  total  about  8  MeV  per  fission. 
These  prompt  gamma  rays,  like  prompt  neutrons,  are 
released  whenever  the  compound  nucleus  formed  by  an 
incident  neutron  and  a  target  nucleus  splits  into  two  (or 
more)  fission  fragments  (see  Section  2.6.1).  The  energy 
spectrum  of  such  gamma  rays  is  known  reasonably  well 
only  for  ^  ^  ’  U,  measurements  having  been  made  for  this 
isotope  by  Motz^  *  and  by  Maienschein  et  al.  ^  ^ 

2.8.2.  NEUTRON  INELASTIC  SCATTERING 

In  the  neutron  inelastic-scattering  process  described  in 
Section  2.S.2  the  energy  of  the  target  nucleus  is 
elevated  to  an  excited  state,  from  which  it  returns  to 
the  ground  state  by  the  emission  of  gamma  rays.  In 
weapons  detonations  the  target  nuclei  can  be  weapons 
materials  or  environmental  materials,  such  as  those 
comprising  the  atmosphere  and  the  ground.  Investiga¬ 
tions  of  the  production  of  such  gamma  rays  have  been 
carried  out  by  Dickens  and  Perey^*  and  by  Orphan  and 
Hoot.^^  (See  also  Chapter  6  of  this  Handbook.) 

2.8.3.  NEUTRON  CAPTURE 

The  neutron  radiative-capture  process,  also  described 
in  Section  2.5.2,  results  in  the  de-excitation  of  the 
compound  nucleus  via  the  emission  of  gamma  rays.  And 
again,  the  target  nuclei  can  be  weapons  materials  or 
environmental  materials  such  as  those  comprising  the 


atmosphere  and  the  ground.  The  gamma  rays  emitted 
following  radiative  capture  can  be  given  off  immedi¬ 
ately,  in  which  case  they  are  called  capture  gamma  rays, 
or  they  can  be  emitted  at  later  times,  in  which  case  they 
are  called  activation  gamma  rays.  Activated  materials 
can  contribute  to  fallout  or  can  prevent  access  to 
regions  in  the  vicinity  of  a  detonation. 

The  radiative  capture  process  is  most  important  for 
low-energy  neutrons.  I'arly  data  on  capture  gamma  rays 
produced  by  neutron  captures  in  various  nuclei  were 
published  by  Hartholomew  and  ('ampion,''’  Keller  ct 
and  Iroubetzkoy  and  Goldstein,^ “  and  later 
information  has  been  presented  by  Maerker  and 
Muckenthaler,"  Uartholomew  a  and  (Jroshevtv 
long-term  emitters  of  activation 
gamma  rays  have  been  published  by  (Tucker  and 
Conners^*  and  by  (Tucker  and  Wong.^*’ 

2.8.4,  FISSION-PRODUCT  DECAY 

The  fission  fragments  produced  in  the  fission  process 
undergo  radioactive  decay  through  the  emission  of 
gamma  (and  beta)  radiations.  These  gamma  rays  consti¬ 
tute  an  important  source  from  the  moment  of  detona¬ 
tion  to  long  periods  following  the  detonation,  .Studies 
of  the  early-time  (<  I  min)  fission-product  gamma  rays 
have  been  made  by  Engle  and  Fisher.”  Data  on  the 
total  radiative  decay  of  fission  products  fur  limes  up  to 
70  yeais  have  been  reported  by  Crocker  and  l  uo'er."' 
Tlie  total  ganmia-ray  energy  released  throiigli  fission- 
product  decay  is  about  7  MeV  per  fission. 


2.8.5.  BRE.MSSTRAHLUNG 

Electromagnetic  radiation,  called  bremsstrahlung,  is 
emitted  when  a  free  charged  particle  experiences  an 
acceleration  or  deceleration  by  the  electrostatic  field  of 
a  nucleus.  The  total  bremsstrahlung  produced  varies  as 
the  square  of  the  atomic  number  of  the  recoil  nucleus 
and  inversely  witli  the  square  of  the  mass  of  the 
incident  particle.  Bremsstrahlung  is  characterized  by  a 
continuous  spectrum,  with  the  energy  emitted  by  a 
particular  particle  ranging  from  zero  up  to  its  total 
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kinetic  energy.  Because  of  the  strong  dependence  on 
the  incident  particle’s  mass,  charged  particles  other  than 
electrons  can  be  completely  neglected  insofar  as  the 
production  of  bremsstrahlung  is  concerned. 

Large  numbers  of  electrons  are  produced  as  second¬ 
aries  during  the  attenuation  of  gamma  rays;  for  ex¬ 


ample,  photoelectrons,  Compton  recoils,  and  the 
electron-positron  pairs.  All  these  electrons  will  produce 
numerous  bremsstrahlung  radiations;  however,  their 
energies  will  be  relatively  low  and  for  all  present 
shielding  calculations  they  are  assumed  to  be  absorbed 
at  the  point  of  formation. 


2.9.  Quantities  Used  to  Describe  Responses  to  Radiation 


Although  calculations  of  the  interactions  of  radiation 
with  matter  deal  with  the  fluxes  and  currents  described 
in  Section  2.1,  such  quantities  are  not  easily  measured 
in  the  laboratory,  especially  as  a  function  of  energy  and 
direction.  However,  the  property  of  a  radiation  field 
called  dose,  which  is  related  to  the  energy  deposited  in 
a  medium  by  the  radiation  interacting  within  it,  is 
determinable,  and  this  is  the  quantity  that  is  used 
to  relate  a  radiation  field  to  the  biological  damage  it 
causes  in  a  human  body. 

In  an  ideal  situation  a  given  dose  would  always 
produce  the  same  biological  effect  irrespective  of  the 
nature  and  energy  of  the  radiation  or  of  the  body  organ 
being  irradiated.  But  studies  by  radiobiologists  in  which 
doses  and  biological  damage  to  specific  organs  are 
correlated  for  known  radiation  fields  have  established 
that  nature  is  not  that  simple.  In  addition  to  being  a 
function  of  absorbed  energy,  biological  responses  are 
functions  of  the  irradiated  organ  and  of  the  type,  rate, 
and  energy  of  the  radiation.  From  this  discovery  has 
evolved  the  concept  of  a  relative  biological  effect 
(RBE),  which  is  a  weighting  factor  that  can  be  used  to 
compare  the  biological  effects  produced  by  the  same 
physical  dose  (same  amount  of  energy  deposited)  due 
to  radiation  of  a  different  type  and/or  energy.  When  the 
physical  dose  is  multiplied  by  the  RBE  it  becomes  a 
biological  dose. 

The  foregoing  gives  some  indication  of  the  problems 
associated  with  establishing  doses  for  shield  design  that 
are  related  to  biological  hazards.  The  shield  designer’s 
task  is  further  complicated  by  the  perturbation  of  the 
radiation  field  by  the  human  body,  an  effect  not 
usually  included  in  a  typical  shielding  calculation. 
Obviously  the  dose  received  at  a  particular  location 
within  the  body  is  not  the  same  as  the  dose  in  a  small 
detector  at  the  same  location  in  space  with  the  body 
absent.  It  is  apparent  that  shielding  studies  should 
include  anthropomorphic  phantoms  as  part  of  the 
shield  configuration,  a  theoretically  possible  but  usually 
impractical  consideration. 

As  an  alternative,  slab  or  cylindrical  phantoms  with 
compositions  resembling  that  of  the  human  body  have 
been  used,  and  the  doses  have  been  calculated  or 


measured  as  a  function  of  depth  in  the  phantom  for  a 
given  incident  radiation  field.  The  results  are  then  used 
as  response  functions  to  relate  an  unperturbed  radiation 
field  (usually  called  a  free  field)  to  the  dose  in  a  human 
body  if  it  had  been  present. 

It  will  be  apparent  in  the  following  discussion  on  the 
various  quantities  used  to  define  physical  and  biological 
doses  that  confusion  has  developed  in  the  terminology 
and  in  the  definition  of  units,  both  because  of  the 
burgeoning  nuclear  science  and  technology  and  because 
radiobiologists  and  sliield  designers  basically  differ  in 
their  viewpoints.  The  International  Commission  on 
Radiation  Cnits  and  Measurements*  recognized  this 
confusion  and  in  an  effort  to  mitigate  it  recom¬ 
mended*’^'  a  consistent  set  of  definitions  and  units. 
While  shield  designers  have  accepted  most  of  their 
recommendations,  they  have  continued  to  use  terms 
not  included  in  the  ICRU  recommendations  because 
these  terms  are  so  ingrained  in  llie  shielding  field  and 
because  the  ICRU  did  not  include  all  the  concepts 
needed  in  shield  design.  In  tJie  descriptions  given  below 
the  ICRU  recommendations  and  the  traditional  view¬ 
points  arc  contrasted. 

In  addition  to  descriptions  of  tlie  quantities  them¬ 
selves,  this  section  includes  data  for  some  of  the 
quantities  considered  most  pertinent  fur  shield  design. 


2.9.1.  ABSORBED  DOSE 

Absorbed  dose  is  the  energy  imparted  by  radiation  to 
a  unit  mass  of  matter  and  as  such  is  a  physical  quantity 
as  opposed  to  a  biological  effect.  A  fonnal  definition 
for  the  energy  imparted,  I'p,  is:  The  energy  imparted 


*Thc  International  Commission  on  Radiation  Units  and 
Measurements  has  as  its  principal  objective  the  development  of 
internationally  acceptable  recommendations  regarding  quanti¬ 
ties  and  units  of  radiation  and  radioactivity,  procedures  suitable 
for  the  measurement  and  application  of  these  quantities,  and 
physical  data  needed  in  the  application  of  these  procedures. 
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by  ionizing  radiation*  to  the  matter  in  a  volume  is  the 
difference  between  the  sum  of  the  energies  of  all  the 
directly  and  indirectly  ionizing  particles  which  have 
entered  the  volume  and  the  sum  of  the  energies  of  all 
those  which  have  left  it,  plus  the  energy  equivalent  (Q) 
of  any  decrease  in  rest  nmss  that  took  place  in  nuclear 
or  elementary  particle  reactions  within  the  volume.  For 
all  practical  purposes  the  energy  imparted  is  usually 
equal  to  the  heating  effect,  but  in  some  cases  part  of 
the  energy  deposited  may  result  in  changes  in  inter¬ 
atomic  bond  energies. 

As  proposed  by  the  ICRU,'  the  energy  imparted  can 
be  expressed  as 

where 

2]f-in  =  the  sum  of  the  energies  (excluding  rest 
energies)  of  all  those  directly  and  indirectly 
ionizing  particles  which  have  entered  the 
volume, 

2] =  the  sum  of  the  energies  (excluding  rest 
energies)  of  all  those  directly  and  indirectly 
ionizing  particles  which  have  left  the  volume, 

the  sum  of  all  the  energies  released,  minus 
the  sum  of  all  Uie  energies  expended,  in  any 
nuclear  reactions,  transformations  and  ele¬ 
mentary  particle  processes  which  have 
occurred  within  tJie  volume. 


The  special  unit  of  absorbed  dose  is  the  rad,  which  is 

1  rad  =  ^J/kg=  l(X)ergs/g.  (2.130) 

The  absorbed  dose  rate  is  the  quotient  of  AD  divided 
by  At,  where  AD  is  the  increment  in  absorbed  dose  in 
time  Af,  i.e.. 

Absorbed  dose  rate  =  — .  (2.131) 

At 

A  special  unit  of  absorbed  dose  rate  is  any  quotient  of 
the  rad  or  its  multiple  or  submultiple  by  a  suitable  unit 
of  time  (rads/sec,  mrads/hr,  etc.). 

Energy  can  be  imparted  to  a  volume  of  matter  by 
ntany  different  particle  reactions,  and  particles  having 
the  same  initial  energy  do  not  necessarily  deposit  the 
same  amount  of  energy.  This,  of  course,  is  because  the 
energy  deposition  depends  not  only  on  the  initial 
energy  but  also  on  the  type  of  radiation  and  the  kinds 
of  interactions  tliat  occur.  Schematic  representations  of 
energy  imparted  to  the  matter  contained  within  a 
volume  element  AV  for  two  particles  having  initial 
energies  /-.o  and  El  are  shown  in  Fig.  2.16  (superscripts 
1  and  2  identify  the  reactions  produced  by  particles  I 
and  2  respectively).  The  corresponding  equations  for 
the  energy  imparted  are 

Al-l^^El  F;‘  E\  Q\  (2.132fl) 


The  absorbed  dose  ID)  is  the  quotient  of  AE^^  by 
AM,  where  is  the  energy  imparted  by  ionizing 

radiation  to  the  matter  in  a  volume  element  and  AM  is 
the  mass  of  the  matter  in  tliat  volume  element: 


*loni/ing  radiation  is  a  radiation  consisting  of  directly  or 
indirectly  ionizing  particles  or  a  mixture  of  both.  Directly 
ionizing  particles  are  charged  particles  (electrons,  protons,  alpha 
particles,  etc.)  having  sufficient  kinetic  energy  to  produce 
ionization  by  collision.  Indirectly  ionizing  particles  are  un¬ 
charged  particles  (neutrons,  gamma  rays,  etc,)  which  can 
liberate  directly  ionizing  particles  or  can  initiate  a  nuclear 
transformation. 

^The  notation  AEp  implies  that  the  volume  clement  ^il' 
associated  with  the  clement  of  mass  XM  be  of  an  appropriate 
size  such  that  the  limiting  process  D  =  yield  a 

meaningful  estimate  of  the  absorbed  dose.  The  difficulties  arise 
because  the  interactions  between  radiations  and  atoms  are 
macroscopic  and  the  limiting  process  cannot  proceed  to  the 
usual  mathematical  limit  -*0. 


Al-;j,  =  El  El^Q\  (2.132/1) 

where  AEp  is  the  energy  imparted  to  the  matter  within 
the  volume  element  AV,  Q  is  the  energy  equivalent 
of  any  change  in  the  rest  mass  due  to  nuclear  or 
elementary  particle  reactions  within  tlie  volume,  and 
subscripts  0  and  7  refer  to  the  type  of  particle  leaving 
the  volume.  In  reaction  I,  an  incident  gamma  ray 
undergoes  a  Compton  scattering  within  the  volume 
element  AV.  Q'  is  the  binding  energy  of  the  Compton 
electron  which  is  usually  of  negligible  magnitude.  The 
Compton  electron  loses  some  of  its  energy  through 
ionization  within  tlie  volume  and  then  departs  with  the 
energy  E^.  The  scattered  (degraded)  gamma  ray  also 
leaves  the  volume.  In  the  second  reaction,  a  neutron 
undergoes  radiative  capture  and  a  gamma  ray  is  pro¬ 
duced  which  leaves  the  volume.  is  the  binding 
energy  associated  witli  the  neutron  capture  which 
appears  as  excitation  energy.  Note  that  no  secondary 
collisions  within  the  volume  element  were  included  in 
this  schematic  representation.  This  was  purposely  done 
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Fig.  2.16.  Schematic  Representation  of  Relationship  Between 
Absorbed  Dose  (Energy  Deposition)  and  Kerma. 


to  illustrate  that  the  volume  element  AV  should  be 
large  enough  for  many  primary  interactions  to  occur 
but  sufficiently  small  so  that  a  primary  particle  and/or 
secondary  neutrons  or  gamma  rays  usually  do  not  suffer 
collisions  subsequent  to  the  iriitial  collision  by  the 
incident  particle  within  the  volume  element  AV. 

The  average  energy  imparted  to  the  matter  within  the 
volume  element  AK  by  all  reactions  of  particles  of  a 
particular  type  and  energy  can  be  written  as 


2  Ri(f)Ei{E) 


(2.133) 


where 


Ep(E)  =  the  average  energy  imparted  in  ergs  or  MeV, 

Rj(E)  =  rate  of  the  /th  type  of  reaction  for  particles 
of  energy  E  within  the  volume  element  AV, 

Ej(E)  =  energy  imparted,  in  ergs  or  MeV,  associated 
with  a  particle  of  energy  E  undergoing 
reaction  /  within  the  volume  element  AV. 


2.9.2.  FIRST-COLLISION  DOSE  AND  KERMA 

The  quantity  absorbed  dose  discussed  in  the  previous 
section  is  a  physical  variable  that  is  closely  related  to 


the  potential  hazard  of  a  radiation  field  but  for  a  real 
situation  is  sometimes  very  difficult  to  calculate  accu¬ 
rately  or  to  relate  to  the  response  in  a  human  being. 
Consequently,  other  concepts  have  often  been  used. 

One  of  tlie  concepts  used  is  the  first-collision  dose 
(also  called  tlie  single-collision  dose).  This  quantity  has 
been  subject  to  several  interpretations  that  differ  from 
one  another  in  subtle  ways.  One  of  the  most  compre¬ 
hensive  reviews  of  the  meanings  of  first-collision  dose 
has  been  written  by  Attix,^®  who  says  that  the 
first-collision  dose  due  to  a  radiation  field  (neutrons 
and/or  gamma  rays)  can  be  interpreted  in  each  of  the 
following  three  ways: 

1.  The  first-collision  dose  is  the  absorbed  dose 
contributed  by  all  particles  incident  on  an  isolated  small 
mass  (unless  the  term  is  explicitly  modified  to  include 
only  one  component,  e.g.,  neutrons)  in  which  charged- 
particle  equilibrium*  exisis. 

2.  The  fiist-collision  dose  is  the  absorbed  dose 
contributed  by  primary  radiation  only  in  a  small  mass 
that  can  be  located  anywhere  and  in  which  charged 
particle  equilibrium  exists. 

3.  The  first-collision  dose  is  the  energy  transferred 
from  uncharged  incident  radiation,  regardless  of  origin, 
to  the  charged  particles  formed  in  a  limitingly  small 
probe  that  may  be  located  anywhere.  Charged-particle 
equilibrium  in  the  probe  is  not  a  requirement. 

The  small  mass  specified  in  interpretations  I  and  2 
means  that  the  probe  is  small  enough  to  leave  the 
radiation  field  unperturbed  and  the  probability  is 
negligible  that  the  incident  particles  will  interact  with 
the  probe  more  than  once  or  that  secondary  gamma 
rays+  produced  within  the  mass  will  be  absorbed 
within  it.  The  small  probe  specified  in  interpretation  3 
has  the  same  requirements  but  it  can  be  smaller  tlian 
that  required  for  the  other  interpretations  because 
charged-particle  equilibrium  is  not  necessary. 

It  can  be  stated  at  this  point  that  interpretation  2  is  a 
special  and  seldom  used  definition  tliat  will  not  be 
elaborated  on  here.  Interpretation  1 ,  on  the  other  hand, 
is  a  widely  used  definition  -  particularly  by  experi¬ 
mentalists  —  and  is  often  referred  to  as  the  free-field 
dose,  the  air  dose,  the  firee-air  dose,  etc.  But  in  an 


*Charged-paiticle  equilibrium  may  be  viewed  as  that  condi¬ 
tion  when  on  the  .sverage  as  many  charged  particles  such  as 
electrons  enter  the  volume  element  AF  as  leak  out,  thereby 
resulting  in  an  essentially  zero  net  transfer  of  energy  by  the 
flow  of  electrons. 

^  As  used  here,  secondary  gamma  nys  refer  to  those  gamma 
rays  produced  by  the  interaction  of  the  incident  radiation  (e.g., 
capture  gamma  rays  and  inelastic-scattering  gamma  rays). 
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attempt  to  alleviate  confusion,  the  ICRU  chose  inter¬ 
pretation  3  to  use  as  a  basis  for  an  exact  definition  and 
called  the  quantity  so  defined  the  kerma  (A:inetic  energy 
released  in  material).  In  fact,  the  ICRU  encourages  the 
exclusive  use  of  the  term  kerma'  with  the  implication 
that  the  term  first-collision  dose  and  its  other  interpre¬ 
tations  should  be  relegated  to  the  archives. 

The  definition  of  kerma,  K,  as  recommended  by  the 
ICRU'  i!, 

(2.134) 

where  =  sum  of  the  initial  kinetic  energies  (ergs) 
of  all  the  charged  particles  liberated  by  indirectly 
ionizing  particles  in  a  volume  element  of  material  AV, 
and  AM  =  mass  (g)  of  the  material  contained  in  the 
volume  element  AV. 

The  kerim  rate  is  the  quotient  of  AK  by  At,  where 
AK  is  the  increment  in  kerma  in  time  At: 

Kenna  rate  =  AA'/A/.  (2.135) 

Since  is  the  sum  of  the  initial  kinetic  energies  of 
all  charged  particles  liberated  by  indirectly  ionizing 
particles,  it  includes  not  only  the  kinetic  energy  these 
charged  particles  expend  in  collisions  but  also  the 
bremsstrahlung  they  radiate  and  the  energy  of  any 
charged  particles  that  arc  produced  in  secondary 
processes. 

The  kenna  or  kenna  rate  for  a  specified  material  in 
free  space  or  at  a  point  inside  a  different  material  will 
be  that  which  would  be  obtained  if  a  small  quantity  of 
the  specified  material  were  placed  at  the  point  of 
interest. 

With  the  definition  of  kenna  in  mind,  tlie  reactions 
shown  in  Fig.  2.16  are  interpreted  in  tenns  of  the  initial 
kinetic  energy  imparted  to  charged  particles; 

Ah:l;=k:'o  0'.  (2.1 36a) 

=/;■?  (2.l36ii) 

It  is  noted  tliat  the  expression  for  At'j^  and  the  ex¬ 
pression  given  in  Hq.  2.132A  for  AEJj  are  the  same, 
while  the  expression  for  Mj^-  differs  from  that  given  in 
Eq.  2.132a  for  This  is  because  the  Compton 

electron  produced  by  reaction  1  did  not  deposit  all  of 
its  energy  within  the  volume  element  and  there  was  no 
compensating  in-leakage  of  an  electron;  i.e.,  charged- 
particle  equilibrium  did  not  exist. 

Auxier*"  strongly  criticizes  the  ICRU  for  not  for¬ 
mally  recognizing  first-collision  dose  as  given  by  inter¬ 


pretation  1.  His  objections  are  based  on  the  fact  that 
the  first-collision  doses  obtained  from  measurements 
with  ionization  chambers,  proportional  counters,  etc., 
are  absorbed  (imparted)  energy  and  not  transferred 
kinetic  energy.  Therefore,  since  measured  doses  cannot 
be  interpreted  as  kerma,  except  by  deduction,  some 
dosimetrists  are  insisting  on  retaining  the  term  first- 
coUision  dose.  On  the  other  hand,  if  charged-particle 
equilibrium  exists,  which  it  may  or  may  not,  and  if 
bremsstrahlung  radiation  is  ignored,  then  first-collision 
dose  and  kerma  are  equivalent.  It  can  be  seen  that,  in 
bulk  shields  or  large  tissue  masses,  first-collision  dose 
and  kerma  have  nearly  the  same  magnitude  but  that  in 
thin  layers,  such  as  clothing  or  skin,  they  can  be  quite 
different. 

The  early  calculations  of  fust-collision  dose  per  unit 
fiuence  were  performed  by  Snyder^'  for  mono- 
energetic  neutrons  incident  on  a  four-element  tissue 
model.  Although  the  term  kerma  had  not  yet  been 
introduced,  the  doses  obtained  by  Snyder  are  those 
described  by  interpretation  3  above  and  thus  are  kerma 
per  unit  fiuence.  His  calculations  considered  only 
neutron  captures  and  elastic  scatterings,  with  the  latter 
assumed  to  take  place  isotropically  in  the  center-of- 
mass  system.  Later,  Henderson**  made  similar  calcula¬ 
tions  for  n.utrons  for  a  four-element  tissue  model  and 
for  ethylene.  He  reported  his  results  as  absorbed  dose 
(in  rads/hr  per  unit  fiux),  but  they  are  actually  kerma 
rate  per  unit  flux. 

The  most  recent  and  most  comprehensive  calculations 
of  neutron  kemia  factors  (kemia  per  unit  fiuence)  were 
made  by  Ritts  et  al.*^  They  included  essentially  all 
reactions,  considered  anisotropic  scattering  in  those 
reactions  for  which  angular  cross-section  data  were 
avaUable,  and  in  all  cases  used  the  most  recent 
cross-section  data.  They  used  several  models,  one  of 
which  was  an  11 -element  “standard  man’’  model  in 
which  all  tissue  organs  and  bone  were  homogenized  and 
the  other  models  were  for  specific  body  organs.  The 
elemental  compositions  used  in  the  calculations  are 
shown  in  Table  2.8,  and  all  the  reactions  considered  are 
listed  in  Table  2.9. 

The  calculations  by  Ritts  et  al.  were  made  for  816 
discrete  neutron  energies  corresponding  to  the  super- 
group-subgroup  structure  of  the  05R  Monte  Carlo 
code**  over  an  energy  range  0.023  eV  to  19.2  MeV. 
The  kerma  factor  for  a  particular  irradiated  material 
was  found  by  summing  the  averages  of  the  kinetic 
energies  imparted  to  the  struck  nuclei  and  the  energies 
associated  with  charged  particles  that  were  emitted. 
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Table  2.9.  Reactions  Considered  in  Calculations  of  Neutron  Kerma  Calculations^ 


Type  of  Reaction 

H 

0 

c 

N 

Ca 

p 

s 

K 

Na 

Cl 

Mg 

Radiative  capture 

X 

X 

X 

X 

X 

x» 

X 

X 

Elastic  scattering 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

Inelastic  scattering 

X 

X 

X 

X 

X 

X 

X 

X 

X 

(n,p)  reaction 

X*- 

X 

X* 

X 

X 

X 

(n.Q)  reaction 

X 

X 

X 

x* 

X 

x* 

X 

(n,n'3a)  reaction 
(fi,2n)  reaction 

X 

X*^ 

x'^ 

(n,d)  reaction 

X 

X 

(n,r)  reaction 

X 

(/i,2o)  reaction 
(n.n'p)  reaction 

X 

X 

X 

X 

(n,n'a)  reaction 

X 

‘'l-rom  tef.  43. 

^Includes  beta  decay  after  reaction. 
^Includes  positron  decay. 


Therefore  the  kenna  per  unit  fluence  niay  be  expressed 
as 

=  LI  CA', .(/•-■)  (2.137) 

where 

Nj  =  number  density  of  element  /  (atoins/g), 

o,y  =  microscopic  cross  section  of  reaction  /  for 
element  »(cm^/atom), 

h'jj-  =  average  energy  deposited  by  reaction  j  for 
element  /  (MeV), 

C=  1.602  X  I  O'*  (ergs/MeV). 

Since  the  kemia  values  are  given  per  unit  fluence, 
they  are  fluence-lo-kerma  conversion  factors  (usually 
referred  to  simply  as  kerma  factors).  Such  kerma 
factors  for  the  standard  man  model  are  plotted  in  Fig. 
2.17  and  tabulated  in  Table  2A.1  in  Appendix  2A, 
along  with  kerma  factors  for  the  individual  elements 
making  the  greatest  contributions  to  the  total  factor. 
Kerma  factors  for  specific  body  components  (lung, 
muscle,  bone,  brain,  and  red  marrow)  are  also  tabulated 
in  Appendix  2A,  as  well  as  unweighted  kerma  factors 
for  the  important  elements  in  tlie  standard  man.  With 
these  values  a  total  kemia  factor  for  any  arbitrary 
composition  of  these  materials  can  be  easily  calculated. 

Kenna  rate  factors  for  monoenergetic  gamma  rays 
incident  on  carbon,  air,  and  a  four-element  tissue  model 
were  calculated  by  Henderson'*^  and  are  plotted  in  Fig. 


2.18.  (Henderson’s  gamma-ray  results,  like  his  neutron 
results,  were  actually  reported  as  absorbed  dose,  the 
term  kerma  not  having  been  introduced  yet.)  Prior  to 
Henderson’s  work,  exposure  factors  for  air  published  by 
Rockwell*®  were  generally  used  as  the  equivalent 
flux-to-dose  rate  conversion  factors  for  tissue. 


:.9.3  EXPOSURE 

Exposure  is  a  term  that  should  be  used  only  for 
gamma  rays.  Formerly  referred  to  as  exposure  dose,  it 
came  into  common  usage  after  the  problem  of  speci¬ 
fying  biological  dose  associated  with  X  rays  was  first 
encountered.  As  recommended  by  the  ICRU,  exposure 
(20  describes  the  deposition  of  energy  in  air  by 
electromagnetic  radiation  and  is  defined  by 

X  =  AQIAM,  (2.138) 

where  AG  =  sum  of  the  electrical  charges  of  all  the  ions 
of  one  sign  produced  in  air  when  all  electrons  liberated 
by  photons  in  a  volume  of  air  AV  are  completely 
stopped,*  and  AM  =  mass  of  the  material  contained  in 
the  volume  element  AV.  Exposure  rate  is  given  by 


*Thc  ionization  that  would  be  produced  by  the  brems- 
strahJung  associated  with  secondary  electrons  is  not  included  in 
^Q.  Lxeept  for  this  small  difference,  which  is  insignificant  for 
photon  energies  less  than  IS  MeV,  exposure  is  equivalent  to 
kerma. 
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(From  Rittset  al.,  ref.  43.) 
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Fig.  2.18.  Gamma-Ray  Kerma-Rate  Factors  for  Air,  Carbon,  and  Tissue.  (From  Henderson,  ref.  42.) 
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AA7A/,  where  AX  is  the  increment ‘of  exposure  during 
the  increment  of  time  At. 

The  unit  of  exposure  is  the  roentgen  (r).  It  is  the 
quantity  of  X  or  gamma  radiation  that  produces,  in  air, 
ions  carrying  1  electrostatic  unit  of  charge  per 
0.001793  g  of  dry  air  (or  2.58  X  10"^  coulombs/kg).  In 
terms  of  ergs,  the  roentgen  is  equal  to  87.7  ergs  per 
gram  of  air  or  96  ergs  per  gram  of  tissue.  For  gamma 
rays  above  3  MeV  in  energy,  the  range  of  secondary 
electrons  in  low-atomic-number  materials  becomes  com¬ 
parable  with  the  relaxation  length  of  the  gamma  rays. 
Consequently,  the  ionization  produced  in  a  small 
volume  is  no  longer  a  sole  measure  of  the  intensity  of 
the  radiation  at  that  point.  The  ICRU  does  not 
recommend  the  use  of  the  roentgen  above  3  MeV,  but 
in  practice,  it  is  still  used  above  3  MeV  with  instruments 
calibrated  in  terms  of  energy  absorption  in  air. 

Since  I  rad  is  equal  to  100  ergs  per  gram  of  irradiated 
material  (see  Section  2.9.1),  the  roentgen  and  the  rad 
are  frequently  interchanged  when  tissue  exposure  is 
referred  to,  the  difference  of  4%  being  no  real  conse¬ 
quence  in  shielding  calculations.  Strictly  speaking, 
however,  the  rad  should  be  reserved  for  describing  an 
absorbed  dose. 

Away  from  boundaries,  the  distinction  between 
kerma  and  exposure  loses  relevance  if  only  the  three 
principal  processes  (Compton  scattering,  photoelectric 
effect,  and  pair  production)  are  considered  and  if  the 
photoelectric  effect  and  pair  production  are  assumed  to 
be  absorptions.  For  this  simple  and  widely  used  model* 
the  fluence-to-exposure  conversion  factor  is  propor¬ 
tional  to  the  product  of  the  photon  energy  and  the 
energy-deposition  coefficient,  in  which  case  exposure, 
absorbed  dose,  first-collision  dose,  and  kemia  in  air  are 
all  equal  in  magnitude. 

In  some  calculations,  however,  a  sliglitly  more  com¬ 
plex  model  is  used  in  which  only  the  photoelectric 
process  is  treated  as  an  absorption.  As  in  the  above 
model,  the  total  kinetic  energy  of  the  electron  pair 
created  in  the  pair-production  process  is  assumed  to  be 
absorbed  at  tlie  point  of  their  emission.  However,  the 
two  0.51 -MeV  gamma  rays  produced  by  the  annihila¬ 
tion  of  the  positron  are  treated  as  scattered  ganuna 
rays,  and  pair  production  is  viewed  as  a  scattering 
process.  This  is  described  in  terms  of  a  modified 
pair-production  cross  section  which  is  given  by 

|Upp^^=trpp  (^1  £>  1.022  MeV.  (2.1.39) 


*l-or  this  und  the  subsequent  model  it  is  assumed  that  the 
bremsstrahlung  radiation  that  would  be  produced  by  the 
secondary  electrons  (and  positrons)  is  immediately  absorbed. 


where  £  is  the  energy  of  the  incident  photon  (MeV). 


2.9.4  RBE  DOSE;  DOSE  EQUIVALENT 

Whenever  the  biological  response  of  a  human  organ  to 
radiation  exposure  is  of  concern,  merely  knowing  the 
absorbed  dose  (energy  deposited)  is  not  enough.  Bio¬ 
logical  responses  vary  both  with  the  nature  and  energy 
of  the  radiation  and  with  the  composition  and  function 
of  the  irradiated  organs.  Thus  when  an  organ  is  exposed 
to  a  mixed  radiation  field  or  to  a  field  of  radiation 
comprised  of  one  type  of  particle  with  various  energies, 
the  energy  deposited  by  each  type  of  particle  of  a  given 
energy  must  be  weighted  by  some  factor  before  the 
total  biological  dose  received  by  the  organ  can  be  de¬ 
termined  by  a  summing  of  the  individual  contributions. 

For  a  specific  biological  effect  in  a  particular  mass  the 
weighting  factor  is  mainly  a  function  of  the  linear  rate 
of  energy  transfer  (/.£T)  to  the  system  by  charged 
particles  set  in  motion  by  the  interactions  of  incident 
radiation.  The  LET  of  a  material  is  related  to  the  linear 
stopping  power  but .  the  concepts  are  different.  The 
linear  stopping  power  is  the  average  energy  loss  per  unit 
path  length  by  a  charged  particle  in  traversing  a  medium 
regardless  of  where  the  energy  is  absorbed.  The  LET, 
however,  refers  to  energy  imparted  within  a  limited 
volume.  These  two  quantities  are  equal  only  in  the 
special  case  when  the  LET  includes  the  absorption  of  all 
secondary  particles,  in  which  case  it  is  called  LET.. 
These  quantities  are  described  more  completely  in  ref. 
21(p.51). 

The  LET-dependent  weighting  factor  is  called  the 
relative  biological  effectiveness  {RBE)  and  is  defined  as 

rads  of  standard  radiation  producing 
_  a  given  biological  effect  (2  140) 

rads  of  another  type  of  radiation 
producing  the  same  effect 

The  standard  radiation  referred  to  is  X  rays  of  250-keV 
energy,  and  thus  by  definition  the  RBE  for  250-keV  X 
rays  is  1 . 

When  a  dose  given  in  rads  is  weighted  by  an  RBE 
value,  tlie  resulting  dose,  sometimes  referred  to  as  the 
RBE  dose,  is  given  in  rems,  a  unit  derived  from  the  term 
roentgen  equivalent  man: 

RBE  dose  in  rems  =  dose  in  rads  X  RBE  .  (2.1.41) 

It  follows  that  for  250-keV  X  rays  the  rad  dose  and  the 
rem  dose  are  numerically  equal. 
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The  rem  is  quantitatively  defined  as;  ihe  absorbed 
dose  due  to  any  ionizing  rmiiation  that  has  the  same 
bioiogical  effectiveness  as  I  rad  of  X  rays  with  an 
average  specific  ionization  of  WO  ion  pairs  per  micron 
of  water  in  terms  of  its  air  equivalent.  It  turns  out  that 
X  rays  and  ganuna  rays  generally  do  not  exceed  tliis 
specific  ionization  (LET  <  3.5  keV/p),  which  is 
considered  to  be  the  boundary  between  low-LET  and 
high'LET  radiation  and  thus  to  be  the  limiting  condi¬ 
tion  for  RBE  *  1. 

In  the  past  RBE  has  been  used  in  the  fields  of  both 
radiobiology  and  radiation  protection,  but  this  gener¬ 
ated  concern  because  of  the  differences  in  application 
and,  to  sorr>e  extent,  in  concept.  Consequently  ICRU 
recommends  that  RBE  be  used  only  for  correlating 
radiobiological  experiments  and  that  a  new  term,  the 
quality  factor  (QF),  be  used  in  the  field  of  radiation 
protection.  Quality  factors  are  actually  those ‘values  of 
RBE  that  are  intended  to  embrace  all  effects  that  are 
hazardous  to  human  beings.  In  other  words,  QF  values 
are  not  related  to  specific  organs  of  the  body  as  are 
some  of  the  RBE  values.  V.'hen  a  QF  value  is  used  to 
weight  absorbed  dose,  the  resulting  dose  in  rems  is 
identified  as  the  dose  equivalent  (DE).  It  is  assumed 
that  dose  equivalents,  like  RBE  doses,  are  additive. 

Although  QF  is  defined  in  purely  physical  terms  (that 
is,  as  a  function  of  LET  only,  which  in  turn  is  a 
function  of  the  number  of  ion  pairs  produced  per 
centimeter  of  travel  by  charged  particles),  the  basis  of 
the  legislated  value  is  biological.  Table  2.10  shows  the 
recommendations^®’^’  of  the  International  Committee 
on  Radiological  Protection  (ICRP)  for  the  relation 
between  QF  and  both  LET  and  ion  pair  formation. 

The  ICRP  recommends  a  QF  of  1  for  X  rays,  gamma 
rays,  electrons,  and  positrons  of  any  specific  ioniza¬ 
tion;®®  QF’s  for  neutrons,  protons,  and  heavy  recoil 
nuclei  vary  with  energy.  The  ICRP  recommendations 
for  QF  for  neutrons  between  0.01  and  1000  MeV  are 
plotted  in  Fig.  2.19.  Tabulated  values  for  lower  energy 
neutrons  are  also  shown.  These  values  are  based  on  the 
calculated  results  of  Snyder  and  Neufeld®®  and  of 
Neary  and  Mulvey.®*  As  a  rough  and  conservative 
approximation,  a  QF  of  10  is  applicable  to  protons  and 
a  QF  of  20  can  be  used  for  heavy-recoil  nuclei. 

The  preceding  QF’s  are  assumed  to  be  applicable  for 
whole-body  irradiation.  When  specific  o.jans  are  being 
considered,  additional  modifying  factors  may  be 
needed.  For  example,  ICRP  specifies  that  when  the  lens 
of  the  eye  is  being  irradiated  a  modifying  factor  of  3 
should  be  used  if  the  QF  is  equal  to  or  greater  than  10 
and  a  modifying  factor  of  1  should  be  used  when  QF  = 

1 .  Strangely  enough  no  recommendations  are  given  for 


I  <  QF  <  10.  Presumably,  linear  interpolation  is  to  be 
used  between  these  limits. 

When  the  dose  in  an  organ  is  nonuniform,  a  dose 
distribution  factor  (DF)  should  also  be  applied.  These 
factois  cannot  be  established  with  certainty  and  only 
one  is  in  current  use.  It  is  the  relative  damage  factor 
applied  in  calculating  the  dose  equivalent  in  bone  from 
internal  radiation. 

In  general  the  shielding  community  has  accepted  the 
use  of  QF  and  DE,  as  is  reflected  by  the  latest  shielding 
literature.  From  a  practical  viewpoint  the  change 
merely  amounts  to  replacing  RBE  by  QF  and  calling  the 
product  of  the  quality  factor  and  absorbed  dose  the 
dose  equivalent  (DF)  instead  of  the  RBE  dose.  Conse¬ 
quently  no  confusion  should  exist  when  tlie  older 


Table  2.10.  QF  as  a  Function  of 
LET  and  Ion  Pair  Fonnation" 


QF 

Average  Specifle  Ionization 
(ion  pairs/ji*) 

Average  LET 
(keV/p*) 

1 

<100 

<3.5 

2 

200 

7.0 

5 

650 

23 

10 

1500 

53 

20 

2000 

175 

"From  refs.  46  and  47. 
*Of  water. 
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Fif.  2.19.  Quality  Factors  for  Neutrons.  (From  ref.  46.) 


53 


literature  is  consulted.  It  should  be  kept  in  mind, 
however,  that  QF  values  are  chosen  by  such  groups  as 
the  ICRP  and  thus  are  more  subject  to  change  than 
would  be  an  RBE  value  based  on  experimental  data. 

It  should  also  be  pointed  out  that  all  QF's  and  other 
modifying  factors  are  intended  solely  for  chronic 
exposure  to  low-level  radiation  Fields,  with  genetic 
damage  being  the  hazard  considered.  High  accidental 
exposures  must  be  assessed  on  the  basis  of  particular 
circumstances.  The  ICRP  has  no  recommendations  for 
such  exposures,  but  in  many  cases  the  absorbed  dose 
will  give  a  better  indication  of  the  biological  risks  than 
will  the  dose  equivalent.  In  other  words,  for  acute 
effects  due  to  massive  exposures,  the  QF  should  be 
taken  as  unity,  which  partially  results  because  energy 
may  be  wasted  from  the  standpoint  of  the  production 
of  biological  effects.  Consequently  for  military  applica¬ 
tions  during  war,  when  early  death  or  incapacitation  is 
the  effect  to  be  considered,  the  dose  in  rems  may  be 
taken  as  numerically  equal  to  the  absorbed  or  physical 
dose  until  better  information  becomes  available. 

2.9.S  MAXIMUM  ABSORBED  DOSE; 

MAXIMUM  DOSE  EQUIVALENT 

As  is  pointed  out  in  the  introduction  to  this  section, 
the  usual  shielding  calculation  provides  a  detailed 
description  of  the  unperturbed  radiation  field  (free 
field).  That  is,  perturbations  of  the  field  by  a  human 
body  are  not  considered.  This  greatly  simplifies  the 
analysis  without  seriously  compromising  the  overall 
accuracy.  But  the  results  must  somehow  be  related  to 
the  hazards  to  a  human  being  if  he  were  exposed  to 
such  a  radiation  field. 

A  correlation  between  the  unperturbed  radiation  field 
and  the  dose  in  a  human  body  has  been  accomplished 
through  the  use  of  slab  and  cylindrical  phantoms,  which 
have  dimensions  and  compositions  resembling  those  of 
the  human  body.  Calculations  (or  measurements)  of  the 
absorbed  dose  or  the  biological  dose  were  made  as  a 
function  of  depth  in  the  phantom  for  given  mono- 
energetic  neutrons  or  gamma  rays  incident  on  the 
phantom.  The  incident  angular  distribution  was  a 
monodirectional  beam  or  an  isotropic  flux.  (Calcula¬ 
tions  for  neutrons  include  the  contribution  to  the  dose 
by  secondary  gamma  rays.)  The  maximum  values  in  the 
depth-dose  distributions  were  identified  as  the  max¬ 
imum  absorbed  dose  and  maximum  dose  equivalent  for 
the  absorbed  (physical)  and  biological  doses  respec¬ 
tively.  Since  these  values  were  obtained  and  reported  on 
the  basis  of  a  unit  particle  entering  the  phantom,  they 
can  be  used  as  response  functions  to  relate  an  unper¬ 


turbed  field  of  mixed  radiation  to  the  maximum  dose 
that  would  be  received  by  some  part  of  the  body.  Use 
of  these  maximum  values  is  dictated  by  a  conservative 
design  philosophy  that  does  not  allow  the  permissible 
dose  to  be  exceeded  at  any  point  in  the  body. 

Neutron  Doiet  in  Phantomi 

The  first  dose  calculations  in  a  phantom  were 
performed  by  Snyder  and  Neufeld**  for  a  beam  of 
monoenergetic  neutrons  normally  incident  on  a  phan¬ 
tom  represented  by  a  slab  of  tissue.  The  slab,  infinite  in 
the  transverse  directions,  was  assumed  to  be  30  cm 
thick.  Snyder  and  Neufeld  determined  the  distribution 
of  the  absorbed  dose  through  the  slab,  finding  that  the 
maximum  dose  occurred  at  nearly  the  surface  or  within 
a  few  centimeters  of  the  surface.  In  these  calculations, 
only  neutron  captures  and  elastic  scatterings  were 
considered,  which  limited  the  contribution  by  second¬ 
ary  gamma  rays  to  capture  gamma  rays.  The  resulting 
maximum  absorbed  doses  and  biological  doses  as  a 
function  of  the  incident  neutron  energy  are  shown  in 
Table  2.11.  These  values  have  been  widely  used  in 
reactor  shield  design  to  convert  neutron  fluences  to 
doses  and  have  become  a  virtual  standard  since  these 
dose  equivalents  have  been  stipulated  for  use  by  the 
Federal  Register.*® 


Tabic  2.1 1.  Maximum  Abioibed  Dote  and  Maximum 
Doae  Equivalent  for  Monoenergetic  Neutrons 
Incident  on  a  Slab  Tissue  Phantom" 


Neutron 

Energy 

(MeV) 

Maximum 

Absorbed 

Dose® 

(mrad 

neutron”*  cm”^) 

Maximum 

Dose 

Equivalent® 

(mrem 

neutron”*  cm”*) 

Effective 

QF 

Thermal 

3.2(-7)" 

1.04(-6) 

3.25 

0.0001 

6.9(-7) 

1.39(-6) 

2.01 

0.005 

5.7(-7) 

1.22(-6) 

2.14 

0.02 

5.7(-7) 

2.35(-6) 

4.12 

0.1 

1.10(-6) 

8.3(-6) 

7.55 

0.5 

2.4(-6) 

2.30(-5) 

9.61 

1.0 

3.8(-6) 

3.80(-5) 

10.0 

2.5 

4.3{-6) 

3.41(-5) 

7.93 

5.0 

5.8(-6) 

3.80(-5) 

6.55 

7.5 

7.1(-6) 

4.16(-5) 

5.85 

10 

7.0(-6) 

4.16(-5) 

5.94 

"From  ref.  48. 

^Multiply  by  3600  to  convert  to  (miad/hT)/(neutron/cm*‘sec) 
or  to  (mrem/hr)/(neutton/cm*‘sec). 

"Read:  3.2  X  10’’,  etc. 
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In  similar  calculations,  perfunned  later  by  Auxicr  vl 
al,^^  the  inOnite  slab  was  replaced  by  a  Hnite 
cylindrical  phantom  that  more  nearly  mocked  up  a 
human  body.  A  schematic  of  tlie  phantom,  60  cm  high 
and  30  cm  in  diameter,  is  given  in  Fig.  2.20.  The 
maximum  doses  (averaged  fur  a  volume  element) 
occurred  in  the  outermost  volume  element  at  the 
midplane  on  the  side  facing  the  beam.  These  maximum 
doses  are  given  in  Table  2.12,  and  the  doses  fur  each 
volume  element  are  presented  in  Appendix  2B. 

The  effective  quality  factors  shown  in  Tables  2.1 1 
and  2.12  are  the  ratios  of  the  maximum  dose  equivalent 
to  the  maximum  absorbed  dose  for  a  given  incident 
neutron  energy  and  phantom  configuration.  The  magni¬ 
tude  of  the  effective  quality  factor  is  very  close  to  but 
usually  slightly  less  than  the  QF’s  recommended  by  the 
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Fig.  2.20.  Cylindiical  Phantom  Used  in  Neutron  Dose  Cat- 
culations.  (From  Auxier  et  al.,  ref.  51.) 


Table  2.12.  Maximum  Abturbed  Dose  and  Maximum 
Duse  Equivalent  for  Munoenergetiv  Neutrons 
Invideni  on  a  Cylindrical  Tissue  Phantom'' 


Neutron 

l^ncrgy 

(MeV) 

Maximum 

Absorbed 

Dose^ 

( mrad 

neutron"*  cm"*) 

Muximum 

Dose 

Equivalent*’ 

(mrem 

neutron"'  em"^) 

Effective 

QE 

Thermal 

4.68(  If 

l.l5(-6) 

2.46 

0.000001 

5.89(  7) 

1.34(-6) 

2.27 

0.00001 

5.l8(-7) 

1.21(-6) 

2.34 

0.0001 

4.45(-7) 

1.01(-6) 

2.27 

0.001 

4.32(-7) 

8.85(-6) 

2.04 

0.01 

4.34(-7) 

9.92(-6) 

2.29 

0.1 

8.02(-7) 

4.86(-6) 

6.06 

0.5 

l.81(-6) 

1.89(-5) 

10.4 

1 

3.01(-6) 

3.26(-5) 

10.8 

2.5 

3.99(-6) 

3.50(-5) 

8.77 

5 

5.72(-6) 

4.41(-5) 

7.70 

7 

5.70(-6) 

4.03(-5) 

7.07 

10 

7.25(-6) 

4.31(-5) 

5.94 

14 

8.31(-6) 

6.15(-5) 

7.40 

‘'i•'rom  ref.  5 1 . 

^Multiply  by  3600  to  convert  to  (mrad/hr)/(ncutron/cm*'scc) 
or  to  (mrcm/hr)/(neutron/cm*’scc). 

''Read;  4.6S  X  lO’’,  etc. 


ICRP  (see  Fig.  2.19).  The  differences  in  the  values 
would  be  expected  since  the  QF’s  recommended  by  the 
ICRP  consider  only  the  initial  collision  of  the  neutron 
whereas  the  effective  QF’s  represent  an  average  of  all 
collisions  experienced  by  the  neutron  within  the  phan¬ 
tom. 

The  Snyder-Neufeld  calculations  of  maximum  ab¬ 
sorbed  dose  have  recently  been  improved  upon  by  Ritts 
et  al.*^  They  increased  the  number  of  elements  in  the 
infinite-slab  phantom  from  4  to  II,  considered  a  much 
larger  number  of  neutron  source  energies,  allowed  for 
neutron  inelastic  scattering  (and  the  associated  ine¬ 
lastic-scattering  gamma  rays)  and  for  a  number  of  other 
noneiastic  events,  corrected  the  elastic  scattering  to 
linear  anisotropy,  used  more  accurate  secondary- 
gamma-ray  spectra,  used  the  latest  cross  sections  avail¬ 
able,  and  handled  the  neutrons  with  energies  below  1 
eV  more  accurately.  The  percentages  of  the  1 1  elements 
in  the  slab  were  the  same  as  those  used  for  the  standard 
man  kerma  calculations  (see  Table  2.8),  and  the 
reactions  considered  were  also  the  same  as  for  the 
kerma  calculations  (see  Table  2.9).  The  resulting  maxi¬ 
mum  absorbed  doses  are  presented  in  Fig.  2.21  (solid 
line)  and  in  the  second  column  of  Table  2.13.  A 
comparison  with  the  Snyder-Neufeld  results  of  Table 
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Tible  2.13.  Neutfon  Maximum  Abiorbad  Dom  Within 
Stindaid-Man  Slab  Phantoma^ 


Dom  ((rad!i/g)/(ncutron/cm^)|^ 


Source 

Energy 

(eV) 

Hearn  Source 

Isotropic  Flux 

Inflititc 

Slab 

Finite 

Slab‘' 

Infinite 

Slab 

Finite 

Slab‘' 

1.4208(+7)‘' 

7.9598(-9) 

7.4605(-9) 

1.46l6(-8) 

l.3638(-8) 

1.28S6(-f7) 

7.3683(-9) 

6,8683(-9) 

1.3436(-8) 

l.2S53(-8) 

1.1633(+7) 

7.2273(-9) 

6.771 2(-9) 

1.3077(-8) 

1.2274(-8) 

1.0526(+7) 

6.8462(-9) 

6.4060(-9) 

1.2342(-8) 

l.lS92(-8) 

1.0000(+7) 

6.7111(-9) 

6.2789(-9) 

1.2070(-8) 

1.1341(-8) 

9.5242(+6) 

6.574S(-9) 

6.144S(-9) 

1.1798(-8) 

1.1091(-8) 

8.6178(+6) 

6.3415(-9) 

5.91S3(-9) 

1.1308(-8) 

1.064 1(-8) 

7.7977(+6) 

6.3362(-9) 

5.9281(-9) 

1.1207(-8) 

1.0574(-8) 

7.0557(+6) 

6.081 1(-9) 

5.6837(-9) 

1.0736(-8) 

1.0136(-8) 

6.3843(+6) 

5.8372(-9) 

5.4736(-9) 

1.0319(-8) 

9.7712(-9) 

5.7767(+6) 

5.7003(-9) 

5.3445(-9) 

1.0002(-8) 

9.2060(-9) 

5.2270(+6) 

5.4689(-9) 

5.1380(-9) 

9.6108(-9) 

8.8640(-9) 

S.0000(+6) 

5.4174(-9) 

9.4514(-9) 

4.7296(+6) 

5.3667(-9) 

5.0331(-9) 

9.0074(-9) 

8.S870(-9) 

4.2795(+6) 

5.2731(-9) 

4.9584(-9) 

8.8140(-9) 

8.4324(-9) 

3.5038(+6) 

5.2205(-9) 

4.9916(-9) 

8.S412(-9) 

8.21 10(-9) 

3.1703(+6) 

4.7108(-9) 

4.4959(-9) 

7.7798(-9) 

7.4682(-9) 

2.5956(+6) 

4.1653(-9) 

3.9630(-9) 

6.8302(-9) 

6.5576(-9) 

2.0000(+6) 

3.8389(-9) 

6.1462(-9) 

1.9229(+6) 

3.7840(-9) 

3.6060(-9) 

6.0154(-9) 

5.7972(-9) 

1.5743(+6) 

3.4752(-9) 

3.3079(-9) 

5.4326(-9) 

5.2392(-9) 

1.0553(+6) 

3.1381(-9) 

2.9991(-9) 

4.6764(-9) 

4.5326(-9) 

5.2405(+5) 

2.1666(-9) 

2.0390(-9) 

2.9434(-9) 

2.8266(-9) 

5.0000(+5) 

2.1724(-9) 

2.0494(-9) 

2.9096(-9) 

2.7980(-9) 

4.7418(+5) 

2.1766(-9) 

2.0577(-9) 

2.8652(-9) 

2.7590(-9) 

4.2906(+5) 

2.1845(-9) 

2.6712(-9) 

1.9279(+5) 

1.3014(-9) 

1.1931(-9) 

1.5449(-9) 

1.4506(-9) 

9.8803(+4) 

9.2925(-10) 

8.2811(-10) 

1.0393(-9) 

9.5232(-10) 

4.667 1(+4) 

6.5317(-10) 

5.5685(-10) 

6.6372(-10) 

5.7892(-10) 

2.8308(+4) 

5.3582(-10) 

4.4115(-10) 

5.1846(-10) 

4.3284(-10) 

1.0414(-^4) 

4.5827(-10) 

3.5333(-10) 

4.2532(-10) 

3.2792(-10) 

4.9191(+3) 

4.6004(-10) 

3.5522(-10) 

4.2712(-10) 

3.2960(-10) 

1.0976(+3) 

5.2192(-10) 

4.0294(-10) 

4.3888(-10) 

3.4036(-10) 

3.1447(+2) 

5.3477(-10) 

4.1505(-10) 

4.4894(-10) 

3.4994(-10) 

9.0097(+l) 

S.4829(-10) 

4.2781(-10) 

4.5880(-10) 

3.5930(-10) 

3.3145(+1) 

5.5774(-10) 

4.3691(-10) 

4.6500(-10) 

3.6674(-10) 

9.4962(0) 

5.6353(-10) 

4.4345(-10) 

4.7006(-10) 

3.7620(-10) 

4.4857(0) 

5.6082(-10) 

4.4543(-10) 

4.7066(-10) 

3.7836(-10) 

2.1189(0) 

5.5112(-10) 

4.4173(-10) 

4.6484(-10) 

3.7528(-10) 

1.0009(0) 

5.2998(-10) 

4.2684(-10) 

4.4830(-10) 

3.6328(-10) 

6.0707(-l) 

5.0478(-10) 

4.0796(-10) 

4.2706(-10) 

3.4692(-10) 

4.7279(-l) 

4.8934(-10) 

3.9628(-10) 

4.1224(-10) 

2.1800(-1) 

3.5443(-10) 

2.9505(-l()) 

1.9709(-10) 

1.6567(-10) 

"From  ref.  43. 

^These  values  correspond  to  one  neutron  per  cm^  entering  the  slab. 

‘^The  calculations  for  the  finite  slab  include  a  buckling  term  to  compensate  for 
transverse  leakage;  the  transverse  dimensions  for  the  Tinite  slab  phantom  were 
taken  as  40  by  61  cm. 

‘'Read:  1.4208  X  lO’.etc. 
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Fig.  2,21,  Neutron  Maximum  Absoibed  Doaea  in  Standard-Man  Slab  Phantom  Due  to  a  Beam  Source.  (From  Ritts  et  aL,  ref.  43.) 


2.1 1  reveals  that  the  new  values  are  consistently  lower, 
as  would  be  expected  since  the  Snyder-Neufeld  calcula¬ 
tions  were  purposely  conservative. 

In  order  to  determine  the  effect  that  transverse 
leakage  would  have  on  the  calculated  doses,  Ritts  era/, 
also  performed  one-dimensional  calculations  for  a 
30-cm-thick  slab  with  a  geometric  buckling*  value  that 
was  equivalent  to  a  length  of  61  cm  and  a  width  of  40 
cm.  In  general  they  found  that  a  S  to  10%  reduction  in 
the  dose  could  be  attributed  to  transverse  leakage.  The 
results  for  the  finite  slab  are  included  in  Fig.  2.21 
(dashed  line)  and  in  the  third  column  of  Table  2.13. 
These  values  for  the  finite  slab  are  10  to  20%  lower 
than  the  doses  calculated  by  Snyder  and  Neufeld  for 
the  slab  with  no  transverse  leakage. 

Ritts  et  al.  also  made  calculations  with  and  without  a 
buckling  value  for  an  isotropic  neutron  flux  incident  on 
a  slab.  The  results  are  presented  in  Fig.  2.22  and  in  the 
last  two  columns  ofTable  2.13. 


*The  geometric  buckling,  is  the  lowest  eigenvalue  that 
results  from  solving  the  wave  equation  for  a  system  of  particular 
size  and  shape  subject  to  the  condition  that  the  flux  vanish  at 
the  extrapolated  boundary.  In  one-dimensional  calculations, 
leakage  in  the  infinite  directions  can  be  approximated  by  using 
DB'^  (£>  =  diffusion  coefficient)  as  a  fictitious  macroscopic 
absorption  cross  section. 


Gamma-Ray  Doses  in  Phantoms 

The  doses  delivered  to  various  parts  of  a  man-like 
phantom  by  0.027-  to  1.24-MeV  gamma  rays  were 
measured  by  Jones,’  ^  who  later  reported  his  results  in 
terms  of  calibration  factors”  for  dosimeters.  The 
calibration  factors  were  recommended  tt  use  in  the 
design  and  calibration  of  gamma-ray  exposure  meters  to 
permit  a  direct  reading  of  the  maximum  dose  (in  rads) 
delivered  to  the  critical  organ  for  the  particular  environ¬ 
mental  condition.  The  composition  of  the  phantom 
matched  soft  tissue  with  respect  to  density  and  effec¬ 
tive  atomic  number  and  included  a  full  set  of  bones  for 
the  head  and  trunk. 

For  each  of  1 1  photon  energies  used  in  the  experi¬ 
ment  (see  Table  2.14),  measurements  were  made  of  the 
free-field  dose^  (the  dose  with  the  phantom  removed) 
at  a  distance  from  the  source  corresponding  to  a 
particular  dose  position  in  the  phantom.  Measurements 
within  the  phantom  were  made  at  positions  corre¬ 
sponding  to  the  location  of  the  brain  (front),  eye  (lens), 
skin  (front  and  back  of  chest  and  forehead),  gut 
mucosa,  ovarie:.,  testes,  and  bone  marrow  (seven  repre¬ 
sentative  sites).  For  each  source  energy  the  measure- 


i'This  is  effectively  a  measurement  of  the  first-collision  dose 
as  described  in  Section  2.9.2. 
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Fig.  2,22.  Neutron  Maximum  Absorbed  Doses  in  Standard-Man  Slab  Phantom  Due  to  an  Isotropic  Source.  (From  Ritts  et  aL,  ref. 
43.) 


ments  were  made  with  the  phantom  being  irradiated 
from  in  front,  from  behind,  and  while  being  rotated 
about  the  vertical  axis.  Table  2.14  shows  the  location  at 
which  the  maximum  dose  was  received  within  the 
phantom  for  each  exposure  to  gamma  rays  of  different 
energies.  The  results  are  given  in  terms  of  the  ratio  of 
the  absorbed  dose  to  the  free-field  exposure  (rads/r). 

Recent  Monte  Carlo  calculations  by  Sidwell  et  al.^* 
for  monoenergetic  gamma  rays  incident  on  a  phantom 
assumed  to  have  the  characteristics  of  water  and  to  be 
shaped  like  an  elliptic  cylinder  yielded  dose  results  that 
are  reported  to  be  in  good  agreement  with  Jones’ 
measurements. 

Other  recent  calculations  of  the  dose  delivered  to  a 
phantom  by  monoenergetic  gamma  rays  were  per¬ 
formed  by  Claiborne  and  Trubey.®*  Using  the  discrete 
ordinates  method,  with  some  checks  by  the  Monte 
Carlo  method,  they  calculated  the  dose  distributed  in  a 
30-cm-thick  infinite-slab  phantom  having  the  standard 
man  composition  shown  in  Table  2.8.  Their  maximum 
dose  rates,  which  occurred  at  nearly  the  surface  or 
within  the  first  2  cm  of  the  surface,  are  compared  in 
Fig.  2.23  with  Jones’  experimental  results  (converted  to 
the  same  units).  First-collision  dose  rates  for  the 
standard-man  composition  based  on  the  use  of  f1ux-to- 
kerma-rate  conversion  factors  are  also  plotted  in  Fig. 
2.23.  These  values  show  that  the  customary  practice  of 
using  kerma-rate  factors  underestimates  the  maximum 
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Fig.  2.23.  Maximum  Gamma-Ray  Dom  Rate  in  Slab  Phantom 
of  Standaid-Man  Compoiition.  Comparison  with  kemia  calcula¬ 
tions  and  with  phantom  measurements.  (From  Claiborne  and 
Trubey,  ref.  55.) 
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Table  2.14.  Ganiina-Riy  Maximum  Abiurbad  Dorm  al  Vaiiuui 
Locationi  In  a  Man>Llke  Phantom'* 


Source 

Gonima-P.uy 

Fnergy 

(MeV) 

Location  in  Phantom 

Ratio  of 
Maximum  Dose 
to  Exposure 
(rads/r) 

‘®Co 

1.24 

Skin  on  back^ 

1.06 

‘”Cs 

0.660 

Skin  on  chest 

1.11 

•”ln 

0.400 

Skin  on  forehead  and  testes 

1.19 

Bremsstrahlung 
Nanow  band 

0.240 

Skin  on  forehead 

1.30 

0.165 

Gut  mucosa 

1.47 

0.125 

Te.stes 

1.57 

Filtered  by  U 

o.too 

Skin  on  forehead 

1.65 

Filtered  by  Pb 

0.075 

Gut  mucosa 

2.03 

Filtered  by  Ta 

0.060 

Skin  on  chest 

1.75 

Filtered  by  Nd 

0.038 

Skin  on  chest 

1.34 

Filtered  by  Sn 

0.027 

Skin  on  chest 

1.28 

‘*From  ref,  52. 

^Maximum  occurred  when  phantom  was  irradiated  from  behind;  all  other 
maxima  occurred  when  phantom  was  irradiated  from  in  front. 


dose  delivered  to  a  position  within  a  body,  especially 
for  the  lower  gamma-ray  energies.  This  is  because  the 
energy  deposition  due  to  collisions  subsequent  to  the 
first  collisions  is  not  accounted  for  in  the  kerma-rate 
factois.* 

The  excellent  agreement  between  the  Claiborne- 
Trubey  calculations  and  the  Jones  experiment  for 
energies  greater  than  0.07  MeV  seems  surprising  in  view 
of  the  large  differences  in  the  conflgurations.  However, 
since  most  of  the  energy  deposition  apparently  comes 
from  the  first  collision,  the  effect  due  to  the  geometry 
difference  will  be  small. 


2.9.6.  MULTICOLLISION  DOSE 

The  term  multicollision  dose,  which  until  recently  has 
been  used  solely  for  neutron  exposure,  is  not  recog¬ 
nized  by  the  ICRU.  The  multicollision  dose  concept 
was  developed  by  health  physicists  in  an  effort  to  relate 
the  neutron  flux  or  fluence  incident  upon  a  configura¬ 
tion  of  tissue  representing  the  human  body  to  the 
absorbed  dose  or  to  the  dose  equivalent.  The  multi¬ 
collision  dose  could,  in  turn,  be  related  to  a  maximum 


*The  dose  rate  in  a  slab-tissue  phantom  due  to  high-energy 
photons  has  been  caicuiated  by  Alsmiller  and  Moran, but 
except  for  the  values  for  10-MeV  gamma  rays  their  results  are 
for  photon  energies  not  encountered  in  weapons  radiation  or 
reactor  shielding. 


permissible  dose  (MPD)  established  by  shield  design 
criteria.  The  term  has  tacitly  been  defined  by  common 
usage  in  shield  design  as  the  sum  of  the  maximum  dose 
equivalents  (or  maximum  absorbed  doses)  delivered  by 
each  neutron  energy  group  to  a  slab  of  tissue  of  a 
particular  composition.*  This  includes  the  dose  due  to 
secondary  gamma  rays  produced  by  the  neutrons  in  the 
slab.  This  dose  is  calculated  by 

Multicollision  dose  =  2]  'I'g <  (2.142) 

g 

where  yjjg  is  the  group  fluence  as  defined  in  Section  2.1 
and  D^.^^{Eg)  is  the  maximum  dose  equivalent  (or 
maximum  absorbed  dose)  in  the  slab  for  an  incident 
neutron  of  the  energy  corresponding  to  the  ^h  neutron 
group. 

Gamma-ray  multicollision  doses  are  similarly  calcu¬ 
lated  with  Eq.  2.142,  with,  of  course,  \pg  and 
being  the  appropriate  values  for  gamma  rays.  Because 
the  gamma-ray  maximum  absorbed  doses  reported  in 
Section  2.9.5  have  only  recently  become  available,  most 
multicollision  doses  for  gamma  rays  have  in  the  past 
been  determined  by  using  the  product  of  the  gamma- 
ray  energy,  the  mass  energy  deposition  coefficient  for 
tissue  or  air,  and  a  factor  to  convert  to  rads  for  the 


*See  slab  tissue  phantom  described  in  Section  2.9.5  for 
Snyder-Neufeld  calculations.** 
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The  total  dose  for  a  mixed  radiation  Held  incident  on 
a  tissue  stub  is  determined  by  adding  the  neutron  and 
gamma-ray  multicollision  doses  obtained  separately  for 
the  slab. 

Although  the  concept  of  multicollision  dose,  or  dose 
rate,  is  not  recognized  by  the  ICRU,  it  seems  clear  that, 
to  fully  implement  the  trend  to  preciseness  in  the  fields 
of  radiobiology  and  radiation  protection,  the  dose 


calculated  in  shield  design  should  be  given  recognition 
by  some  official  body.  Perhaps  the  term  ntt4lllcollisi()n 
dose  can  be  retained  for  the  neutron-induced  contri¬ 
bution  and  the  term  total  multicolllsion  dose  can  be 
used  when  the  incident  gamma-ray  contribution  is  also 
included.  Or  perhaps  an  entirely  new  term  could  be 
introduced,  such  as  the  term  man  equivalent  dose 
(MED)  proposed  by  Claiborne  and  Trubey  .*  * 


2.10.  Radiation  Exposure  Limits* 


Any  materiul  will  undergo  some  alteration  on  ex¬ 
posure  to  directly  or  indirectly  ionizing  radiation,  hut 
organic  materials  are  most  susceptible  to  damage, 
particularly  living  organisms.  The  earth  and  its  inhabit¬ 
ants  are  continuously  exposed  to  radiation  from 
cosmic-ray  bombardment  and  naturaily  occurring  radio¬ 
isotopes.  In  addition  to  these  natural  sources,  the 
human  body  may  be  exposed  to  X  rays  used  for 
medical  diagnostics  and  treatments,  luminous  wrist 
watch  dials,  television,  etc.  To  date  the  effects  of  these 
exposures,  if  any,  are  essentially  unknown. 

The  International  Commission  on  Radiologicai  Pro¬ 
tection  (ICRP)  has  taken  the  position  that  radiation  in 
any  amount  absorbed  by  the  gonads  can  cause  genetic 
damage  and  increase  the  natural  mutation  rate.  Mostly 
as  a  result  of  the  fear  that  small  chronic  dosages  will 
significantly  increase  the  mutation  rate  in  humans  or 
induce  unknown  deleterious  long-range  effects  (such  as 
leukemia),  the  iCRP  has  recommended  that  workers  in 
the  nuclear  industry  be  limited  to  a  maximum  dose  of  5 
rems/year,  which  amounts  to  2.5  mrems/hr  based  on  a 
40-hr  work  week  for  50  weeks/year.  This  so-called 
maximum  permissible  dose  (MPD)  is  about  25  times  the 
dose  that  the  average  American  will  receive  each  year 
from  environmental  radiation.  The  recommended  MPD 
appears  to  represent  an  acceptable  risk  since  con¬ 
servative  estimates  indicate  that  5  rems/year  may 
increase  the  natural  gene  mutation  rate  by  only  10  to 
20%,  which  when  related  to  the  probable  increase  in  the 
number  of  potentially  defective  children  becomes  a 
statistical  exercise.  It  is  interesting  to  note  that  no 
measurable  defects  attributable  to  radiation  were  found 
in  studies  of  the  children  conceived  after  the  parents 
were  exposed  to  intense  radiation  fields  from  the 
atomic  weapons  exploded  over  Hiroshima  and  Nagasaki 
and  to  the  Marshall  Island  natives  exposed  to  fallout 
from  the  Bikini  test  shot. 

Much  higher  doses  are  required  to  produce  significant 
somatic  damage  or  direct  deleterious  effects  in  the 


♦This  section  is  iargely  quoted  from  Sections  2-8.4  and  2-8.5 
of  Volume  I  of  the  Office  of  Civil  Defense's  manual  for  shelter 
design  and  analysis  (see  ref.  57). 


body.  And  since  the  body  cun  repair  itself,  dose  rule  is  u 
factor,  a  larger  total  radiation  dose  being  required  to 
produce  a  given  degree  of  injury  if  the  dose  is  spread 
over  a  long  period  of  time.  But  the  rate,  or  fraction¬ 
ation,  effect  is  of  no  consequence  for  acute  doses  that 
could  be  received  in  a  nuclear  war  or  from  an  accident 
in  the  nuclear  industry.  Military  considerations  during 
war  require  emphasis  on  minimizing  acute  exposure  and 
the  attendant  probability  of  an  incapacitation  that 
would  prevent  a  necessary  mission.  Consequently  the 
dose  delivered  to  mucosal  layers  of  the  gut  becomes  the 
most  important  consideration. 

Althougli  data  on  the  biological  effects  that  might  be 
expected  for  human  beings  subjected  to  acute  doses  of 
varying  amounts  are  rather  limited,  a  sufficient  amount 
of  knowledge  is  available  for  drawing  useful  conclusions 
for  planning  purposes.  Such  data  have  been  derived 
from  the  study  of  Japanese  casualties,  a  few  nuclear 
accidents,  and  large  animal  experiments.  The  following 
paragraphs  give  a  general  description  of  the  effects  that 
might  be  expected  as  the  result  of  whole-body  exposure 
to  an  acute  dose  received  by  healthy  adults  over  a 
period  of  up  to  two  to  four  days.  For  a  more  complete 
and  detailed  discussion  see  ref.  58. 

Doses  up  to  too  rems.  -  Exposures  resulting  in  doses 
up  to  100  rems  would  give  no  apparent  immediate 
effects.  In  some  cases  nausea  or  vomiting  may  occur, 
but  no  other  symptom  of  radiation  sickness  would  be 
apparent.  No  clinical  surveillance  or  therapy  other  than 
reassurance  would  be  required. 

Doses  from  100  to  200  rems.  -  Doses  between  100 
and  200  rems  would  cause  nausea  and  vomiting  within 
about  3  hr  in  from  5  to  50%  of  those  exposed.  A 
moderate  decrease  in  the  number  of  white  blood  cells 
would  result,  making  clinical  and  hematologic  surveil¬ 
lance  desirable.  Therapy  would  consist  mainly  of 
reassurance.  Prognosis  would  be  excellent  with  a  con¬ 
valescent  period  of  several  weeks.  There  should  be  no 
deaths. 

Doses  from  200  to  600  rems.  -  Doses  in  the  range 
200  to  600  rems  would  result  in  a  100%  incidence  of 
vomiting,  which  would  begin  within  about  2  hr  for 
doses  greater  than  300  rems.  Severe  loss  of  white  blood 
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cells  (Icukupeniu)  would  occur.  Purple  or  livid  spots 
(purpuru)  would  appear  on  the  skin.  There  would  be 
some  hemorrhaging  and  an  incidence  of  infection.  Loss 
of  hair  (epilation)  would  result  for  exposures  above  .100 
reins.  Therapy  would  include  blood  transfusions  and 
the  administration  of  antibiotics.  A  convalescent  period 
of  from  I  to  12  months  would  he  required.  Deaths 
from  hemorrhage  or  infection  within  two  months 
would  probably  begin  at  about  2S0  rcms  and,  for  the 
dose  range  of  200  to  600  rems.  would  range  up  to 
about  KO';; .  A  value  of  around  .100  rems  delivered  to  the 
mid-line  of  the  body  is  usually  ta'<en  as  LD;,gy  ,(,  (lethal 
dose  for  .*50%  within  .10  days)  for  man. 

Doses  from  600  to  1000  rems.  Doses  approaching 
1000  rems  would  result  in  all  the  symptoms  of 
radiation  sickness  described  in  the  preceding  paragraph 
but  at  a  more  severe  intensity.  Therapy  would  be  only 
promising  and  would  involve  not  only  blood  transfusion 
and  antibiotics  but  also  the  consideration  of  bone 
marrow  transplants.  The  incidence  of  deaths  would 
range  between  80  and  100%. 

Doses  Greater  Than  1000  rems.  -  Doses  greater  than 
1000  rems  would  almost  certainly  be  lethal.  Symptoms 
of  radiation  sickness  would  be  very  severe  and  would 
begin  within  .10  min.  Those  symptoms  described  in  the 
preceding  paragraphs  would  be  augmented  by  diarrhea, 
fever,  and  disturbance  of  electrolyte  balance  as  the 
result  of  changes  in  the  gastrointestinal  tract.  Doses 
over  5000  rems  would  affect  the  central  nervous  system 
with  attendant  symptoms  of  convulsion,  tremors,  loss 
of  muscular  coordination,  and  lethargy.  Prognosis 


would  be  hopeless  and  therapy  would  consist  mainly  of 
the  administration  of  sedatives.  For  doses  of  1000  to 
5000  rems,  death  in  alnxist  100%  of  the  cases  would 
result  within  2  weeks  from  circulatory  collapse.  For 
doses  above  5000  rems  death  in  100%  of  the  cases 
would  occur  within  48  hr  as  a  respiratory  failure  or 
edenu  of  the  brain,  or  both.  Doses  of  10,000  rems  and 
up  directly  affect  the  nerve  cells  to  the  point  that 
incapacitation  could  occur  almost  immediately. 

To  summarize,  if  a  person  is  exposed  to  a  dose  of  less 
than  about  200  rems,  he  should  not  become  incapaci¬ 
tated  nor  should  his  ability  to  work  be  seriously 
affected.  Persons  exposed  to  doses  exceeding  about  200 
rems  will  suffer  increasing  radiation  sickness  with 
increasing  dosage,  and  the  probability  of  death  is 
extremely  higlt  if  the  dose  exceeds  600  rems.  For  doses 
greater  than  10,000  rems  immediate  incapacitation  can 
be  expected. 

The  MPD  for  military  purposes  during  war  will 
depend  on  the  assigned  mission  an !  the  acceptable 
casualty  rate.  Obviously  the  low  MPD  set  for  peacetime 
occupational  work  cannot  be  used  as  the  criterion  in 
military  situations.  Based  on  the  available  information, 
an  MPD  =  100  rems  (absorbed  over  a  few  days)  seems 
reasonable  since  the  effect  of  such  a  dose  would  be 
hardly  noticeable  and  would  not  interfere  with  military 
operations.  Since  radiation  shelter  design  will  require  an 
estimate  of  the  free-field  dose,  the  maximum  probable 
exposure  in  a  specified  time  (which  includes  multiple 
strikes)  should  be  used  as  the  design  criterion  for  each 
shelter  in  a  given  area. 
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Appendix  2A.  f  ieutron  Fluence-to-Kerma  Conversion  Factors 
for  Standard  Man  and  Some  Specific  Body  Components 


This  appendix  presents  neutron  fluence-to-kerma  con¬ 
version  factors  calculated  by  Ritts  el  al.*^  for  a 
standard-man  model  and  for  specific  body  components 
(see  Section  2.9.2).  The  neutron  energies  range  from 
approximately  0.022  eV  to  approximately  IS  MeV,and 
the  compositions  for  the  models  are  those  shown  in 
Table  2.8.  Table  2A.1  gives  the  conversion  factors  for 
the  standard  man,  together  with  those  for  the  elements 
making  the  greatest  contributions  in  the  model.  Con¬ 
version  factors  for  specific  body  components  are 
presented  in  Table  2A.2. 

Table  2A.3  gives  the  conversion  factors  on  an  atom 
density  basis  for  the  most  common  elements  contained 
in  man.  These  values  can  be  used  to  obtain  the  kerma 
factors  for  any  arbitrary  combination  of  these  elements 
by  multiplying  the  individual  value  of  each  component 


by  the  number  of  atoms  of  each  element  per  gram  of 
mixture.  For  example,  the  kerma  factor  for  lIjO  ( 1 1 .9 
wt  %  H  and  88.81  wt  %  0)  in  a  field  of  15-MeV 
neutrons  is 


0.1 19 
1.008 

0.881 


X  (0.6023  X  10"'’ )  (7.75  X  lO'") 


16.00 


X  (0.6023  X  10"“)  (6.67  X  10  '’) 


X  10’ 


=  7.41  X  10’’ (ergs/g  of  n2  0)/(neutron/cm" )  .  (Al) 


where  the  numbers  1 .008  and  16  are  the  atomic  weights 
of  hydrogen  and  oxygen  respectively  and  0.6023  X 
I0"“  is  Avogadro’s  number. 


Table  2A.I.  Neutron  Fluence-to-Kerma  Factors  lor  Standard  Man"  * 


£nergy  Fluence-to-Kerma  Factor  (ergs/g)/(neutron/cm") 

(MeV)  Total  Hydrogen  Oxygen  Carbon  Nitrogen  Phosphorus  Sodium  Chlorine 


15.14 

7.0791-7)" 

4.633(-7) 

l.S3(  7) 

8.08(-8) 

9.35(-9) 

14.84 

6.948(-7) 

4.629(-7) 

1.46(-7) 

7.53(-8) 

9.l3(-9) 

14.54 

6.811(-7) 

4.623(-7) 

1.40(  7) 

6.89(-8) 

8.87(-9) 

14.24 

6.652(-7) 

4.617(-7) 

1.29(  7) 

6.47(-8) 

8.6  7( -9) 

13.93 

6.501(-7) 

4.61 1( -7) 

l.l8f-7) 

6.16(-8) 

8.41(  9) 

13.63 

6.276(-7) 

4.604(-7) 

9.90(-8) 

5.93(-8) 

7.99(  9) 

13.33 

6.1 19(-/) 

4.598(-7) 

8.63(-8) 

5.72(-8) 

7.74(-9) 

13.03 

6.007(-7) 

4.591(-7) 

7.87(  8) 

5.44(  8) 

7.63(-9) 

12.73 

5.966(-7) 

4.584(-7) 

8.03(-8) 

4.96(  8) 

7.44(  9) 

12.43 

5.949(-7) 

4.577(-7) 

8.30(-8) 

4.61(-8) 

7.29(-9) 

12.13 

5.936(-7) 

4.568(-7) 

8.45(-8) 

4.45(-8) 

7.04(-9) 

11.83 

5.964(-7) 

4.559(-7) 

9.28(  8) 

4.01(-8) 

6.77(  9) 

11.52 

5.877(-7) 

4.550(-7) 

8.90(-8) 

3.64(-8) 

6.47(  9) 

11.22 

5.818(-7) 

4.540(  -7) 

8.48(  8) 

3.59(-8) 

6.23(  9) 

10.92 

5.691(-7) 

4.531(-7) 

7.49(-8) 

3.43(-8) 

5.95(  9) 

10.62 

5.554(-7) 

4.522(-7) 

6.75(-8) 

2.91(-8) 

5.85(-9) 

10.32 

5.483(-7) 

4.512(-7) 

6.81{-8) 

2.:6(-8) 

5.73(  9) 

10.02 

5.421(-7) 

4.50U-7) 

6.46(-8) 

2.12(  8) 

5.52(  9) 

9.72 

5.353(-7) 

4.478(-7) 

5.98(-8) 

2.18(-8) 

5.28(  9) 

9.60 

5.313(-7) 

4.468(-7) 

5.70(  -8) 

2.16(-8) 

5.21(-9) 

9.45 

5.271(-7) 

4.456(-7) 

5.40(-8) 

2.19(-8) 

4.97(~9) 

9.30 

5.222(-7) 

4.443(-7) 

5.00(-8) 

2.23(-8) 

4.91(-9) 

9.15 

5.173(-7) 

4.430(-7) 

4.74(-8) 

2.14(-8) 

4.85(  9) 

9.00 

5.123(-7) 

4.417(-7) 

4.78(-8) 

1.75(-8) 

4.70(-9) 

8.85 

5.078(-7) 

4.410(-7) 

4.54(-8) 

1.61(-8) 

4.64(-9) 

8.70 

5.048(-7) 

4.402(-7) 

4.37(-8) 

1.56(-8) 

4.60(-9) 

8.62 

5.007(-7) 

4.399(-7) 

4.01(-8) 

1.56(-8) 

4.53(-9) 

8.47 

4.988(-7) 

4.391{-7) 

3.90(-8) 

1.57(-8) 

4.40(-9) 

8.32 

5.036(-7) 

4.383(-7) 

4.38(-8) 

1.66(-8) 

4.28(-9) 

8.17 

5.062(-7) 

4.376(-7) 

4.24(-8) 

2.15(-8) 

4.30(-9) 
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Table  2A.I.  (continued) 


Eneijy  Fhience-to-Kenna  Factor  (ergs/g)/( neutron/ cm^) 

(MeV)  7otal  Hydrogen  Oxygen  Carbon  Nitrogen  Phosphorus  Sodium  Chiorine 


8.02 

4.977(-7) 

4.368(-7) 

2.94(-8) 

2.63(-8) 

4.74(-9) 

7.87 

5.038(  7) 

4.345(-7) 

3.81(  8) 

2.62(-8) 

4.38(-9) 

7.72 

5.060(  7) 

4.319(-7) 

3.89(-8) 

3.09(-8) 

3.88(-9) 

7.57 

4.935(-7) 

4.292(-7) 

3.65(  8) 

2.34(-8) 

3.84(-9) 

7.42 

4.970(  7) 

4.266(-7) 

4.11{-8) 

2.49(-8) 

3.98(-9) 

7.27 

4.995(  7) 

4.239(-7) 

5.32(-8) 

1.80(-8) 

3.87(-9) 

7.12 

4.773(  7) 

4.21 1(-7) 

4.13(-8) 

1.08(-8) 

3.52(-9) 

6.97 

4.650(  7) 

4.187(-7) 

3.13(-8) 

1.10(-8) 

3.49(-9) 

6.82 

4.745(  7) 

4.17j(-7) 

4.14(-8) 

1.21(-8) 

3.27(-9) 

6.67 

4.707(  7) 

1.159(-7) 

3.94(-8) 

l.i7(-8) 

3.27(-9) 

6.52 

4.605(  7) 

4.145(-7) 

2.97(-8) 

1.25(-8) 

3.38(-9) 

6.37 

4.731(  7) 

4.130(-7) 

3.89(-8) 

1.71(-8) 

3.57(-9) 

6.21 

4.61 1(-7) 

4.115(-7) 

2.46(-8) 

2.09(-8) 

3.65(-9) 

6.06 

4.572(  7) 

4.100(-7) 

2.77(-8) 

1.53(-8) 

3.77(-9) 

5.91 

4.633(-7) 

4.077(-7) 

3.73(-8) 

1.42(-8) 

3.67(-9) 

5.76 

4.457(-7) 

4.048(-7) 

2.32(-8) 

1.38(-8) 

3.51(-9) 

5.61 

4.502(  7) 

4.018(-7) 

3.08(-8) 

1.37{-8) 

3.55{-9) 

5.46 

4.312(-7) 

3.987(-7) 

1.47{-8) 

1.42(-8) 

3.25(-9) 

5.31 

4.294(-7) 

3.956(-7) 

1.70(-8) 

1.30(-8) 

3.39(-9) 

5.16 

4.449(-7) 

3.925(-7) 

3.46(-8) 

1.33(-8) 

4.17(-9) 

5.01 

4.344(  7) 

3.892(-7) 

2.76(  8) 

1.32(-8) 

4.01(-9) 

4.69 

4.180(-7) 

3.810(-7) 

1.79(-8) 

1.50(-8) 

3.73(-9) 

4.50 

4.116(-7) 

3.760(-7) 

1.66(-8) 

1.42(-8) 

4.34(-9) 

4.31 

4.181(-7) 

3.709(-7) 

2.91(-8) 

1.33(-8) 

4.48(-9) 

4.16 

4.245(-7) 

3.666(-7) 

3.37(-8) 

1.89(-8) 

4.98(-9) 

4.01 

4.102(-7) 

3.623(-7) 

2.24(-8) 

2.05(-8) 

4.78(-9) 

3.82 

4.121(  7) 

3..562(-7) 

2.91(-8) 

2.26(-8) 

3.88(-9) 

3.67 

4.099(-7) 

3.S11(  7) 

3.05(-8) 

2,42(-8) 

3.67(-9) 

3.52 

o 

o 

1 

3.459(-7) 

3.09(-8) 

2.45(-8) 

4.48(-9) 

3.33 

3.981(  7) 

3.392(-7) 

3.46(-8) 

2.07(-8) 

3.32(-9) 

3.18 

3.337(-7) 

2.06(-8) 

1.54(-8) 

2.83(-9) 

3.03 

3.555(  7) 

3.280(-7) 

1.53(-8) 

9.44(-9) 

2.47(-9) 

2.84 

3.499(  7) 

3.203(-7) 

1.14(-8) 

1.60(-8) 

2.01(-9) 

2.69 

3.396(-7) 

3.138(-7) 

1.10(  8) 

1.26(-8) 

1.89{-9) 

2.50 

3.279(  7) 

3.054(-7) 

1.01(-8) 

1.09(-8) 

1.34(-9) 

2.40 

3.150(  7) 

3.000(-7) 

3.40(-9) 

1.03(-8) 

1.07(-9) 

2.31 

3.136(  7) 

2.950(-7) 

7.36(  9) 

9.85(-9) 

1.21(-9) 

2.21 

3.1 19(  7) 

2.898(-7) 

1.08(-8) 

9.51(-9) 

1.61(-9) 

2.10 

3.057(-7) 

2.835(-7) 

1.14(-8) 

9.68{-9) 

9.62(-10) 

2.00 

2.995(-7) 

2.780(-7) 

1.13(-8) 

9.16(  9) 

9.20(-10) 

1.91 

2.952(-7) 

2.721(-7) 

1.31(-8) 

8.82(-9) 

9.82(-10) 

1.80 

2.883(-7) 

2.648(-7) 

1.31(-8) 

8.56(-9) 

1.72(-9) 

1.70 

2.800(  7) 

2.585(-7) 

1.21{-8) 

8.36(-9) 

8.76(-10) 

1.59 

2.717(-7) 

2.507(-7) 

1.18(-8) 

8.12(-9) 

9.84(-10) 

1.50 

2.644(-7) 

2.440(-7) 

1.15(-8) 

7.93(-9) 

9.06(-10) 

1.40 

2.577(-7) 

2.370(-7) 

1.13(-8) 

7.73(-9) 

1.66(-9) 

1.31 

2.602(-7) 

2.297(  7) 

2.23(-8) 

7.52(-9) 

6.44(-10) 

1.20 

2.426(-7) 

2.209(-7) 

1.39(-8) 

7.25(-9) 

5.11(-10) 

1.10 

2.371(-7) 

2.121(-7) 

1.63(-8) 

7.00(-9) 

7.45(-10) 

1.00 

2.432(  7) 

2.037(-7) 

3.21(-8) 

6.70(-9) 

6.33(-10) 

9.0(  1) 

2.122(  7) 

1.939(-7) 

1.15(-8) 

6.35(-9) 

3.39(-10) 

8.0(  1) 

1.980(-7) 

1.834(  7) 

7.95(  9) 

5.97(-9) 

3.90(-10) 

7.0(  1) 

1.857(  7) 

o 

1 

6.6U-9) 

5.56(  9) 

4.97(-10) 

6.0(-l) 

1.719(-7) 

1.609(-7) 

5.61(  9) 

5,07(-9) 

2.80(-10) 

5.0(-l) 

1.592(-7) 

1.48.3(-7) 

5.74(-9) 

4.52(-9) 

5.32(-10) 

4.5(-l) 

1.639(-7) 

1.403(-7) 

1.90(-8) 

4.20(-9) 

2.90(-10) 
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Table  2A.1.  (continued) 


Energy 

(MeV) 

Fhjence-to-Kerma  Factor  (etgs/g)/(neutron/cm*) 

Total 

Hydrogen 

Oxygen 

Carbon 

Nitrogen 

Phosphorus 

Sodium 

Chlorine 

4.0(-l) 

1.486(-7) 

1.318(-7) 

1.18(-8) 

3.89(  9) 

2.96(  10^ 

3.5(-l) 

l.339(-7) 

1.236(-7) 

6.52(-9) 

3.5U-9) 

2.46(  10) 

3.0(-l) 

1.223(-7) 

1.142(  7) 

4.73(  9) 

3.I2(  9) 

2.13(-10) 

2.5(-l) 

1.104(-7) 

1.037(-7) 

3.74(-9) 

2.70(-9) 

1.93(  10) 

2.0(-l) 

9.649(-8) 

9.1l6(-8) 

2.89V-9) 

2.23(  9) 

1.82(-10) 

1.5(-1) 

8.165(-8) 

7.758(-8) 

2.16(-9) 

1.75(-9) 

1.47(-10) 

l.O(-l) 

6.335(-8) 

6.056(-8) 

1.45(-9) 

1.22(  9) 

1.08(-10) 

7.0(-2) 

4.973(-8) 

4.774(-8) 

1.02(-9) 

8.76(-10) 

9.36(-ll) 

5.0(-2) 

3.830(-8) 

3.705(  8) 

7.30(  9) 

6.37(  10) 

7.62(  11) 

3.0(-2) 

2.560(-8) 

2.472(-8) 

4.43(-9) 

3.92(-10) 

4.97(  11) 

2.0(-2) 

1.807{-8) 

1.747(-8) 

2.98(-9) 

2.66(  10) 

3.73(-ll) 

1.5(-2) 

l.361(-8) 

1.316(-8) 

2.18(-9) 

1.95(-10) 

3.14(  ID 

1.0(-2) 

9.550(-9) 

9.287(-9) 

1.47(-9) 

1.34(-10) 

2.29(-ll) 

7.0(-3) 

6.776(-9) 

6.558(-9) 

1.03(-9) 

9.24(  11) 

1.87(-11) 

5.0(-3) 

4.880(-9) 

4.832(-9) 

7.53{-9) 

6.74(-ll) 

l.S6(  11) 

4.03(-3) 

3.951(-9) 

3.820(-9) 

5.94(-ll) 

5.32(-ll) 

1.44(-11) 

2.88(  12) 

3.00(-3) 

2.964(-9) 

2.855(-9) 

4.43{-ll) 

3.96(-ll) 

1.32(  11) 

3.02(  12) 

2.01(-3) 

2.002(-9) 

1.928(-9) 

2.98(-ll) 

2.67(-ll) 

1.27(-11) 

3.12(  12) 

1.50(-3) 

l.495(-9) 

1.436(-9) 

2.22(-ll) 

1.99(-11) 

1.30(  11) 

3.20(-12) 

1.03(-3) 

1.030(-9) 

9.832(-10) 

1.52(-ll) 

l.36(-ll) 

1.40(-11) 

3.25(-12) 

7.68(-4) 

7.774(-10) 

7.368(-10) 

l,14(-ll) 

1.02(-11) 

1.52(-11) 

3.27(-12) 

2.72(  13) 

S.13(-4) 

5.299(-10) 

4.925(-10) 

7.60(-12) 

6.80(-12) 

1.77(  11) 

3.30(  12) 

l.84(-12) 

4.02(-4) 

4.244(-10) 

3.867(-10) 

S.96(-12) 

5.i4(-12) 

1.95(-ll) 

3.31(  12) 

3.28(  12) 

3.I0(-4) 

3.347(-l0) 

2.986(-10) 

4.60(-12) 

4.12(-12) 

2.19(-U) 

3.32(-12) 

).97(-12) 

2.01(-4) 

2.299(-10) 

I.939(-10) 

2.98(-I2) 

2.67(-12) 

2.67(-ll) 

3.33(-12) 

5.91(  14) 

1.42(-4) 

1.764(-10) 

1.372(-10) 

2.n(-12) 

1.88(-12) 

3.15(-11) 

3.33(-12) 

1.12(  13) 

1.006(-4) 

1.417(-10) 

9.7S3(-11) 

1.49(-12) 

1.33(-12) 

3.73(-ll) 

3.34{-12) 

3.87(  13) 

7.1(-5) 

I.192(-10) 

6.920(-ll) 

1.05(-12) 

9.42(-13) 

4.42(-lI) 

.3.34(  12) 

I 

o 

9.84(-14) 

S.03(-5) 

1.071(-10) 

4.94(-ll) 

7.46(-13) 

6.67(-13) 

5.237(-ll) 

3.34(  12) 

4.11(-13) 

1.10(-13) 

4.1I(-S) 

l.036(-10) 

4.07(-ll) 

6.09(-13) 

5.45(-13) 

5.786(-ll) 

3.34(-12) 

4.54(-13) 

1.31(-13) 

3.43(-5) 

l.025(-10) 

3.43(-U) 

5.09(-l3) 

4.56t-13) 

6.323(-ll) 

3.34(-12) 

4.99(  13) 

1.57(-13) 

2.74(-5) 

!.033(-10) 

2.78(-ll) 

4.07(-13) 

3.64(  13) 

7.065(-ll) 

3.35(-12) 

5.58(-13) 

1.95(  1.3) 

2.05(-5) 

1.077(-10) 

2.I3(-11) 

3.05(-l3) 

2.73(-l3) 

8.1SS(  11) 

3.35(-12) 

6.44(  13) 

2.59(  13) 

1.60(-5) 

1.143(-10) 

1.7U-11) 

2.38(-l3) 

2.13(-13) 

9.226(-ll) 

3.35(-12) 

7.34(-13) 

3.14(-13) 

1.26(-5) 

l.230(-10) 

1.40(-ll) 

I.86(-I3) 

1.67(-I3) 

1.041(-10) 

3.35(  12) 

8.34(  13) 

3.79(  13) 

8.0(-6) 

l.454(-10) 

l.Ol(-ll) 

1.I9(-13) 

1.06(-13) 

1.302(-10) 

3.35(  12) 

1.04(-12) 

5.08(-13) 

4.0(-6) 

1.965(-10) 

7.20(-12) 

5.94(-14) 

5.32(  14) 

1.836(-10) 

3.35(-12) 

1.46(-12) 

7.2»(-13) 

2.0(-6) 

2.723(-10) 

6.63(-12) 

2.97(-14) 

2.66(-14) 

2.591(-10) 

3.35(-12) 

2.09(-12) 

l,a3(  12) 

1.0(-6) 

3.811(-10) 

7.59(-12) 

1.49(-14) 

1.33{-14) 

3.656(-10) 

3.35(  12) 

2.97(-12) 

1.45(-12) 

5.0(-7) 

5.354(-10) 

9.87(--12) 

5.159(-10) 

3.35(-12) 

4.15(-12) 

2.07{-12) 

2.5(-7) 

7.529(-10) 

1.35(-11) 

7.280(-10) 

3.35(  12) 

5.89(-12) 

2.12(  12) 

1.7(-7) 

I.I44(-9) 

2.03(-ll) 

1.109(-9) 

3.35(  12) 

9.00(  12) 

2.12(-12) 

7.2(-8) 

1.398(-9) 

2.48(-ll) 

1.356(-9) 

3.35(  12) 

l.lO(-ll) 

2.12(  12) 

4.5(-8) 

1.766(-9) 

3.14(-11) 

1.715(-9) 

3.35(-12) 

1.39(  11) 

2.12(  12) 

2.23(-8) 

2.489(-9) 

4.18(-11) 

2.424(-9) 

3.35(-12) 

1.83(  11) 

2.12(-12) 

“From  ref.  43. 

^he  standard-man  model  also  includes  the  elements  calcium,  sulfur,  potassium,  and  magnesium  but  their  contributions  to  the 
kerma  are  negligible. 

‘'Read:  7.079  X  lO*’,  etc. 
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Table  2A.2.  Neutron  Ruence-lo-Keima  Factors  forSpecinc  Body  Components" 


|.;nergy  Fluence-to-Kerma  Factors  (ergs/g)/(neutron/cm^) 


(MeV) 

tung 

Muscle 

Bone 

Brain 

Red  Marrow 

15.00 

7.04(  7)'’ 

6,97(-7) 

5.60(-7) 

6.98(-7) 

7.01(-7) 

14.54 

6.85(  -7) 

6.78(-7) 

5.41(-7) 

6.78(-7) 

6.81(  7) 

13.94 

6.52(  7) 

6.46(  7) 

5.17(-7) 

6.46(-7) 

6.48(-7) 

13.48 

6.14(-7) 

6.08(-7) 

4.94(-7) 

6.10(-7) 

6.11(-7) 

13.03 

5.96(-7) 

5,90(-7) 

4.81(-7) 

5.92(-7) 

5.93(-7) 

12.43 

5.96(  7) 

5.90(-7) 

4.73(-7) 

S.91(-7) 

5.92(-7) 

11.98 

6.00(-7) 

5.94(-7) 

4.72(-7) 

5.94(-7) 

5.96(-7) 

11.52 

5.96(  7) 

590f-7) 

4.63(-7) 

5.88(-7) 

5.91(-7) 

11.07 

5.85(-7) 

5.79(  7) 

4.55(-7) 

5.78(-7) 

5.80(-7) 

10.47 

5.61(  7) 

5,54(-7) 

4.34(-7) 

5.53(-7) 

5.55(-7) 

10.00 

5.53(  7) 

5.46(-7) 

4.25(-7) 

5.44(-7) 

5.47(-7) 

9.53 

5.38(  7) 

5.31(-7) 

4.15(-7) 

5,30(-7) 

5.32(-7) 

9.00 

5.21(-7) 

5.14(-7) 

4.01(-7) 

5.13(-7) 

5.15(-7) 

8.55 

,  .  5.05(  7) 

4.99(  7) 

3.90(-7) 

4.98(-7) 

5.00(  7) 

8.02 

'  4.96(  7) 

4.90(-7) 

3.91(-7) 

4.92(-7) 

4,92(-7) 

7.49 

4.94(-7) 

4.88(-7) 

3.85(-7) 

4.89(-7) 

4.98(-7) 

6.97 

4.73(-7) 

4.67(-7) 

3.60(-7) 

4.66(-7) 

4.67(-7) 

6.52 

4.67(-7) 

4.61(-7) 

3.S7(-7) 

4.60(-7) 

4.61{-7) 

6.00 

4.58(-7) 

4.53(~7) 

3.52(-7) 

4.52(-7) 

4.53(-7) 

5.54 

4.41(  7) 

4.36(-7) 

3.40(-7) 

4.36(-7) 

4.37(-7) 

5.00 

4.39(-7) 

4.34(-7) 

3.36(-7) 

4.33(-7) 

4.35(-7) 

4,50 

4.13(-7) 

4.08(-7) 

3.20(-7) 

4.08(-7) 

4.09(-7) 

4.00 

4.09(-7) 

4.04(-7) 

3.19(-7) 

4.05(-7) 

4.06(-7) 

3.52 

4.05(-7) 

4.00(-7) 

3.15(-7) 

4.01(-7) 

4.02(-7) 

3.00 

3.59(-7) 

3.54{-7) 

2.77(-7) 

3.55(-7) 

3.5S(-7) 

2.50 

3.28(  7) 

3.24(-7) 

2.51(-7) 

3.25(-7) 

3.2S(-7) 

2.00 

3.01(-7) 

2.97(  7) 

2.29(-7) 

2.98(-7) 

2.98(-7) 

1.80 

2.90(-7) 

2.87(-7) 

2.20(-7) 

2.87(-7) 

2.87{-7) 

1.60 

2.74(  7) 

2.71{-7) 

2.08(-7) 

2.71(-7) 

2.71(-7) 

1.40 

2.59(  7) 

2.56(-7) 

1.97(-7) 

2.56(-7) 

2.57(-7) 

1.20 

2.45(-7) 

2,42(-7) 

1.84(-7) 

2.42(-7) 

2.42(-7) 

1.00 

2.51(-7) 

2.4  7( -7) 

1.83(-7) 

2.46(-7) 

2.47(-7) 

8.05(-l) 

2,00(  7) 

1,97(  7) 

l,Sl(-7) 

l.98(-7) 

1.98(-7) 

6.0(  1) 

1.73(  7) 

I.71(-7) 

1.31(-7) 

1.7l(-7) 

1.71(-7) 

5.0(-l) 

l.60(-7) 

1.58(  7) 

1.21(-7) 

1.58(-7) 

1.58(-7) 

4,0(  1) 

l,51(-7) 

1.48(-7) 

l.ll{-7) 

1.48(-7) 

1.50(-7) 

3.0(-l) 

1.23(  7) 

1.22(-7) 

9.27(-8) 

1.22(-7) 

1.22(-7) 

2.0(  1) 

9.73(-8) 

9,67(  8) 

7.33(-8) 

9,61(-8) 

9.63(-8) 

1.0(  1) 

6.39(-8) 

6.30(-8) 

4.80(-8) 

6.31(-8) 

6.31(-8) 

5.03(  2) 

3.88(-8) 

3.83(-8) 

2.92(-8) 

3.84(-8) 

3.84(-8) 

2.04(  2) 

l.82(-8) 

l,80(-8) 

1.38(-8) 

1.80(  8) 

l.80(  8) 

1.02(  2) 

9.69(  9) 

9.57(-9) 

7.54(-9) 

9.57(-9) 

9.57(-9) 

4.97(-3) 

4.90(-9) 

4.84(  9) 

4,02(-9) 

4.84(-9) 

4,84(-9) 

3.00(  3) 

2,99(-9) 

2,95(-9) 

2.61(-9) 

2,96(-9) 

2.94(-9) 

2.014(  3) 

2.02(  9) 

2.00(-9) 

l,90(-9) 

2,00(-9) 

1.99(-9) 

1.025(  3) 

1.04(  9) 

1.03{-9) 

l,18{-9) 

1.03(-9) 

1.025(-9) 

4.94(  4) 

5.16(  10) 

5,11(  10) 

7.94(-10) 

5.12(-10) 

5.08(-:10) 

2.01(-4) 

2.30(-10) 

2.29(-IO) 

5.94(-10) 

2.27(  10) 

2.29(-10) 

1.00(  4) 

l.39(-10) 

1.39(-10) 

5.38(-10) 

1.36(-10) 

1.41(-10) 

5.03(  5) 

1.03(  10) 

l.04(  10) 

5.26(  10) 

9.59(-ll) 

1.09(  10) 

3.43(  5) 

9.72(  11) 

9.78(-l)) 

5.32(-10) 

8.76(  11) 

I.OS(-IO) 

2.05(  5) 

I.OK  10) 

l.Ol(-lO) 

5.52(-10) 

8.65(  11) 

1,12(-10) 

1.03(  -5) 

l.22(  10) 

1.22(  10) 

6.00(  10) 

9.93(  11) 

1.40(-10) 

5.I4(^6) 

1.61(  -10) 

1.61(-10) 

6.75(  10) 

1.27(  10) 

1.88(  10) 

2,00(  6) 

2.49(  10) 

2,47(  10) 

8.33(-10) 

1.90(  10) 

2.93(  10) 

l.OOt  -6) 

3.48(-10) 

3.45(-10) 

1.01(-9) 

2.62(-10) 

4.12(  10) 

5.0(-7) 

4.89(  10) 

4,83(10) 

1.26(-9) 

3.6S(  10) 

5.81(-10) 

2.5{  7) 

6.87(-10) 

6.79(  10) 

1.61(  -9) 

5.09(  10) 

8.18(-10) 

1.07(  7) 

l.04(  9) 

l.03(  -9) 

2.24{-9) 

7.68(-10) 

1.25(-9) 

4.5(  8) 

1.61(  9) 

1.59(-9) 

3.24(  9) 

1.18(-9) 

1.93(-9) 

2.23(  8) 

2.27(  9) 

2.24(  9) 

4.4  2( -9) 

1.66(-9) 

2.72(-9) 

"From  ref.  43. 

'’Read:  7,04  X  10  ’.etc. 
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6.97  7.00(-6)  1.39(-6)  9.13(-7)  2.8(K-6)  3.82(  7)  2.93(-6)  4.46(  7)  6.54(-6)  7.57(-7)  2.16(-6)  1.31(-6) 

6.82  6.98(-6)  1.83(-6)  1.00(-6)  2.62(  6)  3.78(-7)  2.9U-6)  4.43{  7)  6.37(-6)  7.61(  7)  2.11(-6)  1.24(-6) 

6.67  6.96(-6)  1.74(-6)  9.76(-7)  2.62(-6)  3.74(-7)  2.77(-6)  4.39(-7)  6.21(-6)  7.69(  6)  2.07(-6)  1.07(-6) 

6.52  6.93(-6)  1.31(-6)  1.04(  6)  2.71(-6)  3.70(-7)  2.65(-2)  437(  1)  6.06(  6)  7.81(-6)  2.02(-6)  9.97(-7) 

6.37  6.91(-6)  1.72(-6)  1.42(-6)  2.86(-6)  3.66(-7)  2.58(-6)  4.34(-7)  5.90(  6)  8.06(-6)  1.98(  6)  8.69(-7) 
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Appendix  2B.  Spatial  Distributions  of  Absorbed  Dose  and 
Dose  Equivalent  in  Cylindrical  Phantom  Due  to 
Incident  Monoenergetic  Neutrons 


This  appendix  presents  calculated  absorbed  doses 
and  dose  equivalents  in  a  four-element  cylindrical  tissue 
phantom  due  to  a  broad  beam  of  0.025-eV  to  14-MeV 
monoenergetic  neutrons  incident  on  the  phantom.  The 
calculations,  performed  by  Auxier  et  al.^'  for  the 
phantom  schematic  shown  in  Fig.  2.20,  considered 
dose  contributions  by  both  the  incident  neutrons  and 
the  capture  gamma  rays  they  produced.  The  neutron 
beam  was  assumed  to  be  incident  normal  to  the  axis  of 


the  cylinder,  and  the  absorbed  doses  and  dose  equiva¬ 
lents  given  in  Table  2B.1  each  represent  the  average 
value  over  a  radial  volume  element  in  layer  3,  which 
is  locatid  at  the  midplane  of  the  cylinder.  These  values 
are  multicollision  doses  per  unit  fluence  and  thus  are 
fluence-to-dose  conversion  factors  for  cylindrical  phan¬ 
toms. 

Table  2B.2  gives  the  absorbed  doses  due  to  the  sec¬ 
ondary  (capture)  gamma  rays  only. 
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